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Abstract 

Photonic  crystals  (PhC’s)  are  periodic  structures  of  differing  dielectrics  that 
create  a  photonic  band  gap  (PBG).  A  PBG,  in  turn,  inhibits  the  propagation  of  elec¬ 
tromagnetic  waves  of  a  specific  frequency  range.  This  thesis  focuses  on  the  fabrication 
and  characterization  of  triangular-structured,  two-dimensional  PhC’s  with  a  PBG  de¬ 
signed  for  visible  wavelengths  and  with  applications  in  visible  integrated  photonic 
systems.  A  three-dimensional  PhC  with  a  PBG  in  the  infrared  is  also  studied  for  its 
characteristics  in  regard  to  its  PBG. 

The  two-dimensional  fabrication  processes  pursued  were:  focused  ion  beam, 
electron  beam  lithography  and  holographic  photo-polymerization/lithography.  The 
fabrication  techniques  and  materials  used  to  create  the  PhC  in  part  determined  the 
characterization  technique  required  to  investigate  the  PBG.  Characterization  tech¬ 
niques  include:  the  coupling  of  a  beam  by  means  of  a  prism  into  a  wave-guiding 
medium  in  which  the  PhC  has  been  fabricated,  Fourier  transform  infrared  spectrom¬ 
eter,  spectrophotometer,  and  edge  firing  techniques.  Analysis  of  the  transmission  and 
reflectance  properties  of  a  PhC  for  various  incident  angles  (within  the  two  dimen¬ 
sional  plane  of  the  PhC)  confirms  the  presence  of  a  PBG.  The  design  of  the  PhC  was 
based  on  a  program  created  to  display  gap  maps  for  triangular  structures.  The  PBG 
of  the  structure,  designed  from  the  gap  maps,  was  modeled  using  another  program 
created  to  display  band  diagrams  for  triangular  structures.  Finite  difference  time 
domain  (FDTD)  software  was  then  applied  to  compare  the  experimental  results  with 
the  modeling  and  simulation  of  the  PhC  fabricated.  Two  different  FDTD  softwares 
were  also  evaluated. 

The  three-dimensional  PhC  sample  obtained  was  successfully  tested  with  the 
Fourier  transform  infrared  spectrometer  and  the  spectrophotometer.  The  results  of 
fabricating  two-dimensional  PhC’s,  on  the  order  of  visible  wavelengths,  is  shown  for 


IV 


the  holographic  photo-polymerization/lithography  process  and  the  focused  ion  beam. 
Modeling  and  simulation  of  the  these  designs  were  accomplished  and  compared  to 
some  of  the  experimental  results.  A  design  that  incorporates  the  electron  beam  lithog¬ 
raphy  process  was  created  and  modeled. 
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CHARACTERISTICS  OF  TWO-DIMENSIONAL 


TRIANGULAR  AND  THREE-DIMENSIONAL 
FACE-CENTERED-CUBIC  PHOTONIC  CRYSTALS 

I.  Introduction 

Photonics  has  come  to  be  known  as  the  relationship  between  electronics  and 
optics,  optics  being  grounded  in  the  “generation,  propagation  and  detection  of 
light”  [18].  Thus,  photonics  is  the  control  of  photons  [18].  A  photonic  crystal  (PhC) 
is  one  way  to  control  these  photons.  More  importantly,  a  PhC  with  a  photonic  band 
gap  (PBG)  prevents  the  propagation  of  certain  wavelengths  within  the  PhC.  This 
thesis  focuses  on  the  capabilities  of  fabrication  and  characterization  of  triangular- 
structured,  two-dimensional  PhC’s  with  a  PBG  designed  for  visible  wavelengths  and 
with  applications  in  visible  integrated  photonic  systems.  A  three-dimensional  PhC 
with  a  PBG  in  the  infrared  is  also  studied  for  its  characteristics  in  regard  to  its  PBG. 
The  three  different  techniques  for  fabricating  these  PhC’s  are:  focused  ion  beam, 
electron  beam  lithography  and  holographic  photo-polymerization/lithography.  An 
evaluation  of  two  different  finite  difference  time  domain  (FDTD)  software  programs 
will  be  discussed. 

PhC’s  can  be  of  one-,  two-  or  three-dimensional  structures.  A  quarter-wave 
stack  is  an  example  of  a  one-dimensional  PhC;  examples  of  two-  and  three-dimensional 
PhC’s  will  be  covered  in  this  work.  The  history  and  theory  of  PhC’s  is  essential  to 
the  understanding  of  the  fabrication  techniques  used  and  is  presented  in  this  chapter. 

1.1  History  Behind  Photonic  Crystals 

There  are  many  examples  throughout  history  of  man  manipulating  and  using 
the  observable  properties  of  optics.  It  is  documented  that  Euclid,  in  300  B.C.,  knew 
of  the  Law  of  Reflection,  and  Plato,  in  50  A.D,  documented  the  bending  of  light  in 
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water.  Evidence  shows  that  the  Romans  possessed  burning  glasses,  thought  to  be 
similar  to  magnifying  glasses  [9].  By  the  end  of  the  thirteenth  century,  Bacan  (1215- 
1294)  began  inventing  lenses  for  the  use  of  correcting  vision  [9] . 

However,  in  the  late  1800’s,  discoveries  made  by  a  few  people  in  a  short  period 
of  time  created  a  giant  leap  forward  in  understanding.  Planck,  in  1899,  published  his 
theory  on  the  quantization  of  light,  which  directly  affected  the  development  of  Ein¬ 
stein’s  photoelectric  effect  in  1905.  Concurrently,  Lord  Rayleigh  began  studying  wave 
propagation  in  a  periodic  media  (1887)  [14],  later  progressing  to  more  direct  studies  of 
a  one-dimensional  PBG,  in  which  he  showed  systems  that  produced  such  effects  [14]. 
Later  in  1904,  he  discovered  the  phenomenon  now  called  Rayleigh  scattering.  In 
1922,  the  publication  and  study  of  Bohr’s  model  of  the  atom  furthered  the  under¬ 
standing  of  the  nature  of  photons  on  the  atomic  scale.  In  1928,  Bloch  studied  wave 
propagation  in  a  three-dimensional  periodic  material  proving,  by  adding  Maxwell’s 
equations  and  Schrodinger’s  equations,  that  waves  in  such  a  media  can  propagate 
without  scattering  [14],  thus  paving  the  way  for  three-dimensional  PhC’s  and  PBG’s. 
Ohtaka,  in  1979,  derived  the  vector-spherical-wave  expansion  method,  allowing  the 
calculations  of  dispersion  and  transmittance  for  arrays  of  dielectric  spheres  [1].  But, 
the  most  directly  related  advancement  was  made  in  1987,  when  two  independent 
researchers,  Sajeev  John  and  Eh  Yablonovitch,  with  two  different  goals,  submitted 
proposals  within  months  of  each  other  on  the  probability  and  possibility  of  PBG’s  [25]. 
Since  then,  three-dimensional  structures  have  been  predicted  [12]  and  fabricated  by 
John  and  Yablonovitch  in  1987.  Since  John  and  Yablonovitch  were  the  Erst  to  fab¬ 
ricate  a  three-dimensional  PhC  with  a  PBG,  the  terms  that  they  chose  to  describe 
the  phenomena  (photonic  band  gap  and  photonic  crystal)  were  adopted,  and  they  are 
credited  with  being  the  fathers  of  this  field. 

1.2  Theory  of  Photonic  Crystals 

A  PhC  is  a  periodic  structure  of  differing  dielectrics  that  creates  a  PBG.  A 
crystal-like  material  is  one  that  has  a  periodic  arrangement  of  atoms  or  molecules  [12], 
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periodic  in  periodic  in 

one  direction  two  directions 

Figure  1.1:  Examples  of  one-,  two-  and  three- 

dimensional  crystal  structures.  The  red  and  yellow  col¬ 
ored  regions  are  regions  of  different  dielectric  constants 

[14]- 

The  geometric  placement  of  these  periodic  arrangements  is  crucial  to  the  creation  of 
PhC’s  -  photonic  because  they  scatter  light,  crystals  because  they  are  of  a  periodic 
arrangement  [16].  Thus,  the  fundamental  principals  of  a  PBG  are  the  same  as  that  of 
the  energy  band  gap  of  a  semiconductor  and  analogous  qualities  exist  between  the  two. 
The  energy  band  gap  of  a  semiconductor  is  a  complete  band  gap  of  that  semiconductor 
and  restricts  the  existence  of  electrons  of  the  energy  band  gap  states  to  exist  in  the 
material.  A  three-dimensional  PhC,  with  a  PBG,  has  the  same  characteristics  as  that 
of  the  energy  band  gap  of  a  semiconductor.  PhC’s  prevent  the  propagation  of  certain 
frequencies  (frequencies  on  the  order  of  the  periodic  arrangement);  thus,  if  such  a 
frequency  is  incident  upon  the  PhC,  it  will  be  reflected.  The  full  three-dimensional 
PBG  is  a  PBG  in  the  crystal  and  restricts  the  ability  of  any  photon  of  that  wavelength, 
propagating  at  any  angle,  to  exist  in  the  material.  Note  that  a  one- dimensional  or  two- 
dimensional  PhC  that  contains  a  PBG  does  not  have  a  complete  PBG  based  on  the 
previous  definition;  depending  on  the  angle  of  incidence,  a  different  periodic  structure 
will  be  realized  and  there  will  be  a  momentum  vector  in  the  reciprocal  lattice  space 
(Brillouin  zone  or  k-space)  in  which  the  frequency  can  exist  within  the  crystal.  Some 
basic  examples  of  one-,  two-  and  three-dimensional  periodic  structures  are  shown  in 
Figure  1.1.  Understanding  the  details  of  a  PBG  requires  solving  Maxwell’s  equations 
within  a  periodic  dielectric  material.  Maxwell’s  equations  govern  the  propagation  of 
light  [12]. 
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V  ■  B  =  0 


V  •  D  =  Anp 


(1.1) 


where  E  is  the  electric  held,  H  is  the  magnetic  held,  D  is  the  displacement  held,  B  is 
the  magnetic  induction  held,  J  is  free  current,  p  is  the  free  charge,  c  is  the  speed  of 
light  and  t  is  time. 

If  there  are  no  free  charges  and  no  free  currents  in  the  material,  then  J  =  0 
and  p  =  0.  The  material  is  then  a  linear,  lossless,  constant  dielectric  material  and 
Maxwell’s  equations  reduce  to: 


V  •  //(r,  t)  =0 


V  ■  e(r)E(r,t)  =  0 


(1.2) 


where  r  is  the  wave  vector.  With  a  few  more  algebraic  manipulations,  an  equation  in 
regards  to  the  magnetic  held  ( H(r ))  [12]  is  observed. 


v  x  ^ v  x  H{r)]  =  02R(r) 


(1.3) 


Solving  for  H( r),  the  electric  held,  E( r)  [12],  is  found  by: 


—ic 


)V  x  H (r) 


(1.4) 
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By  rewriting  equation  (1.3)  in  the  following  form,  it  can  be  seen  that  an  eigenvalue 
problem  exists  [12]. 

QH(  r)  =  (~)2H(  r)  (1.5) 

c 

where  0  takes  the  curl,  divides  by  e(r)  and  takes  the  curl  again  [12].  This  operator, 
0,  is  a  linear  operator,  and  thus,  in  defining  the  modes,  since  0  is  linear,  any  linear 
combination  of  solutions  is  also  a  solution.  So,  any  fields  that  differ  only  by  a  linear 
multiplier  are  of  the  same  mode  [12].  In  order  to  solve  for  the  specific  frequencies 
that  are  either  allowed  or  forbidden  to  propagate  through  the  crystal,  the  previously 
derived  equations  will  be  used.  If  the  fields  are  represented  as  a  set  of  plane  waves  and 
require  that  H{ r)  has  the  same  periodicity  of  e(r),  by  using  the  crystal’s  reciprocal 
lattice  vectors,  then  each  mode  can  be  represented  with  a  wave  vector,  k.  Also,  by 
expanding  e(r)  as  plane  waves,  these  equations  of  plane  waves  can  be  inserted  into  the 
original  equation  (1.3).  Solving  this  Hcrmitian  eigenvalue  problem  [12]  is  the  basic 
equation  behind  photonic  crystals. 

There  are  several  computer  software  programs  available  that  make  these  com¬ 
putations;  however,  each  one  can  be  slightly  different  in  the  calculations  based  on 
assumptions  made  to  solve  Maxwell’s  equations.  Regardless,  solving  these  equations 
by  hand  would  be  a  long  and  tedious  endeavor  at  best  and  in  most  cases  impossible 
since  most  cases  can  only  be  solved  numerically. 

Since  these  solutions  are  periodic  with  regards  to  the  wave  vector,  k,  they  will 
be  analyzed  with  respect  to  the  Brillouin  zone.  The  PBG  is  obtained  due  to  the  fact 
that  the  fields  of  the  lowest  mode  try  to  satisfy  the  condition  of  obtaining  low  energy 
by  having  a  low  curl  (slow  oscillations).  Also,  this  mode  desires  to  have  the  electric 
field  concentrated  in  the  high  dielectric  region.  Since  the  next  highest  mode  has  to 
be  orthogonal  to  the  lower  mode,  it  will  concentrate  its  energy  in  the  low  dielectric. 
This  gives  rise  to  a  frequency  band  gap  [16]. 

These  two  modes  are  also  analogous  to  the  semiconductor  band  gap  with  re¬ 
spect  to  electrons;  electrons  in  the  conduction  band  (like  modes  concentrated  in  low 
dielectric)  feel  fewer  restrictions  to  them  than  those  in  the  valence  band  (like  modes 
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Figure  1.2:  One-dimensional  crystal,  the  blue  and 

green  represent  the  different  materials  of  different  dielec¬ 
tric  constants  [12]. 


concentrated  in  the  high  dielectric  region).  With  these  concepts  in  mind,  a  one¬ 
dimensional  crystal  with  a  PBG  can  be  examined.  A  quarter-wave  dielectric  stack  is 
a  good  example  of  a  one-dimensional  PhC  containing  a  PBG  (Figure  1.2).  As  long 
as  the  incident  wave  is  perpendicular  to  the  layers  (and  has  a  wavelength  that  is  four 
times  the  optical  thickness  of  each  layer),  it  will  be  totally  reflected,  due  to  Bragg 
reflection.  However,  if  the  incident  wave  is  no  longer  perpendicular  to  the  layers,  it 
will  not  be  perfectly  reflected.  Looking  at  the  Brillouin  zone  for  this  case,  a  PBG  at 
the  edges  (k  =  it /a,  where  ‘a’  is  the  lattice  constant)  is  observed  [12];  the  modes  at 
this  point  are  standing  waves,  i.e.,  the  periodic  structure  is  180°  out  of  phase  with 
the  others  in  the  crystal  at  this  frequency.  The  modes  will  tend  to  be  located  in  high 
dielectric  material  for  the  low  frequency  mode  and  the  low  dielectric  material  for  the 
next  higher  frequency  mode  (Figure  1.3)  [12].  In  a  one- dimensional  photonic  crystal, 
there  will  always  be  a  PBG  in  the  direction  of  propagation  (kz)  for  some  frequency, 
but  not  in  the  off  axis  {e.g.  ky)  direction  for  that  same  frequency. 

For  a  two-dimensional  structure  (Figure  1.4),  PBG’s  occur  at  some  frequency 
for  any  direction  within  the  two-dimensional  plane.  Consider  the  two-dimensional 
plane  to  be  the  x-y  plane,  and  therefore,  the  desire  is  to  create  a  PBG  for  any  di¬ 
rection  in  that  plane.  This  can  be  accomplished  with  a  specific  periodic  structure 
(correlating  to  the  wavelength  of  choice)  and  differing  dielectric  constants  in  the  x-y 
plane.  However,  now  there  will  be  a  Brillouin  zone  that  has  two  distinct  polarizations 
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(a)  E-field  for  mode  at  top  of  band  1 


Figure  1.3:  E-field  distribution  on  high  dielectric  (blue) 
and  low  dielectric  (green)  [12]. 

to  it  that  is  not  typically  shown  in  the  one-dimensional  case,  since  the  one- dimensional 
case  is  usually  dealing  with  the  incident  wave  perpendicular  to  the  surface  (Figure 
1.5)  [12].  These  polarizations  are  the  transverse  electric  (TE)  and  transverse  magnetic 
(TM)  waves.  TE  waves  are  defined  as  the  case  where  the  electric  component  is  in  the 
plane  of  incidence  and  the  magnetic  component  is  normal  to  the  plane  of  incidence. 
The  TM  case  is  where  the  magnetic  component  is  in  the  plane  of  incidence  and  the 
electric  component  is  normal  to  the  plane  of  incidence.  Figure  1.5  is  a  graph  of  the 
Brillouin  zone  where  the  TM  modes  have  a  PBG  but  the  TE  modes  do  not.  Thus,  the 
PBG’s  differ  for  different  polarizations  [12].  After  several  experimentations  and  cal¬ 
culations,  it  was  found  that  isolated  high  dielectric  regions  make  good  TM  PBG’s  and 
connected  lattices  make  good  TE  PBG’s  [12].  Therefore,  a  compromise  is  necessary 
to  get  a  complete  PBG  for  both. 

Three-dimensional  lattice  structures  are  disco uragingly  difficult  to  construct. 
There  are  a  few  basic  designs  that  can  be  produced  accurately;  however,  one  of  the 
challenging  tasks  is  to  construct  a  three-dimensional  structure  on  the  scale  necessary 
to  be  used  in  most  applications.  Two  common  types  are  shown  in  Figure  1.6.  The  two- 
dimensional  examples  provide  building  blocks  to  use  for  assembling  three-dimensional 
structures.  Figure  1.6  illustrates  the  use  of  spheres  placed  in  a  periodic  sequence  in 
which  the  spheres  are  touching,  and  then  arranged  in  a  diamond  lattice  structure  to 
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Figure  1.4:  Illustration  of  a  two-dimensional  PhC  of 

columns  of  dielectric  in  air,  where  ’a’  is  the  lattice  con¬ 
stant  or  periodicity  and  V  is  the  radius  of  the  columns 
[12]- 


Figure  1.5:  Two-dimensional  crystal  Brillouin  zone  representation  of 
frequency  versus  k-space.  T  corresponds  to  k  =  0,  x  corresponds  to 
k  =  (^)x,  M  corresponds  to  k  =  (^)x  +  (-)y,  The  inset  on  the  left 
shows  the  wave  vector  with  respect  to  the  reciprocal  lattice  space.  The 
inset  on  the  right  shows  the  reciprocal  lattice  of  the  crystal  [12]. 
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(a) 


Figure  1.6:  Two  examples  of  three-dimensional  pho¬ 

tonic  crystal  structures,  (a)  Spheres  in  a  diamond  lattice 
structure  and  (b)  the  same  structure  where  the  spheres 
were  replaced  with  rods  [12]. 
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Figure  1.7:  Brillouin  zone  band  diagram  of  a  three- 

dimensional  crystal  of  frequency  versus  k-space  of  incor¬ 
rect  geometry  and  large  dielectric  contrast  of  13  to  1. 

provide  the  same  periodic  look  in  all  three  directions  [12].  This  figure  also  shows  taking 
a  diamond  lattice  structure  and  connecting  the  points  with  columns  of  dielectric  [12], 
For  clarification,  it  should  be  noted  that  the  larger  the  dielectric  contrast  between  the 
two  materials  that  make  up  the  PhC,  the  larger  the  PBG.  However,  this  is  not  the  only 
requirement  for  a  PBG  in  a  one-,  two-  or  three-dimensional  crystal.  Figure  1.7  depicts 
a  Brillouin  zone  for  a  three-dimensional  crystal  with  a  dielectric  contrast  of  13  to  1; 
yet,  there  is  no  band  gap  due  to  the  incorrect  geometry  of  the  crystal  [12],  However, 
Figure  1.8  shows  a  three-dimensional  crystal  structure  arranged  in  interpenetrating 
diamond  lattice  spheres  that  have  produced  a  band  gap  [12],  Many  other  two- 
and  three-dimensional  crystal  structures  have  and  are  currently  being  designed  and 
fabricated. 

1.3  Conclusion 

It  was  shown  in  this  chapter  how  the  culmination  of  works  from  several  great 
men  in  history  lead  to  the  eventual  understanding  and  creation  of  PhC’s  as  discovered 
by  Sajeev  John  and  Eli  Yablonovitch,  the  fathers  of  PhC’s.  A  brief  introduction  of  the 
theory,  explaining  PhC’s,  was  presented,  illustrating  one  of  the  most  important  facts 
of  PhC’s;  the  requirement  for  the  proper  dielectric  contrast  and  periodic  geometry 
of  the  PhC  is  detrimental  in  obtaining  a  PBG.  In  the  next  chapter,  an  overview  of 
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Figure  1.8:  Three-dimensional  photonic  crystal  Bril- 

louin  zone  representation  of  frequency  versus  k-space  of 
a  complete  PBG  for  a  interpenetrating  diamond  lattice 
of  spheres  [12], 


the  different  types  and  fabrication  techniques  of  two-  and  three-dimensional  PhC’s  is 
covered,  illustrating  the  variety  of  fabrication  techniques  possible  depending  on  the 
wavelength  band  of  the  PBG  and  the  application  required  of  the  PhC’s. 
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II.  Background  of  Two-  And  Three-Dimensional  Photonic 

Crystal 


lie  utility  of  PhC’s  is  becoming  increasingly  popular  in  all  areas  of  engineering; 


_L  some  applications  include  telecommunications,  lasers,  electro-optical  comput¬ 
ing,  sensors  and  wave  guides.  As  the  process  to  fabricate  two-dimensional  PhC’s  with 
complete  in-plane  PBG’s  becomes  easier  and  more  industrially  available,  this  need 
will  increase  exponentially.  Along  with  this  advancement  in  fabrication  will  arise  the 
advancement  of  applications  and  designs  of  two-  and  three-dimensional  PhC’s.  Ulti¬ 
mately,  applications  involving  a  three-dimensional  PhC’s  with  complete  PBG’s  will 
be  inevitable  and  required  for  technological  advancement  in  this  field.  In  order  to 
reach  these  goals,  a  building  block  approach  must  be  accomplished.  Two-dimensional 
PhC’s  must  be  able  to  be  fabricated  and  characterized  for  a  broader  range  of  frequen¬ 
cies  in  order  to  heighten  the  knowledge  in  this  field,  thus,  paving  the  way  for  better 
three-dimensional  PhC’s. 

In  Chapter  I,  the  history  and  theory  of  PhC’s  was  discussed,  where  it  was  shown 
that  a  complete  PBG  can  be  created  for  a  specific  frequency  range  depending  on  the 
geometry  and  dielectric  contrast  of  the  features  in  the  PhC.  This  chapter  focuses  on 
the  considerations  of  fabrication  and  testing  of  PhC’s,  as  well  as  introducing  some 
recent  work  in  two-  and  three-dimensional  PhC’s. 

2.1  Fabrication  and  Characterization  of  Photonic  Crystals 

The  visible  frequency  regime  has  had  very  little  successful  experimentation  con¬ 
ducted  and  needs  more  exploration  to  further  PhC  research  and  design  in  that  area. 
Fabrication  techniques  for  two-dimensional  PhC’s  vary  depending  on  the  size  and 
type  of  PhC  desired.  Attempts  to  fabricate  PhC’s  with  complete  two-  and  three- 
dimensional  PBG’s  have  been  met  with  some  success.  This  thesis  examines  the  ca¬ 
pabilities  to  fabricate  two-dimensional  PhC’s  with  complete  PBG’s  designed  for  the 
visible  range  and  to  characterize  two  and  three-dimensional  PhC’s.  Different  analy¬ 
sis  processes  can  be  used  to  determine  the  presence  of  a  two-  and  three-dimensional 
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PBG  within  a  PhC,  including:  measuring  the  reflection  amplitude  of  the  incident 
beam  on  the  PhC,  measuring  the  transmittance,  recording  the  specular  scattering 
that  occurs  incident  on  the  PhC,  and  obtaining  the  near  field  by  using  scanning  opti¬ 
cal  microscopy,  to  name  a  few.  For  the  fabrication  techniques  used  in  this  experiment, 
characterization  could  be  accomplished  by  measuring  the  transmittance  and  observ¬ 
ing  the  specular  scattering  produced  by  beams  of  different  wavelengths  and  different 
angles  incident  on  the  PhC’s.  In  order  to  better  understand  the  results  obtained, 
each  device  must  be  modeled  and  simulated  under  similar  conditions.  The  use  of  a 
commercial  FDTD  software  will  be  used  to  model  and  simulate  the  designs  of  this 
thesis. 

One  program  that  is  used  for  this  experiment  is  a  two-dimensional  triangular 
structured  PhC  band  solver  created  by  Dr.  Joseph  Haus  of  the  University  of  Day- 
ton  [8].  This  program  is  listed  in  Appendix  C.  Appendix  C  displays  another  program 
that  creates  gap  maps  based  on  the  materials’  dielectric  constants,  the  periodicity  of 
the  triangular  columns  of  air  and  a  specified  range  of  the  air  column  radii.  Commer¬ 
cial  FDTD  software  and  band-solver  software  was  then  applied  to  model  and  simulate 
the  initial  design.  The  modeled  results  would  then  be  compared  to  the  actual  test 
data.  There  are  several  companies  that  create  software  to  evaluate  these  types  of 
designs.  Characteristics  and  examples  of  two  different  commercial  softwares  will  be 
presented  as  part  of  this  study. 

2.2  Recent  Studies  of  Two-Dimensional  PhC’s 

Two-dimensional  PhC’s  have  been  successfully  fabricated  over  the  past  five  years 
and  can  be  produced  (depending  on  wavelength  of  the  PBG)  without  too  much  diffi¬ 
culty.  Some  of  the  latest  work  in  two-dimensional  PhC’s  is  in  the  realm  of  tuneable 
PBG’s.  Tuning,  or  changing,  the  band  gap  in  real  time  can  be  accomplished  by  several 
methods.  One  such  method  would  be  to  change  the  dielectric  constant,  and  another 
would  be  to  change  the  periodic  structure.  Henry  Driel  et  al.  have  investigated  three 
different  ways  of  creating  a  tunable  two-dimensional  PhC  [7].  One  method  was  to 
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fill  the  periodic  columns  of  air  in  the  silicon  substrate  with  a  nematic  liquid  crystal, 
and  by  adjusting  the  temperature,  they  were  able  to  shift  the  PBG.  Another  method 
was  to  optically  inject  electron-hole  pairs  into  the  silicon  and  change  the  dielectric 
constant;  approximately  30-nm  shifts  were  recorded.  As  can  be  seen  by  the  scanning 
electron  microscope  (SEM)  image  in  Figure  2.1,  their  periodicity  is  on  the  order  of 
1.5  /jm  which,  in  turn,  produced  a  PBG  around  3-6  /irn  as  illustrated  in  Figure  2.2. 
Lastly,  by  using  non-resonant,  nonlinear  induced  optical  changes  (Kerr  effect),  they 
were  able  to  achieve  faster  tuning  effects.  By  utilizing  a  pre-patterned  electrochemical 
etching  process,  they  fabricated  a  triangular  lattice  of  air  pores  of  high  geometrical 
accuracy.  Upon  applying  a  temperature  change  to  the  substrate,  a  shift  in  the  band 
gap  was  noticed;  however,  it  was  not  as  much  as  predicted  and  only  affected  the  air 
band  gap  edge  (higher  mode)  with  any  appreciability.  This  can  be  explained  due  to 
the  fact  that  it  is  the  air  band  edge  that  will  change  the  most  since  it  housed  the 
medium  that  was  changing  [7] . 

Christopher  Summers,  Georgia  Institute  of  Technology,  et  al.  merged  the  abil¬ 
ity  to  tune  both  the  refractive  and  dispersion  properties  [4].  “Super  lattice”  struc¬ 
tures,  as  they  call  them,  were  the  result  of  their  work.  They  took  an  ordinary  two- 
dimensional  column  structure  aligned  in  a  triangular  shape  and  made  the  alternating 
rows  of  a  different  dielectric  column  and  resized  it  to  be  smaller  in  diameter  (see  Fig¬ 
ure  2.3).  This  had  the  effect  of  folding  the  Brillouin  zone  into  a  rectangular  zone  when 
liquid  crystal  (LC)  is  deposited  within  the  lattice  structure  and  by  addressing  alter¬ 
nate  rows  (changing  the  refractive  indices).  Due  to  this  unique  property,  depending 
on  the  direction  of  propagation  through  the  crystal,  it  can  be  switched  from  having 
no  PBG  to  having  a  PBG,  thus  either  allowing  light  to  flow  through  or  not  [4], 

Schonbrun  et  al.  have  demonstrated  a  technique  to  change  the  dielectric  of  the 
lattice  by  stressing  the  material.  A  dual  set  of  electrostatic  comb  actuators  were  used 
to  effectively  stress  the  material  as  depicted  in  Figure  2.4  [24],  What  they  ended  up 
proving  from  their  testing  was  an  ability  to  change  the  refraction  angle  by  up  to  75°. 
Their  structure  consisted  of  an  array  of  high  index  dielectric  embedded  in  a  flexi- 
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Figure  2.1:  SEM  images  of  a  triangular  structured  two-dimensional 
PhC  with  a  periodicity  of  1.5  /irn  of  air  columns  in  silicon  [7]. 
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Figure  2.2:  (Top)  Transmission  of  H-polarized  spec¬ 

trum  prior  to  infusion  of  liquid  crystal  into  the  PhC  of 
Figure  2.1  and  (bottom)  after  infusion  of  liquid  crystal 
into  the  PhC  [7]. 
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Figure  2.3:  Super  lattice  structures,  (a)  Triangular  lat¬ 
tice,  (b)  Triangular  lattice  with  rows  of  different  dielec¬ 
tric  columns  (c)  Static  Superlattice  -  alternating  column 
sizes,  (d)  Hybrid  Superlattice  -  combination  of  (b)  and 
(c).  This  picture  shows  the  progression  of  standard  tri¬ 
angular  two-dimensional  PhC  to  the  new  Hybrid  Super¬ 
lattice  two-dimensional  triangular  PhC  [4]. 


Figure  2.4:  MEMS  actuators  used  to  stress  the  poly¬ 
mer  PhC  [24], 


blc  polymer  film  (poly-dimethylsiloxane)  and  polymide.  The  microelectromechanical 
systems  (MEMS)  electrostatic  actuator,  stresses  the  material  when  activated,  causing 
the  crystal  structure  to  change,  thus  changing  the  refractive  index  and  the  shifting 
the  refractive  angle.  They  were  able  to  achieve  up  to  75°  for  incident  angles  between 
5°  and  15°.  A  combination  of  low  pressure  chemical  vapor  deposition  (LPCVD)  and 
reactive  ion  etch  (RIE)  was  used  to  create  this  structure  [24], 

The  use  of  specific  semiconductor  materials  for  the  near  infrared  wavelengths 
was  reported  by  S.  Noda  [15].  These  are  two-dimensional  structures  that  combine 
to  make  three-dimensional  crystals.  He  reports  that  utilizing  either  gallium  arsenide 
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Figure  2.5:  Log-Log  stack,  a  two-dimensional  structure 
used  to  create  a  three-dimensional  structure  [15]. 


(GaAs)  or  indium  phosphide  (InP)  will  give  a  band  gap  in  the  low  infrared  range. 
However,  upon  review  of  the  process,  it  would  be  impractical  for  mass  production.  He 
utilized  a  GaAs  substrate  and  then  grew  an  etch  stop  layer.  On  top  of  this  etch  stop 
layer,  he  patterned  strips  using  either  electron  beam  lithography  or  photolithography. 
The  strip’s  dimensions  and  the  thicknesses  of  the  substrate  are  determined  based  upon 
the  wavelength  of  the  desired  PBG.  The  next  step  is  to  create  the  previous  steps  again 
separately,  then  bond  the  two  strips  together,  etch  away  the  top  substrate  layers  and 
repeat  as  many  times  as  desired  (see  Figure  2.5).  Another  drawback  to  this  technique 
is  the  precision  and  difficulty  of  aligning  the  first  and  third  layers.  The  resultant 
band  gap  is  reported  to  be  in  the  5-10-/im  wavelength  range.  He  reports  the  ability 
to  make  the  gap  at  any  near  infrared  frequency,  depending  on  the  geometry  of  the 
structure  as  stated  earlier.  He  was  also  able  to  create  a  band  gap  of  wavelengths  from 
1-2  /.un  as  well.  The  attenuation  was  only  about  18  dB.  To  test  his  crystal,  he  utilized 
a  globar  lamp  monochromator  and  mercury  cadmium  telluride  (HgCdTe)  detector  to 
measure  the  reflections  of  the  normally  incident  wave  in  the  <001>  direction  (see 
Figure  2.6)  [15]. 


2.3  Recent  Studies  of  Three-Dimensional  PhC’s 

As  of  late,  three-dimensional  PhC’s  are  of  the  most  desired,  but  also  the  most 
elusive  to  obtain.  Many  different  designs  have  been  introduced:  log  design,  square  pile 
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Figure  2.6:  Top  view  of  the  Log-Log  stack  showing  the 
crystal  orientation  and  periodicity  of  4  /ini  [15]. 


structure,  and  face  centered  cubic  (fee)  structures  [21];  however,  the  search  to  discover 
better  and  more  applicable  ways  to  create  three-dimensional  PhC’s  with  complete 
PBG’s  continues.  S.  Johnson  and  J.  Joannopoulos  designed  a  structure  that  gives  a 
complete  three-dimensional  band  gap  and  is  able  to  be  fabricated  using  planer  layer 
stacks  [16].  Their  structure  (Figure  2.7)  is  an  fee  structure  with  air  cylinders  in  a 
dielectric.  This  crystal  has  a  complete  PBG  around  the  1.55-/iin  wavelength;  when 
using  silicon  and  air,  a  12:1  (dielectric  constant  of  silicon  to  the  dielectric  constant  of 
air)  ratio  which  gives  a  21%  band  gap  and  an  8%  gap,  is  obtained  at  a  ratio  of  12:2 
ratio  of  dielectric  constants  [14],  This  structure  continues  to  provide  a  gap  down  to 
a  ratio  of  4:1  [16].  The  percentage  band  gap  is  a  method  of  comparing  the  quality  of 
band  gaps  of  various  PhC’s.  The  percentage  band  gap  is  calculated  by  dividing  the 
frequency  range  of  the  band  gap  by  its  mid-gap  frequency. 

J.  Fleming  et  al.  have  demonstrated  a  MEMS  fabrication  technique  that  en¬ 
ables  a  slightly  broader  range  of  frequencies  [11].  This  structure  is  a  face-centered 
tetragonal  structure  of  ‘logs’  (see  Figure  2.8).  This  new  technique  consists  of  combin¬ 
ing  a  series  of  molding,  filling  and  chemical  mechanical  polishing  (CMP)  steps.  This 
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Figure  2.7:  Three-dimensional  PhC  design  and  associ¬ 
ated  band  diagram  [16]. 


Figure  2.8:  SEM  view  of  the  log- log  stack  [11], 
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process  enabled  a  band  gap  with  the  propagation  vector  normal  to  the  substrate  to 
appear  from  about  10-14  microns  [11].  They  accomplished  this  by  patterning  parallel 
rods  separated  by  half  the  wavelength.  The  next  layer  is  patterned  orthogonal  to  the 
first  layer,  again  spaced  at  half  the  wavelength,  and  separated  from  the  first  layer  by 
half  the  wavelength.  The  third  layer  is  parallel  to  the  first  layer  and  placed  halfway 
between  the  rods  of  the  first  layer,  and  thus  this  pattern  can  continue  to  be  repeated 
as  desired.  New  fabrication  techniques  they  developed  while  creating  this  structure 
have  allowed  them  to  make  the  features  extremely  small.  The  molding,  filling  and 
CMP  process  allowed  them  to  fabricate  down  to  sizes  of  1.2  //in.  The  second  process 
allowed  them  to  fabricate  down  to  0.18  /un;  this  was  done  with  a  stepper  of  0.5-//m 
resolution  and  fillet  processing,  RIE,  and  CMP.  They  tested  their  crystal  with  an 
infrared  measurement  system  that  had  a  frequency  range  from  2.5-25  /un.  In  order 
to  mitigate  the  scattering  of  light  off  the  back  of  the  Si  substrate,  they  polished  it  to 
a  smoothness  of  less  than  0.3  //in.  They  extracted  the  true  attenuation  of  the  design 
by  using  an  identical  Si  substrate  as  a  base  to  reference  against  [11]. 

Other  techniques  of  three-dimensional  fabrications  and  designs  have  been  cre¬ 
ated;  John  et  al.  have  proposed  another  fabrication  technique  that  incorporates  in¬ 
terweaving  polygonal  spiral  posts  and  involves  a  custom-built  deposition  system  [23]. 
Utilizing  non-coplanar  laser  beam  techniques  on  a  photo-resist,  cross-linked  poly¬ 
mer  rods  can  be  produced  according  to  Campbell  et  al,  but  it  is  extremely  compli¬ 
cated  [23].  Another  new  approach  to  fabricating  three-dimensional  PBG  structures, 
outlined  by  Venkataraman  et  al,  uses  a  single  planar  etch  mask  with  timing  sidewall 
deep  anisotropic  RIE  etching  and  an  isotropic  etch,  followed  by  a  sidewall  polymer 
deposition  and  removal  from  the  bottom  of  spheres  so  as  to  repeat  the  process  (Fig¬ 
ures  2.9,  2.10  and  2.12)  [23].  The  key  is  the  ability  to  etch  sphere-like  voids  which 
are  sufficiently  large  so  as  to  perforate  into  the  spheres  next  to  them.  This  would  be 
analogous  to  the  two-dimensional  case  of  having  columns  of  either  air  or  dielectric 
connected  by  a  beam  in  order  to  take  advantage  of  both  the  TM  and  TE  gaps.  A 
polymer  is  deposited  to  protect  those  spheres  from  further  etching;  the  polymer  is 
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Figure  2.9:  Three-dimensional  PhC  computer  design 

of  novel  methodology  [23]. 

able  to  be  removed  from  the  bottom  of  that  layer  of  spheres  by  “smash”  etching, 
allowing  them  to  continue  etching  the  next  layer  of  spheres,  and  then  the  cycle  can 
be  repeated.  Their  calculations  lead  them  to  believe  that  they  will  create  a  complete 
three-dimensional  PBG  with  this  design  [23]. 

Miklyaev  et  al.  improves  on  the  holographic  technique  of  using  directed  laser 
beams  that  form  an  interference  pattern  in  the  photo-resist,  causing  nulls  and  peaks 
in  intensity  whose  structure  is  in  the  form  of  an  fee  lattice  [26].  The  places  within  this 
interference  structure  where  the  intensity  is  high  causes  the  photo-resist  to  effectively 
be  exposed.  This  capability  allowed  them  to  produce  lattice  features  down  to  550 
nm.  Figures  2.11  and  2.13  depict  an  SEM  image  of  a  structure  created  by  Miklyaev 
et  al.  and  the  band  gap  with  the  transmission  of  this  structure,  respectively.  They 
used  optical  spectroscopy  with  a  spot  size  of  100  //in  to  test  the  structure  [26]. 
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(1)  Apply  etch  mask  layer 


(2)  Pattern  etch  mask 


(3)  Etch  top  layer  cf  voids 


(5)  Anisotropic  smash  etch 


(7)  Passivate 


(4)  Passivate 


(6)  Etch  2nd  layer  of  voids 
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(8)  Anisotropic  smash  etch 
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Figure  2.10:  Three-dimensional  PhC  fabrication  tech¬ 
nique  of  novel  methodology.  Steps  one  through  eight  de¬ 
scribe  the  etching  process  of  this  novel  fabrication  tech¬ 
nique  [23]. 


23 


(a) 


Figure  2.11:  SEM  image  of  sample  produced  with  holo¬ 
graphic  techniques  in  photo-resist,  where  (a)  shows  the 
structure  with  a  20-/im  scale,  (b)  shows  the  structure 
with  a  5-/im  scale  and  (c)  shows  the  structure  with  a 
2-/im  scale  and  an  inserted  500-nm  scale  [26]. 
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Figure  2.12:  SEM  image  of  single-cell  and  two-layer 

PhC,  using  method  outlined  in  Figure  2.10.  (a)  Shows 
the  results  of  the  first  four  steps  in  the  novel  fabrication 
technique  of  Figure  2.10.  (b)  Shows  the  results  of  second 
four  steps  in  the  novel  fabrication  technique  of  Figure 
2.10  [23], 
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Figure  2.13:  Band  gap  and  transmission  spectrum  of 

holographic  sample  of  Figure  2.11  [26]. 
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Figure  2.14:  SEM  image  of  the  sample  fabricated  with 
GLAD.  The  left  is  a  computer  generation  of  the  individ¬ 
ual  structure  and  the  inset  is  a  magnified  section  of  the 
sample  [19]. 

A  new  technique  has  been  introduced  by  Kennedy  et  al.  that  incorporates  fabri¬ 
cating  a  three-dimensional  PhC  using  a  thin-hlm  deposition  process  that  allows  large 
PBG’s  to  be  produced  in  a  one-step  process.  The  technique,  glancing  angle  deposition 
or  GLAD,  is  outlined  in  their  paper  [19].  Figure  2.14  is  an  SEM  image  along  with 
a  computer  generated  structure  of  the  device  they  created  using  GLAD.  They  used 
a  Fourier  transform  infrared  (FT1R)  spectrometer  in  reflection  mode  to  measure  the 
properties  of  this  design.  Figure  2.15  is  a  graph  of  the  reflectance  spectrum  at  normal 
incidence,  and  the  inset  is  a  plot  of  the  reflectivity  at  shorter  wavelengths  and  was 
accomplished  using  a  Perkin-Elmer  Lambda  900  spectrophotometer  and  an  optical 
system  with  a  numerical  aperture  of  0.6  [19]. 

Auto-cloning  is  one  technique  that  has  recently  shown  some  promise  to  the 
creation  of  larger  scale  three-dimensional  PhC’s.  This  type  of  three-dimensional 
PhC  is  fabricated  by  a  combination  of  sputter-deposition  and  sputter-etching  simu- 
taneously  [17].  This  process  requires  a  two-dimensional  PhC  as  substrate  on  which 
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Figure  2.15:  Reflectance  spectrum  at  normal  incidence 
of  the  GLAD  fabricated  sample  in  Figure  2.14  [19]. 


Deposition 


Initial  shape 


(a)  (b)  (c) 

Figure  2.16:  Pictorial  of  the  auto-cloning  evolution 

showing  the  conformal  nature  created  [17]. 


to  be  sputtered.  As  the  normal  sputtering  process  is  typically  non-conformal,  with 
the  radio  frequency  (RF)  etching  process  occurring  simultaneously,  it  creates  a  more 
conformal  sputtering  process  as  can  be  seen  in  Figure  2.16.  Auto-cloning  does  not  re¬ 
quire  a  perfect  two-dimensional  structure;  imperfections  are  often  fixed  after  a  couple 
of  layers  of  sputtering  so  as  to  match  the  rest  of  the  structure’s  dimensions.  Figure 
2.17  is  an  SEM  image  of  a  two-dimensional  grating  that  has  been  auto-cloned.  After 
approximately  four  to  five  periods  of  deposition,  the  structure  is  then  resolved,  and 
that  periodicity  holds  for  the  rest  of  the  required  deposition  layers  [17]. 
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Figure  2.17:  SEM  image  of  the  auto-cloning  process 

on  a  two-dimensional  grating  [17]. 


2-4  Conclusion 

This  chapter  looked  at  the  different  types  of  PhC’s  and  the  various  fabrication 
techniques  used  to  create  two-  and  three-dimensional  PhC’s  based  on  their  frequency 
range  of  the  PBG.  In  order  to  accomplish  the  ultimate  goal  for  the  PhC  field  (three- 
dimensional  PhC’s  with  complete  PBG’s  for  frequency  range  of  the  desired  applica¬ 
tion),  a  building  block  approach  must  be  accomplished.  Based  on  the  fabrication 
techniques  accomplished  to  date,  this  thesis  will  focus  on  expanding  this  PhC  knowl¬ 
edge  base  by  examining  some  fabrication  techniques  to  create  a  two-dimensional  PhC 
with  a  PBG  in  the  visible  range.  To  aid  in  the  fabrication  process,  FDTD  software, 
used  to  model  and  simulate  various  PhC  designs,  should  be  used  to  verify  the  designs 
created  before  the  fabrication  process.  The  next  chapter  will  evaluate  two  FDTD 
softwares  for  their  use  as  sufficient  modeling  and  simulation  software  of  PhC’s. 
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III.  Two-  and  Three-Dimensional  Photonic  Crystal 

Modeling  Tools 


The  utility  to  model  and  simulate  a  PhC  design  prior  to  attempting  fabrication 
is  paramount.  The  validation  of  a  specific  PhC  design  for  a  specific  PBG  pre¬ 
vents  wasted  time  and  funds  in  fabricating  a  PhC  that  may  or  may  not  be  designed 
properly.  Also,  the  use  of  FDTD  software  for  modeling  PhC’s  allows  for  more  precise 
designs  based  on  the  application  desired,  which  can  be  simulated  with  the  software. 

Chapters  I  and  II  covered  the  history,  theory  of  PhC’s  and  some  examples  of 
current  two-  and  three-dimensional  PhC  structures  with  differing  fabrication  tech¬ 
niques,  pointing  out  the  fact  that  the  most  important  attributes  to  creating  a  PhC 
are  the  dielectric  contrast,  geometric  structure  and  the  fabrication  technique  based 
on  the  size  of  the  PhC  structure  desired.  This  chapter  evaluates  two  FDTD  softwares 
used  to  design  and  model  PhC’s.  Modeling  the  specific  PhC  design  is  a  necessary 
step  prior  to  the  fabrication  stage. 

3.1  Modeling  and  Simulation  Software  Attributes 

The  two  different  softwares  used  were  Optiwave’s  OptiFDTD  and  Photon  De¬ 
sign’s  Crystalwave.  These  commercial  softwares  solve  the  Maxwell  wave  equations 
problem  by  using  the  FDTD  method.  This  is  a  direct  numerical  solution  of  the 
time-dependent  Maxwell  curl  equations.  These  softwares  have  a  graphical  interface, 
allowing  for  ease  of  design,  simulation  and  analysis.  They  contain  many  analysis 
tools,  such  as  plane  wave  expansion  (PWE)  band  solver,  Fourier  analysis,  and  graphi¬ 
cal  representation  of  the  TE  and  TM  modes  as  they  propagate  through  the  crystal.  A 
few  basic  properties  similar  to  both  programs  will  be  demonstrated  using  Optiwave’s 
OptiFDTD  software.  Figure  3.1  is  a  graphical  user  interface  (GUI)  representation  of  a 
simple  PBG  structure  that  is  one  row  by  six  columns  of  dielectric  cylinders  (refractive 
index  3.1)  in  air.  The  radius  of  the  cylinder  is  0.2  /mi  and  the  separation  (center  to 
center)  is  1  /mi.  This  crystal  will  be  evaluated  by  an  input  wave  from  the  left  of  1.9- 
g m  wavelength  Gaussian  modulated  continuous  wave,  simulating  a  two-dimensional 
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Figure  3.1:  Single  PhC  row  of  cylindrical  columns  created  with  Opti- 
wave. 


Figure  3.2:  Depiction  of  the  Ey  wave  of  the  PhC  columns  of  Figure 
3.1. 


situation  and  analyzing  the  TE  component.  Figure  3.2  shows  the  Ey  (electric  field  in 
the  y  direction)  simulated  wave  propagation. 

Another  analysis  tool  is  the  frequency  analyzer.  Here,  the  transmittance  and 
reflectance  of  the  simulated  wave  is  normalized  to  the  input  plane  (Figure  3.3).  This 
next  example  is  how  Optiwave  simulates  a  defect  in  a  PhC  with  a  PBG.  When  prop¬ 
erly  designed,  the  PhC  is  almost  like  a  perfect  reflector.  When  the  defect  is  arranged 
in  a  line,  it  can  be  used  as  a  channel  to  guide  the  trapped  wave  (Figure  3.4).  An 
important  note  is  that  this  is  much  like  a  fiber  optic  cable  except  sharp  bends  (i.e., 
90°  or  better)  can  be  made  and  have  little  to  no  loss,  unlike  that  of  fiber  optics.  From 
Figure  3.4,  it  is  shown  that  the  defect  is  the  absence  of  the  periodic  dielectric  columns. 
The  columns  have  a  radius  of  0.3  fin i,  and  are  separated  by  1  fin i.  They  are  arranged 
in  a  two-dimensional  hexagonal  pattern.  This  is  simulated  using  a  1.9-fim  Gaussian 
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Figure  3.3:  Fourier  transform  of  the  PhC  columns  depicted  in  Figure 
3.1,  showing  the  incident  wave  (black  line)  and  the  transmitted  wave 
(red  line). 

modulated  continuous  wave  in  two  dimensions  with  respect  to  the  TE  wave.  Figure 
3.5  displays  the  normalized  frequency  analysis  of  this  design.  For  the  frequency  of 
interest,  it  has  a  wide  bandwidth  that  almost  completely  transmits  through  with  little 
loss. 

The  resulting  Poynting  vector  can  be  seen  in  Figure  3.6,  along  with  the  Ey 
field.  This  was  a  sample  of  what  OptiFDTD  can  do  in  two-dimensions.  Optiwave’s 
interface  relies  more  on  the  ability  to  construct  the  lattice,  based  on  mathematical 
formulas,  which  allows  for  an  extremely  versatile  capability  to  construct  any  lattice 
desired.  It  also  requires  a  little  more  advanced  knowledge  of  the  mathematical  for¬ 
mulas  to  implement  more  complicated  designs.  Optiwave’s  PWE  is  extremely  useful 
and  easy  to  implement.  Optiwave  also  features  the  ability  to  design  graded  indices  of 
atoms  and  allowances  for  doping  in  the  material. 

Crystalwave  has  similar  capabilities.  The  GUI  is  very  user  friendly;  however, 
this  may  present  some  limitations  in  the  complexity  of  design  allowed  to  be  created. 
Crystalwave,  when  modeled  for  a  similar  situations  as  Optiwave,  provides  similar  re- 
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Figure  3.4:  Top  view  of  dielectric  crystal  with  defect  channel  created 
in  Optiwave. 


Figure  3.5:  Frequency  component  of  the  PBG  of  the  Y-split  channel 
depicted  in  Figure  3.4. 
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Figure  3.6:  Poynting  vector  (top)  and  Ey  field  (bottom)  of  Y-split 
channel  depicted  in  Figure  3.4. 
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Figure  3.7:  (a)  Graphical  user  interface  for  the  Y-design  in  Crystal- 

wave  software,  showing  the  oxidation  layer  around  the  columns,  (b) 
Intensity  pattern  for  the  Y-design  of  (a). 

suits.  Figure  3.7  (a)  shows  the  waveguided  Y-split  using  a  PhC;  note  the  Y-split 
has  been  optimized  at  the  split  which  accounts  for  the  extra  off-sized  columns  in  the 
Y-split.  The  intensity  pattern  of  this  design  in  Crystalwave  is  modeled  in  Figure  3.7 
(b). 

Incorporated  into  Crystalwave  is  a  common  occurrence,  oxidation.  Crystal- 
wave  allows  for  a  wall  layer  to  be  modeled  on  any  atom  structure  easily;  to  do  the 
same  in  OptiFDTD  can  be  done,  but  not  as  easily.  This  wall  layer  could  be  used  to 
account  for  oxidation  that  occurs  unavoidably  in  most  situations.  Crystalwave  is  more 
confusing  when  trying  to  add  multiple  atoms  to  a  region.  Crystalwave  incorporates 
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an  interesting  feature,  allowing  the  side  cut  away  view  to  be  pictured  anywhere  in  the 
structure;  whereas  OptiFDTD  allows  a  three-dimensional  look  at  the  entire  structure 
that  can  be  rotated  about  all  three  axes.  Crystalwave  allows  not  only  the  FDTD  anal¬ 
ysis  in  two- dimensions  for  a  two-dimensional  structure,  but  also  in  three-dimensions 
for  a  two-dimensional  structure.  This  allowance  accounts  for  the  layer  thickness  of 
the  structure,  instead  of  assuming  infinite  layer  thickness.  Crystalwave  does  not  have 
a  correctly  operating  three-dimensional  band  solver.  Another  unique  feature  that 
Crystalwave  has  is  an  optimization  function.  This  feature  allows  the  computer  to 
solve  for  the  best  size  and  placement  of  designated  atoms  based  on  the  user’s  desired 
results  of  the  PhC  and  the  inputs  of  the  dependent  and  independent  variables,  which 
is  somewhat  confusing  and  requires  some  knowledge  in  C++  programming.  The  K- 
vector  of  the  band  diagram  can  be  easily  manipulated  in  Crystalwave  (i.e.,  drag  and 
drop  feature)  as  compared  to  Optiwave  (which  has  a  more  complex  method  of  doing 
the  same  thing),  allowing  for  any  wavevector  angle  to  be  evaluated  and  not  just  the 
standard  wavevectors  of  high  periodicity.  Also,  Crystalwave  allows  for  a  choice  of 
graphing  techniques  of  the  band  diagrams  where  the  y-value  can  represent  different 
characteristics  (i.e.,  frequency,  normalized  wavelength,  etc.)  Optiwave  has  a  fixed 
band  diagram  y-axis  (normalized  frequency).  Crystalwave  currently  has  no  way  of 
propagating  a  mode  in  the  y-direction.  However,  this  can  be  simulated  by  using  a 
Gaussian  mode  exciter.  To  demonstrate  this  property,  the  following  example  is  shown. 

The  best  way  to  test  this  feature  is  to  build  a  vertical  cavity  surface  emitting 
laser  (VSCEL)  in  the  modeling  software.  Figure  3.8  shows  a  VSCEL  design  that 
incorporates  a  two-dimensional  PhC  in  the  gain  medium  for  the  lateral  confinement. 
This  PhC  and  VCSEL  is  designed  for  the  wavelength  of  632.8  nm.  In  order  to  make 
this  design  easy  to  model,  the  distributed  Bragg  reflector  (DBR)  layers  are  designed 
to  be  of  optimum  refractive  index  to  facilitate  minimum  layers  to  achieve  a  0.999 
reflectance  at  the  bottom  of  the  gain  medium  and  a  0.9  reflectance  at  the  top  of  the 
gain  medium.  The  top  of  the  gain  medium  will  also  act  as  a  DBR  for  the  waveguided 
layer.  The  top  of  the  waveguided  layer  has  a  DBR  with  a  reflectance  of  0.999,  also. 
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Figure  3.8:  Side  view  of  the  VCSEL  and  the  DBR  layers  and  the 

PhC’s  used  for  horizontal  confinement. 


37 


Figure  3.9:  Top-down  view  of  the  VCSEL  design  in  the  GUI  of 

Crystalwave;  the  different  sized  columns  are  the  different  PhC’s  of  the 
waveguide  and  gain  layers  depicted  in  Figure  3.8. 


Both  two-dimensional  PhC’s  in  Figure  3.8  are  set  to  confine  the  wave  laterally 
at  a  wavelength  of  632.8  nm;  the  reason  that  the  two  PhC’s  do  not  seem  to  have  the 
same  periodicity  is  due  to  the  different  dielectrics  of  the  gain  and  waveguide  materials. 
The  waveguide  layer  has  an  index  of  2.0,  where  the  gain  medium  is  aluminum  gallium 
arsenide  (AlGaAs)  with  an  index  of  3.455.  Figure  3.9  shows  the  top-down  GUI  view 
of  the  VCSEL  cavity.  The  physical  layer  view  is  seen  in  Figure  3.10  containing  the 
corresponding  physical  (not  optical)  thickness  required  and  index  necessary  to  obtain 
the  DBR  and  VCSEL  situation  described  above.  The  waveguide  layer  has  a  defect 
channel  built  into  it,  which  can  more  readily  be  seen  in  Figure  3.9.  This  waveguide 
channel  is  a  defect  line  of  missing  columns;  the  defect  line  was  placed  so  that  the  cou¬ 
pled  lateral  mode,  due  to  the  resonance  (standing  wave)  in  the  y-direction,  is  lined 
up  with  the  coupled  mode. 

Figure  3.11  (a)  is  the  intensity  of  the  mode  as  it  propagates  down  the  waveguide 
channel.  Figure  3.11  shows  the  Ez  (b)  and  the  Ex  (c)  fields,  respectively,  propagating 
down  the  waveguide  channel.  A  reference  probe  and  a  transmitted  probe  were  placed 
at  the  opening  of  the  waveguide  from  the  cavity  and  at  the  opposite  end,  respectively. 
Based  on  these  probes,  the  following  data  was  able  to  be  assessed.  Figure  3.12  shows 
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Figure  3.10:  Physical  layers  indicating  the  thickness  in  /irn  and  re¬ 
fractive  index  of  the  VCSEL  model  in  Crystalwave. 
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Figure  3.11:  (a)  Intensity,  (b)  Ez  field  and  (c)  Ex  field  of  the  VCSEL 

simulation  in  Crystalwave  of  the  coupled  wave  propagated  down  the 
waveguide. 
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Figure  3.12:  field  of  VCSEL  simulation  in  Crystalwave  showing 

(a)  the  reference  probe  data  and  (b)  the  transmitted  probe  data. 

the  reference  and  transmitted  data  of  the  Ex  field  for  the  waveguide  channel,  and 
Figure  3.13  shows  the  net  flux  data  of  the  Ex  field  for  the  waveguide  channel.  A 
three-dimensional  image  of  the  intensity  in  the  waveguide  layer  is  shown  in  Figure 
3.14.  Figure  3.15  is  an  image  of  the  confinement  due  to  the  PhC  in  the  gain  medium 
of  this  VCSEL  design. 

In  order  to  determine  if  the  loss  in  the  waveguide  channel  is  average,  a  test 
design  was  set  up  to  model  the  waveguide  layer  as  it  was  designed  in  the  VCSEL 
simulation.  A  mode  exciter  was  placed  in  the  waveguide  layer  at  the  opening  of  the 
channel  and  pointed  down  the  channel.  Again,  probes  were  placed  at  the  beginning 
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Figure  3.13:  Ex  field  of  VCSEL  simulation  in  Crystalwave  showing 
the  net  flux  of  the  transmitted  flux  referenced  to  the  reference  flux. 


Figure  3.14:  Three-dimensional  graph  of  the  waveguide  channel  of 

the  VCSEL  simulation  in  Crystalwave. 
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Figure  3.15:  The  confinement  of  the  gain  medium  of  the  VCSEL 

simulation  in  Crystalwave. 
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Figure  3.16:  Waveguide  channel  intensity  of  the  test  waveguide  sim¬ 
ulation  in  Crystalwave. 


and  end  of  the  waveguide  channel  as  in  the  VSCEL  case.  Figure  3.16  is  the  intensity 
image  of  this  test  scenario  and  Figure  3.17  is  the  net  flux  of  the  transmitted  flux 
referenced  to  the  reference  flux. 

It  can  be  seen  that,  compared  to  the  actual  VCSEL  simulation,  the  net  flux 
is  significantly  less  than  that  of  the  test  case  where  the  mode  exciter  exists  in  the 
waveguide  region.  It  can  also  be  seen  that  the  reference  amplitude  of  the  VCSEL 
case  is  less  than  the  reference  of  the  test  case  (Figures  3.12  (a)  and  3.18).  This  could 
be  due  to  the  fact  that  the  mode  in  the  VCSEL  must  pass  through  a  0.9  DBR  be¬ 
fore  setting  up  a  resonance  in  the  waveguiding  region.  A  comparison  of  the  reference 
data  (Figures  3.12  (a)  and  3.18)  shows  that  the  starting  flux  at  the  beginning  of  the 
waveguide  channel  is  significantly  less  than  that  of  the  test  case. 

This  scenario  shows  that  using  a  Gaussian  mode  exciter  is  a  good  work-around 
in  the  Crystalwave  software  to  simulate  a  mode  in  the  y-direction.  It  also  shows  the 
possibility  of  coupling  off  of  a  VCSEL  laser  and  propagating  the  mode  down  a  hori¬ 
zontal  waveguide. 

Some  of  OptiFDTD’s  and  Crystalwave’s  analysis  results  have  already  been 
shown,  and  more  of  Crystalwave’s  results  will  be  shown  throughout  the  rest  of  the 
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Figure  3.18:  Ex  field  at  the  reference  probe  in  the  waveguide  test 

simulation  in  Crystalwave. 
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document  as  they  explain  some  of  the  modeled  situations  of  this  thesis. 

Overall,  both  softwares  are  very  useful  and  helpful  in  the  modeling  of  PhC’s. 
Optiwave’s  analysis  is  more  of  an  ideal  theory  model  and  it  was  easy  to  represent 
band  diagrams  as  previously  modeled  in  other  works.  Crystalwave’s  features  and  user 
friendly  design  are  very  nice,  and  their  results  tend  to  be  more  the  real-world  approx¬ 
imation.  Both  come  up  with  very  similar  results.  The  expectation  is  that  either  one 
of  these  software  packages  will  provide  results  showing  that  the  experiments  closely 
match  that  of  the  modeled  data.  The  goal  is  to  prove  that  these  softwares  can  produce 
results  for  the  measured  data  that  are  accurate  enough  to  validate  their  use  for  the 
design  of  future  crystal  structures. 

3.2  Conclusion 

The  tools  used  to  model  any  PhC  design  should  be  validated  against  test  results. 
Once  proper  validation  is  concluded,  the  use  of  these  tools  will  be  extremely  powerful 
for  experimentation  in  different  types  and  designs  of  PhC’s,  which  will  lead  to  better 
and  more  practical  PhC  designs  with  applications  reaching  into  more  fields.  This 
software  will  also  be  extremely  useful  in  determining  the  validity  of  a  PhC  design 
before  the  time  intensive  and  costly  fabrication  stage.  Two  different  FDTD  types  of 
software  were  evaluated  and  compared  in  this  chapter.  The  end  result  showed  that 
the  either  software  is  adequate;  however,  it  does  depend  on  the  purpose  of  the  specific 
PhC  design  being  created  and  its  intended  purpose.  This  chapter  shed  some  light 
into  this  aspect  for  these  two  softwares.  Crystalwave  will  be  used  throughout  the 
rest  of  this  thesis  as  a  continuation  of  its  evaluation  on  the  specific  designs  presented 
in  the  next  chapter.  The  next  chapter  will  demonstrate  the  design  and  fabrication 
techniques  for  the  two-dimensional  triangular  PhC  structures.  One  very  important 
aspect  to  consider  when  evaluating  an  FDTD  software  is  to  fully  understand  the 
techniques  used  for  modeling.  For  example,  Crystalwave  currently  can  do  a  three- 
dimensional  FDTD  solve;  however,  it  is  not  a  true  three-dimensional  FDTD  analysis; 
it  is  an  analysis  of  a  two-dimensional  structure  in  three-dimensions.  However,  they 
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are  working  on  the  ability  to  have  a  three-dimensional  structure  solved  using  a  three- 
dimensional  FDTD  analysis.  The  warning  is  to  be  aware  of  the  terminology  used  by 
the  software  under  question  and  fully  understand  the  techniques  employed  by  their 
analysis  techniques. 
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IV.  Experimental  Design  and  Modeling  for  various 
Fabrication  Techniques  of  Two-Dimensional  Triangular  and 
Three-Dimensional  Face-Centered-Cubic  Photonic  Crystals 

lie  creation  of  each  PhC  for  this  thesis  was  designed  using  a  gap  map  program 


_L.  for  two-dimensional  triangular  PhCs.  They  were  then  verified  by  a  band  dia¬ 
gram  program  and  modeled  (as  part  of  the  software  evaluation)  by  commercial  FDTD 
software.  The  fabrication  techniques  explored,  to  determine  viability  of  the  fabrica¬ 
tion  process  for  this  PBG  wavelength  range,  were:  focused  ion  beam,  holographic 
polymerization/lithography,  and  electron  beam  lithography.  The  criteria  for  a  PhC, 
which  was  covered  in  the  theory  as  explained  in  Chapter  I,  and  lessons  learned  from 
previous  experiments,  as  shown  in  Chapter  II,  were  used  and  considered  when  design¬ 
ing  the  two-dimensional  triangular  PhC’s  in  this  thesis.  The  process  of  designing  and 
modeling  PhC’s  based  on  the  fabrication  techniques  listed  above  is  outlined  in  this 
chapter. 

4-1  Focused  Ion  Beam  Technique  of  Two-Dimensional  Triangular  Pho¬ 
tonic  Crystals 

One  of  the  fabrication  techniques  attempted  for  the  two-dimensional  PhC’s  was 
the  focused  ion  beam  (FIB)  process.  For  this  process,  the  desire  was  to  create  a  PhC 
in  a  waveguiding  medium,  allowing  for  the  PhC  to  be  characterized  by  various  inci¬ 
dent  angles  of  a  wavelength  that  resides  in  the  PBG;  in  this  case,  the  PBG  is  designed 
for  a  helium-neon  laser  (632.8-nm)  wavelength.  Silicon  nitride  and  Amorphous  Poly 
Carbonate  (APC)  were  the  two  materials  chosen  in  which  the  PhC  would  be  etched 
using  the  FIB  process.  A  limitation  of  the  FIB  process  is  that  the  size  of  the  PhC  that 
can  be  created,  in  any  reasonable  time  frame,  is  extremely  small,  smaller  than  any 
slab  waveguided  beam.  Due  to  this  constraint,  a  mask  of  air  was  attempted,  as  well 
as  a  mask  of  gold,  in  order  to  prevent  light  from  going  around  the  PhC  and  flooding 
the  detector. 

The  gold  mask  concept  is  a  two-step  process  with  the  FIB.  Once  the  gold  strip 
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Figure  4.1:  SEM  image  of  gold  strip  with  a  10-/nn  wide  channel  cut 
by  the  FIB. 

is  deposited  by  means  of  shadow  masking  across  the  wafer  (in  a  thin  line),  the  FIB 
cuts  a  10-/U1U  channel  in  the  gold  line  as  shown  in  Figure  4.1.  In  terms  of  the  FIB, 
this  gold  line  is  large  and  it  is  extremely  time-consuming  to  cut  a  10- jum  wide  channel 
through  it.  Other  methods  to  shorten  the  time  frame  were  to  only  cut  a  short  chan¬ 
nel  into  the  gold  line.  This  would  allow  enough  room  past  the  PhC  region  that  the 
scattering  of  the  beam,  off  of  the  back  wall  of  the  gold,  could  still  be  observed  if  the 
beam  propagated  through  the  PhC.  Figure  4.2  shows  some  of  these  shorter  channels 
that  were  etched  by  the  FIB  process.  After  the  channel  is  etched  in  the  gold,  the 
wafer  would  have  silicon  nitride  or  APC  deposited  on  it.  The  FIB  is  then  used  to 
etch  the  PhC  structure  in  the  opening  of  the  channel.  This  will  enable  the  PhC  to  be 
tested  at  various  incident  angles. 
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Figure  4.2:  SEM  image  of  gold  stripe  with  a  IO-/1111  wide  channel  cut 
part  way  into  gold  by  the  FIB. 


50 


AFIT-SEM  20.0kV  12.0mm  x400  SE(U)  10/3/2005  13:1 4  '  '  '  lOOum 


Figure  4.3:  SEM  image  of  APC  with  large  rectangular  sections 

blasted  away  leaving  a  15-20-/um  wide  channel  between  them. 

The  air  mask  concept  is  a  one-step  process  for  the  FIB  technique.  After  de¬ 
positing  the  silicon  nitride  or  APC  on  the  wafer,  the  FIB  etches  away  large  pieces  of 
the  material  (silicon  nitride  or  APC)  leaving  a  15-20-/11111  channel  (Figure  4.3)  between 
the  etched  areas.  Then,  the  PhC  can  be  etched  in  the  opening  of  this  channel. 

Another  technique  initially  attempted  was  to  use  the  FIB  to  create  a  mask  to 
be  used  in  a  photolithography  process.  The  FIB  would  etch  the  PhC  structure  in  the 
resist  on  the  back  of  the  mask  so  that  the  chrome  could  then  be  etched  and  the  resist 
stripped  away.  On  the  second  mask,  the  FIB  was  used  to  etch  all  the  way  through  the 
chrome  (Figure  4.4).  This  approach  was  abandoned  due  to  time  constraints;  however, 
it  would  appear  that  diffraction  effects,  due  to  the  wavelength  of  UV  and  feature  size 
of  the  PhC,  would  prevent  this  technique  from  producing  a  PhC  of  this  size 

Before  going  forward  with  the  etching  process,  several  calculations  were  per¬ 
formed  for  each  material  to  ensure  that  light  could  be  waveguided  into  the  material 
(able  to  support  at  least  one  mode),  and  then  determine  if  a  PBG  can  be  obtained 
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Figure  4.4:  SEM  image  of  a  48-period  PhC  etched  into  the  resist  of 
a  photolithography  mask  by  the  FIB. 
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Figure  4.5:  MATHEMATICA®  graph  of  modes  supported  in  2.4  //rn 

thick  silicon  nitride. 

in  the  same  material  based  on  the  construct  of  a  two-dimensional  PhC  for  a  PBG  at 
632.8  nm.  In  order  to  determine  the  waveguiding  and  mode  propagating  capabilities, 
a  series  of  equations  in  Mathematica®  was  performed  (see  Appendix  A).  For  silicon 
nitride  that  has  an  index  of  2.0  with  a  thickness  of  2.4  /mi  deposited  on  top  of  3  //rn 
of  silicon  dioxide  (index  of  1.457)  on  a  silicon  wafer,  approximately  ten  modes  are 
supported  (Figure  4.5). 

The  prism  that  is  used  to  couple  the  beam  into  the  silicon  nitride  must  have  an 
index  greater  than  silicon  nitride  (n  =  2.0).  A  rutile  prism  has  an  index  higher  than 
that  of  silicon  nitride;  however,  it  is  a  birefringent  crystal.  This  property  of  the  rutile 
prism  makes  it  more  difficult  to  use  as  a  waveguide  due  to  the  changing  indices  of 
the  fast  and  slow  axes  based  on  the  orientation  of  the  optic  axis  in  the  prism  and  the 
incident  angle  of  the  beam  on  the  prism.  The  placement  of  the  prism,  in  regards  to 
its  distance  from  the  optics,  is  critical  to  determine  the  proper  incident  angle  into  the 
prism  necessary  to  couple  the  laser  beam  into  the  waveguiding  material;  these  angles 
are  depicted  in  Figure  4.6.  For  a  more  in-depth  look  into  the  principals  of  prism 
coupling,  see  Appendix  B.  The  angle  determined  from  the  equations  in  Appendix 
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Figure  4.6:  Illustration  for  calculating  angles  necessary  for  determin¬ 
ing  proper  prism  coupling. 

A  is  the  Om  of  Figure  4.6.  In  order  to  back  out  the  0  and  0*  of  Figure  4.6,  it  is 
necessary  to  use  the  following  equations  from  Snell’s  Law; 

0j  =  arcsin(np(sin(0m  —  45°))) 

0i  =  45°  -  0  (4.1) 

where  np  is  the  index  of  refraction  of  the  prism.  With  these  angles,  the  distances  and 
heights  required  to  obtain  the  proper  coupling  angles  can  be  calculated. 

Once  the  calculations  show  that  these  materials  can  support  a  mode  and  can 
be  prism-coupled  into,  it  is  necessary  to  determine  if  there  exists  a  band  gap  for  this 
material  in  the  wavelength  range  desired.  A  program  created  for  this  thesis  to  display 
gap  maps  for  various  materials  is  shown  in  Appendix  C.  The  gap  map  for  the  silicon 
nitride,  with  a  periodicity  of  0.35  fin i  and  varying  the  radius  of  the  triangular  columns 
from  0.1  to  0.175  /mi,  can  be  seen  in  Figure  4.7.  This  gap  map  displays  both  po¬ 
larizations’  (E  and  H)  band  gap  positions.  The  two-dimensional  triangular  structure 
with  the  periodicity,  ‘a’,  and  radius,  ‘r’,  as  stated  above,  has  no  overlaping/complete 
PBG  in  two-dimensions  as  shown  in  Figure  4.7.  For  there  to  be  an  overlapping  PBG, 
the  gap  map  would  show  the  PBG  of  each  polarization  (TE  and  TM)  overlapping.  As 
shown  in  Figure  4.7,  the  TM  polarization  (gap  between  the  red  and  blue  lines)  and 
the  TE  polarization  (gap  between  the  green  and  black  lines)  do  not  cross  over  each 
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Figure  4.7:  Two-dimensional  triangular  structured  PhC  gap  map  of 
silicon  nitride  with  a  periodicity  of  0.35  fi m  and  a  radius  range  of  0.1- 
0.175  //m  designed  with  columns  of  air  (e  =  1)  in  the  silicon  nitride 
(e  =  4)  where  the  TM  field  polarization  band  gap  is  shown  between 
the  red  and  blue  lines  and  the  TE  field  polarization  band  gap  is  shown 
between  the  green  and  black  lines. 

other.  Based  on  this  data,  the  E  field  will  be  chosen  on  which  to  base  the  design  of 
the  PhC  for  silicon  nitride.  From  the  gap  map  in  Figure  4.7,  if  u;a/27rc  =  0.5,  then 
correspondingly,  r/a  =  0.45.  By  solving  these  two  equations  for  a  A  of  632.8  nm, 
the  periodicity  (a,  center  to  center  of  the  columns)  is  316.4  nm  and  the  radius  of  the 
columns,  r,  is  142.4  nm.  The  band  diagram  program  of  Appendix  C  can  now  be  used 
to  verify  these  dimensions  for  both  the  E  and  H  fields. 

The  band  diagrams  displayed  in  Figure  4.8  show  the  E  and  H  field  polariza¬ 
tions’  band  diagrams  using  the  periodicity  of  316.4  nm  and  a  radius  of  142.4  nm 
in  silicon  nitride.  The  numbers  of  the  x  axis  in  these  band  diagrams  correspond  to 
the  momentum  vector  direction  in  the  reciprocal  lattice  of  the  crystal.  Going  from 
0  to  20  corresponds  to  going  from  K  to  Y  (in  the  reciprocal  lattice  space),  20  to  40 
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Figure  4.8:  (a)  Band  diagram  of  E  polarization  and  (b)  H  polarization 

for  silicon  nitride  with  a  =  0.3165  // in  and  r  =  0.1424  /irn,  based  on 
columns  of  air  (e  =  1)  in  the  silicon  nitride  (e  =  4).  Notice  there  is  no 
band  gap  for  this  case. 
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Figure  4.9:  Reciprocal  lattice  space  of  the  crystal;  the  left  picture 

indicates  the  wave  vector  direction  observed  and  the  right  picture  in¬ 
dicates  the  reciprocal  lattice  of  the  two-dimensional  triangular  struc¬ 
ture  [12]. 

corresponds  to  V  to  M,  and  40  to  60  corresponds  to  M  to  K.  Figure  4.9  illustrates 
the  reciprocal  lattice  space  vectors  and  image  of  the  crystal.  From  Figure  4.8,  it  is 
evident  (as  shown  by  the  gap  map)  that  there  is  a  band  gap  for  the  E  held  but  none 
for  the  H  held  based  on  these  physical  parameters  in  silicon  nitride. 

This  same  process  and  set  of  equations  hold  true  also  for  APC.  In  this  case, 
the  APC  was  l-/jm  thick  and  had  an  index  of  1.57.  By  plugging  these  numbers  into 
the  Appendix  A  program  and  using  a  prism  that  has  an  index  of  1.77  (SFL11),  the 
number  of  modes  that  can  be  supported  by  the  APC  is  two  (Figure  4.10).  For  APC, 
the  gap  map  shown  in  Figure  4.11,  depicts  E  held  polarization  band  gaps.  No  H  held 
polarization  band  gap  exists  with  the  parameters  used  for  this  gap  map.  Thus,  no 
complete  PBG  can  be  acquired  in  APC  for  these  parameters.  This  gap  map  was 
created  based  on  a  periodicity  of  a  =  1.3  //rn  and  a  radius  that  varied  from  0.5-0.64 
//rn.  To  solve  for  the  periodicity  and  radius  needed  based  on  this  gap  map  and  the 
A  desired,  choose  oja/^nc  =  0.5,  then  correspondingly  read  down  the  gap  map  to  get 
r/a  =  0.3935.  Solving  these  equations  where  A  is  632.8  nm  gives  a  periodicity  of  316.5 
nm  and  a  radius  of  124.5  nm.  The  existence  of  the  band  gap  for  this  geometry  was 
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Figure  4.10:  The  number  of  inodes  supported  in  APC  where  the  APC 

is  1-yUin  thick. 

modified  using  the  band  diagram  program  shown  in  Figure  4.12. 

For  the  photolithography  process,  the  target  material  that  the  etched  PhC 
was  to  reside  in  was  a  silicon  oxynitride  material.  This  was  deposited  by  a  plasma 
enhanced  chemical  vapor  deposition  (PECVD)  process  [2] .  The  index  of  this  material 
is  1.6  and  it  is  1.6-/im  thick.  This  index  is  close  enough  to  APC  that  the  same  PhC 
geometry  should  work  for  it  also. 

The  evaluation  of  the  Crystalwave  software  also  included  the  modeling  of  the 
band  diagrams  of  multiple  materials.  Figure  4.13  considers  the  case  of  a  material  with 
an  index  of  1.57  (like  APC  or  silicon  oxynitride).  The  Crystalwave  band  diagram  in 
Figure  4.13  (a)  can  be  compared  to  that  produced  using  the  gap  map  for  this  material 
and  is  shown  in  Figure  4.12.  Since  Crystalwave  allows  the  user  to  select  the  band 
diagram  to  be  displayed  in  units  of  wavelength  rather  than  frequency,  Figure  4.13  (b) 
shows  this  band  diagram,  also. 
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Figure  4.11:  Two-dimensional  triangular  structured  PhC  gap  map 

showing  the  E  held  of  APC  with  a  periodicity  of  1.3  //m  and  a  radius 
range  of  0.5-0.64  //nr,  designed  with  columns  of  air  (e  =  1)  in  the  APC  (e 
=  2.4649)  where  the  TM  held  polarization  band  gap  is  shown  between 
the  red  and  blue  lines  and  the  TE  held  polarization  band  gap  is  shown 
between  the  green  and  black  lines. 
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Figure  4.12:  Two-dimensional  triangular  PhC  band  diagram  of  APC 
for  a  l-/iin  thick  layer  with  a  =  0.3165  and  r  =  0.1245  based  on  columns 
of  air  (e  =  1)  in  APC  (e  =  2.4649). 
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Figure  4.13:  Band  diagram  of  a  material  with  an  index 
of  1.57  and  a  periodicity  of  316.5  nm  modeled  in  Crys- 
talwave  using  a  triangular  pattern  showing  only  a  TE 
band  gap  in  units  of  (a)  frequency  (TE  gap  of  0.487553- 
0.505439  a/A)  and  (b)  wavelength  (TE  gap  of  0.649161- 
0.626189  (/Wr1). 
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f.2  Electron  Beam  Lithography  Technique  of  Two-Dimensional  Pho¬ 
tonic  Crystals 

The  electron  beam  lithography  technique  is  being  performed  by  the  University 
of  Illinois,  Urbana-Champagne.  The  size  crystal  that  they  have  produced  before  is 
approximately  2  mm  x  2  mm,  and  the  thickest  that  they  have  etched  this  type  of 
structure  into  is  around  500  nm.  The  structure  of  the  wafer  that  they  will  etch  a  PhC 
into  is  a  GaAs  wafer  with  400-nm  growth  of  Alo.ggGao^As  and  a  300-nm  growth  of 
Alg^GaQAAs  followed  by  a  cap  of  2  nm  of  GaAs.  This  wafer  was  grown  by  University 
of  New  Mexico,  where  they  used  a  molecular  beam  epitaxy  (MBE)  process  that  is 
actually  grown  by  a  digital  superlattice  structure  of  aluminum  arsenide  (AlAs)  and 
aluminum  gallium  arsenide  (AlGaAs).  During  this  growth  process,  the  atoms  will  mix 
due  to  the  elevated  temperatures  of  over  600°  C.  The  PhC  structure  will  be  etched 
through  the  Al0  eGa0AAs  layer  into  the  Al098Ga0,02^s  layer.  Then,  the  wafer  will  be 
oxidized  so  that  the  Alg^GagwAs  layer  will  turn  to  aluminum  gallium  oxide  (Al- 
GaO),  which  will  then  be  selectively  etched  away  without  damage  to  the  PhC  in  the 
Al0  6Gao^As  layer.  The  selective  etching  of  this  new  oxidized  layer  is  accomplished 
with  a  solution  of  K0H:H20  (1  :  12  volume  ratio).  This  will  leave  a  structure  that 
has  the  Alg^Gag^As  layer  with  the  PhC  etched  in  it  suspended  in  air  (Figure  4.14). 

This  structure  will  allow  the  PhC  to  better  confine  the  modes  in  the  y-direction 
(the  direction  out  of  the  top  and  bottom  of  the  PhC)  also  providing  for  a  better  PBG, 
based  on  the  light  line  of  the  structure.  The  light  line  goes  as  c/A  and  is  basically 
the  region  that  separates  the  fast  and  slow  wave.  For  a  material  that  is  infinite  in 
the  direction  perpendicular  to  the  direction  of  propagation,  the  light  line  is  not  con¬ 
sidered;  however,  since  this  is  not  the  case,  the  light  line  needs  to  be  considered, 
as  it  will  degrade  the  capability  for  the  light  to  be  contained  in  the  material  at  the 
location  of  the  PhC.  Taking  this  into  account,  the  previously  mentioned  setup  will 
allow  the  PhC  to  be  waveguided  into  by  edge  firing  techniques.  The  technique  of 
edge  firing  will  produce  a  situation  where  the  beam  goes  around  the  material  (over 
the  top  of  the  material)  and  could  saturate  the  detector.  To  prevent  this  saturation, 
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Figure  4.14:  Illustration  of  the  final  product  of  the  AlGaAs  PhC 

design. 


the  AIq£GclqaAs  layer  will  be  etched  into  a  curved  waveguide  region  so  the  detector 
will  not  be  in  the  line  of  sight  of  the  edge  fire  beam  (Figure  4.15). 

The  PhC  structure  is  designed  for  a  PBG  at  a  wavelength  of  632.8  nm.  The 
AlGaAs  has  little  to  no  absorption  at  this  wavelength.  The  index  program  for  Al¬ 
GaAs  in  Appendix  D  provides  the  real  and  imaginary  parts  of  the  index  (imaginary 
corresponds  to  absorption).  For  Al0  eGa0AAs  at  632.8  nm  wavelength,  it  returns  only 
a  real  value  of  3.4554,  indicating  that  there  is  no  absorption  at  this  wavelength.  To 
determine  if  this  size  structure  can  support  any  waveguided  modes,  the  previously 
discussed  set  of  equations  of  Appendix  A  will  be  used  again.  Two  test  cases  must 
be  examined  to  determine  that  this  structure  is  feasible  to  support  a  waveguided 
mode.  The  first  test  case  is  where  the  top  and  bottom  cladding  of  the  waveguide 
material  (Al06Ga0AAs)  are  air  and  y4Z0.98Ga0.02^s,  respectively.  Using  this  design, 
where  Al06GaoAAs  is  300-nm  thick  with  an  index  of  3.4554,  air  has  an  index  of 
1,  and  AIo.qsGclo^As  is  400-nm  thick  with  an  index  of  3.1497,  gives  a  total  of  one 
waveguided  mode  for  this  design  (Figure  4.16  (a)).  The  second  case  is  where  the 
Al06Ga0AAs  layer  is  suspended  in  air.  Modeling  this  situation  shows  that  up  to  two 
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Figure  4.15:  Illustration  of  the  curved  waveguide  for  the  AlGaAs 

PhC  design. 

modes  can  be  supported  (Figure  4.16  (b)). 

Since  the  largest  that  the  University  of  Illinois  can  etch  is  500  nm,  the  300-nm 
design  was  chosen.  The  typical  volume  fraction  of  the  PhC  triangular  pattern  that 
the  University  of  Illinois  can  achieve  is  about  70  percent.  This  will  drive  the  values 
that  are  picked  from  the  gap  map  created  by  the  program  of  Appendix  C  for  this 
material.  The  gap  map  for  Al^^Ga^AAs  shown  in  Figure  4.17  shows  that  there  exists 
a  complete  PBG  for  this  material  and  structure  size.  In  order  to  preserve  this  70 
percent  University  of  Illinois  design  limitation,  the  values  used  are  (jja/2nc  =  0.676 
and  r/a  =  0.3501.  Solving  these  equations,  where  A  =  632.8  nm,  gives  a  periodicity 
of  427.77  nm  and  a  radius  of  149.76  nm.  Plugging  these  values  into  the  band  diagram 
equations  of  Appendix  C  shows  a  complete  band  gap  in  the  TE  polarization  and  a 
small  band  gap  for  the  TM  polarization  that  overlap  for  this  material  and  design, 
producing  a  complete  two-dimensional  PBG  as  shown  in  Figure  4.18. 

The  modeling  of  the  AIq.qGcioaAs  PhC  structure  was  accomplished  with  the 
FDTD  software,  Crystalwave,  as  part  of  the  trial  evaluation  of  this  software.  This 
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Figure  4.16:  (a)  Modes  supported  in  AloeGaoAAs  that  is  300-nrn 

thick  with  air  on  one  side  and  Al^^Gao^As  on  the  other,  (b)  modes 
supported  in  AIq, qGcloaAs  that  is  300-nrn  thick  with  air  on  both  sides. 


65 


r/a 


Figure  4.17:  Two-dimensional  triangular  PhC  gap  map  for 

AIq'qGclqaAs  with  a  =  1  //m  and  r  =  0.34  —  0.38  //in,  based  on  columns 
of  air  (e  =  1)  in  AIq.qGcloaAs  (e  =  11.94)  where  the  TM  held  polar¬ 
ization  band  gap  is  shown  between  the  red  and  blue  lines  and  the  TE 
held  polarization  band  gap  is  shown  between  the  green  and  black  lines, 
showing  a  complete  band  gap  for  this  structure. 
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Figure  4.18:  Two-dimensional  triangular  PhC  band  diagram  for 

AIqsGojqaAs  with  a  =  427.77  nm  and  r  =  149.76  nm,  based  on  columns 
of  air  (e  =  1)  in  AIq.qGcloa As  (e  =  11.94)  with  the  band  gaps  shown, 
for  (a)  E  held  polarization,  and  (b)  H  held  polarization. 
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structure  was  modeled  with  the  layers  as  described  in  the  previous  chapter.  The 
simulation  took  into  account  that  this  structure  was  not  infinite  in  the  y-direction. 
The  band  diagram,  however,  still  assumed  infinite  structure  in  the  y-direction.  The 
band  diagram  that  was  calculated  by  Crystalwave  for  a  periodicity  of  427.77  nm  is 
displayed  in  Figure  4.19  and  can  be  compared  to  that  calculated  using  the  gap  maps 
of  Appendix  C  (Figure  4.18)  for  this  material  and  structure.  Crystalwave  software 
can  also  display  this  band  diagram  as  a  function  of  wavelength  directly,  as  shown  in 
Figure  4.19  (b). 

Figure  4.20  is  the  Crystalwave  GUI  as  it  was  created  to  model  this  PhC  struc¬ 
ture.  At  the  time  of  this  modeling,  the  exact  thickness  of  the  AlogsGaomAs  layer 
had  not  been  determined;  however,  it  was  known  to  be  between  200-500  nm.  In  order 
to  “bracket”  the  modeling  possibilities,  both  of  these  dimensions  were  simulated  in 
Crystalwave.  In  the  Figures  4.21-4.24,  both  of  these  results  will  be  displayed  to  show 
the  bracketed  capabilities  of  the  PhC.  The  actual  thickness  of  the  AIqsGclqaAs  layer 
was  modeled  at  300  nm,  which  was  as  it  was  grown. 

Data  was  obtained  for  incident  waves  at  incident  angles  of  45°  and  0°,  relative 
to  the  input  plane  for  the  AIq^G^aAs  PhC  structure.  From  this  data,  it  is  concluded 
that  for  any  angle  within  the  two-dimensional  plane,  there  exists  a  PBG  at  the  632.8- 
nm  wavelength.  Figures  4.21  and  4.22  are  the  FDTD  results  of  a  45°  incidence  for 
Ex,  Ez,  and  intensity  for  200-nrn  and  500-nm  air  gaps  with  a  632.8-nm  TE  continuous 
wave  mode  exciter,  respectively,  and  Figure  4.23  is  the  FDTD  results  of  0°  incidence 
for  Ex,  Ez,  and  intensity  of  a  200-nm  air  gaps  with  a  632.8-nm  TE  continuous  wave 
mode  exciter. 

The  Fourier  transform  of  the  45°  incident  angle  case  for  the  Ex  and  Ez  fields 
of  the  reference  and  transmitted  wave  for  200-nm  and  500-nm  air  gaps  and  a  632.8-nm 
TE  continuous  wave  mode  exciter  were  also  analyzed.  Figure  4.24  shows  the  net  flux 
ratios  for  the  45°  incidence  case  with  the  200-nm  air  gap  sample  of  the  transmitted 
wave  to  the  reference  wave. 

Figures  4.21  through  4.24  also  show  that  this  design  is  a  good  PhC  for  the  TE 


Figure  4.19:  Band  diagram  of  Al06Ga0^As  modeled 

with  Crystalwave  software  for  a  periodicity  of  427.77  nm 
in  a  triangular  pattern  showing  two  TE  band  gaps  in 
units  of  (a)  frequency  (TE  gap  of  0.221144-0.322734  and 
0.636743-0.655249  a/A)  and  (b)  wavelength  (TE  gap  of 
1.934351-1.325456  and  0.671809-0.652836  (pm)-1). 
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Figure  4.20:  Graphical  user  interface  used  to  model 

the  Al0  &GaQAAs  structure  in  Crystalwave  software. 

polarization  case  in  the  visible  wavelength  range. 

4-3  Holographic  Lithography  Technique  of  Two-Dimensional  Photonic 
Crystals 

The  holographic  lithography  process  involves  a  two-stage  process  to  ensure  good 
structures.  The  first  step,  after  the  laser/optic  setup  is  complete,  is  burning  the 
first  sample  in  a  liquid  crystal /polymer  solution.  This  allows  for  a  quick  and  easy 
determination  of  the  PhC’s  periodicity.  Once  this  process  is  complete,  the  next  step 
is  to  write  this  structure  into  the  photo- resist.  The  holographic  process  is  based 
on  the  principal  of  wave  interference.  Three  laser  beams  of  specific  orientation  and 
polarization  interfere  and  create  an  interference  pattern,  which  is  a  pattern  of  high 
and  low  intensity  areas.  The  equations  for  the  E  fields  and  propagation  vectors 
were  constructed  in  Math  cad®  [20].  The  equations  and  program  are  displayed  in 
Appendix  E.  In  general,  the  three  laser  beams  are  set  at  A-specihc  angles  to  the  x-y 
plane  and  are  referenced  to  the  x  axis.  The  sample  is  in  the  x-y  plane.  The  beams 
are  tilted  off  the  perpendicular  z  axis  by  a  specified  angle,  and  each  beam  is  given 
a  specific  polarization.  One  of  the  polarizations  is  TM  polarized  (in  the  plane  of 
incidence),  while  the  other  is  TE  polarized  (perpendicular  to  the  plane  of  incidence). 
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Figure  4.21:  FDTD  results  of  a  632.8-nm  TE  mode 

exciter  at  45°  incidence  angle  with  a  200-nm  air  gap  in¬ 
cident  on  a  triangular  structure  PhC  with  a  periodicity 
of  427.77  nm  and  a  radius  of  149.76  nm  for  (a)  the  Ex 
held  and  (b)  the  Ez  held  and  (c)  the  intensity. 
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Figure  4.22:  FDTD  results  of  a  632.8-nm  TE  mode 

exciter  at  45°  incidence  angle  with  a  500-nm  air  gap  in¬ 
cident  on  a  triangular  structure  PhC  with  a  periodicity 
of  427.77  nm  and  a  radius  of  149.76  nm  for  (a)  the  Ex 
held  and  (b)  the  Ez  held  and  (c)  the  intensity. 
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Figure  4.23:  FDTD  results  of  a  632.8-nm  TE  mode 

exciter  at  0°  incidence  angle  with  a  200-nm  air  gap  inci¬ 
dent  on  a  triangular  structure  PhC  with  a  periodicity  of 
427.77  nm  and  a  radius  of  149.76  nm  for  (a)  the  Ex  field 
and  (b)  the  Ez  field  and  (c)  the  intensity. 
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Figure  4.24:  FDTD  results  of  a  632.8-um  TE  mode 

exciter  net  flux  at  a  45°  incidence  angle  with  a  200-nm 
air  gap  as  a  ratio  of  the  transmitted  flux  to  the  reference 
flux  on  a  triangular  structure  PhC  with  a  periodicity  of 
427.77  nm  and  a  radius  of  149.76  nm. 

The  third  beam  is  given  an  offset  angle  relative  to  the  TE  polarization.  Figure  4.25 
shows  the  geometry  of  this  setup.  This  pattern  is  also  dependent  on  the  wavelength 
of  the  beam  used  to  write  the  structure.  By  changing  the  angles,  </>,  9  (as  shown  in 
Figure  4.25)  and  the  polarizations  of  the  three  laser  beams,  the  proper  geometry  of 
the  interference  pattern  can  be  achieved. 

The  material  used  for  the  hnal  structure  was  SU-8  2010,  a  negative  photo¬ 
resist.  This  material  has  a  dielectric  constant  at  lOMhz  of  3.2  and  a  dielectric  loss  at 
lOMhz  of  0.007  per  the  manufacturer  of  the  photoresist.  Using  the  following  set  of 
equations,  the  index  at  the  wavelength,  632.8  nm,  can  be  determined  [13]. 

k  =  ~(L)(  2/vr) 

<f)  =  eo(///o)-fc  (4.2) 
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Figure  4.25:  Holographic  laser  beam  setup  for  two-dimensional  struc¬ 
tures  used  in  this  thesis,  with  the  propagation  vectors  of  the  three  laser 
beams  displayed  as  [20]. 
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where  L  is  the  loss  number,  k  is  rate  of  deterioration  constant,  eo  is  the  magnitude  of 
the  dielectric  constant  at  lOMhz,  f0  is  lOMhz  and  f  is  the  target  frequency. 

Plugging  in  the  numbers  to  these  equations  and  solving  for  n  by  the  equation,  n  =  i/e, 
results  in  an  index  of  1.719  at  a  wavelength  of  632.8  nm.  The  prism  coupling  tech¬ 
nique  will  be  used  to  test  this  PhC  design,  and  the  Mathematica®  equation  stated 
in  Appendix  A  will  provide  the  proper  angles  for  prism  coupling  and  the  number  of 
modes  that  can  be  propagated  in  this  material.  From  these  equations,  the  number  of 
modes  that  can  be  propagated  with  a  thickness  of  15  //in  is  over  thirty,  with  9m  =  74°. 
According  to  the  gap  map  for  SU-8  photo-resist,  there  exists  an  E  held  polarization 
gap  and  an  H  held  polarization  gap,  but  not  a  complete  two-dimensional  PBG  (Figure 
4.26).  The  E  held  polarization  band  gap  will  be  used  and  is  shown  in  Figure  4.26, 
where  the  periodicity  and  radius  are  320  nm  and  134  nm,  respectively. 

Now  that  the  periodicity  and  radius  have  been  determined,  the  laser  beam 
angles  for  the  holographic  process  must  be  determined.  From  Appendix  E,  Figure 
4.27  was  generated,  providing  a  periodicity  of  320  nm  and  an  outer  intensity  ring 
that  provides  a  radius  of  columns  to  be  134  nm  that  is  based  on  the  following  angles: 
#1,2,3  =  50°,  02  =  120°,  03  =  240°  and  a  =  10°.  The  angles  6  and  0  are  evident  as 
shown  by  Figure  4.25,  and  a  is  the  polarization  of  one  of  the  laser  beams.  For  a  more 
in-depth  understanding,  reference  Appendix  E. 

For  this  holographic  technique,  a  UV  laser  beam  of  wavelength,  364.8  nm, 
was  used.  The  liquid  crystal/polymer  solution  has  a  25  percent  solution  of  liquid 
crystal.  Thus,  for  the  liquid  crystal/polymer  solution,  the  best  separation  that  could 
be  attained  is  a  volume  fraction  of  25  percent.  This  particular  solution  has  acetone 
in  it  which  evaporates  the  liquid  crystal  once  the  intensity  of  the  interference  pat¬ 
tern  causes  the  liquid  crystal  to  separate  from  the  polymer.  This  reaction  creates 
the  columns  of  air;  however,  due  to  the  low  volume  fraction,  a  PBG  will  not  be  re¬ 
alized.  The  liquid  crystal/polymer  solution  is  simply  a  quick  and  cheap  means  to 
determine  the  periodicity  of  the  structure.  The  desired  volume  fraction  should  be 
able  to  be  accomplished  by  using  this  final  setup  with  SU-8  photo-resist.  The  liquid 
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Figure  4.26:  Two-dimensional  triangular  PhC  gap  map  and  band 

diagram  for  SU-8  negative  photo-resist,  (a)  Gap  map  with  a  periodicity 
of  1  jum  and  radius  range  of  0.3-0. 5  /un,  where  the  TE  field  polarization 
band  gap  is  shown  between  the  green  and  black  lines  and  the  TM  field 
polarization  band  gap  is  shown  between  the  red  and  blue  lines,  (b) 
Band  diagram  of  columns  of  air  (e  =  1)  in  SU-8  (e  =  2.957)  for  a 
specific  periodicity  of  320  nm  and  a  radius  of  134  nm 
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Figure  4.27:  Image  generated  in  Mathcad®  for  a  =  320  nm  and  r  = 
134  nm  in  SU-8  photo-resist,  created  from  equations  in  Appendix  E  for 
the  holographic  laser  setup  used  to  create  a  two-dimensional  triangular 
PhC  design,  the  top  and  bottom  pictures  are  different  perspectives  of 
the  same  graph. 
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Figure  4.28:  Diagram  of  the  optical  setup  used  in  the  holographic 

polymerization/lithography  process  [22], 

crystal /polymer  solution  is  25  percent  acetone,  25  percent  TL213  (liquid  crystal), 
and  50  percent  DPHDA  (penta-acrylate  polymer),  NPG  (a  cross  linking  polymer) 
and  Ciba-geigy  4265  (photo-sensitive  dye).  Furthermore,  upon  exposure  and  time, 
this  material  also  experiences  a  3  percent  shrinkage. 

The  optical  table  setup  for  the  holographic  lithography  process  is  shown  in  Fig¬ 
ure  4.28  [22].  The  ultraviolet  UV  laser  beam  is  basically  split  into  two  beams,  then 
one  of  those  beams  is  split  again;  this  process  creates  a  total  of  three  beams.  The 
UV  laser  beam  is  initially  run  through  an  adjustable  half-wave  plate  that  allows  the 
polarization  of  the  beam  to  be  adjusted  to  any  polarization.  This  will  be  necessary 
to  allow  for  all  three  beams  to  end  up  with  the  same  power  at  the  sample  surface. 
After  the  first  half-wave  plate,  the  beam  is  split  by  a  polarizing  beam  splitter.  The 
‘s’  (horizontal)  beam  is  reflected  back  to  a  mirror,  while  the  ‘p’  (vertical)  beam  is 
transmitted  through  the  splitter.  The  reflected  beam  (this  ends  up  being  beam  2 
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or  k2  from  Figure  4.25)  from  this  optic  is  then  fed  through  a  microscope  objective, 
pinhole  and  collimator,  which  provides  a  nice,  clean,  large-diameter  beam.  After  this, 
it  is  passed  through  an  adjustable  polarizer  (half-wave  plate)  to  allow  for  final  ad¬ 
justment  to  the  polarization  of  the  beam  incident  on  the  sample  surface.  The  rest 
of  the  optics  are  mirrors  used  to  achieve  the  proper  angles  as  designed  earlier.  The 
transmitted  beam  from  the  first  splitter  goes  through  an  adjustable  half-wave  plate 
again.  This  is  offset  to  a  45°  angle  to  provide  two  equal  polarizations  for  the  next 
polarization  splitter.  After  this,  it  is  also  run  through  a  microscope  objective,  pinhole 
and  collimator.  The  clean  large-diameter  beam  is  now  fed  through  the  last  polarizing 
beam  splitter.  The  transmitted  beam  (‘p’  polarized)  becomes  beam  1  or  k±  of  Figure 
4.25.  This  beam  is  fed  through  another  adjustable  half-wave  plate  to  put  fine  tuning 
on  the  polarization  as  it  is  incident  on  the  sample  surface.  Note  that  this  beam  should 
be  polarized  as  a  TE  beam  as  it  hits  the  sample  surface;  the  mirrors  used  after  the 
last  beam  splitter  rotate  the  polarization  so  that  it  is  TE  upon  the  sample  surface. 
The  reflected  beam  (‘s’  polarized)  from  the  last  beam  splitter  becomes  the  third  beam 
or  k3  of  Figure  4.25,  and  is  fed  through  an  adjustable  half-wave  plate,  allowing  the 
unique  offset  polarization  to  be  accomplished  as  it  is  incident  on  the  sample  surface. 

A  45  —  45  —  90  prism  mounted  on  an  adjustable  z-axis  plate  that  is  on  a  rota¬ 
tional  stage  was  used  to  get  the  alignment  of  the  50°  angle,  9  and  0.  The  prism  was 
adjusted  by  5°  with  the  z-axis  control,  and  that  allowed  for  proper  alignment  of  the 
beams  by  aligning  the  back  reflection  directly  with  the  incoming  beam;  this  process 
was  accurate  to  within  2°.  The  polarization  was  checked  by  placing  a  polarizer,  ad¬ 
justed  to  block  the  desired  polarization,  at  the  sample  position,  and  then  adjusting 
the  last-in-line  half-wave  plate  of  the  beam  until  the  beam  was  completely  blocked. 
Finally,  by  placing  a  power  meter  at  the  sample  location,  each  beam’s  intensity  was 
checked  and  equalized  to  each  other  by  adjusting  the  first-in-line  half-wave  plate. 
Each  beam  was  measured  to  be  3.3  mW/cm 2  0  0.1  mW/cm 2 . 

Again,  as  part  of  the  evaluation  of  the  FDTD  software,  Crystalwave  was  used 
to  model  two-dimensional  PhC’s  in  SU-8.  These  calculations  used  a  Gaussian  con- 
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Figure  4.29:  Graphical  interface  display  in  Crystal- 

wave  software  for  the  SU-8  design. 


tinuous  wave  mode  exciter  as  the  source,  and  all  the  calculations  were  accomplished 
assuming  the  boundary  conditions  to  be  infinite  in  the  y  direction  (ie.,  infinite  col¬ 
umn  lengths).  The  GUI  of  the  setup  showing  the  mode  exciter  at  a  45°  incidence  is 
shown  in  Figure  4.29. 

The  results  of  the  modeling  at  0°  and  45°  degree  incidence  were  calculated. 
Figures  4.30-4.34  are  various  FDTD  depictions  and  Fourier  analysis  graphs  for  this 
PhC  in  SU-8.  Due  to  the  boundary  conditions  of  the  simulation,  it  is  shown  in  Figure 
4.31  (a)  and  (b)  that,  for  the  Ex  case,  due  to  extra  bounces  outside  of  the  PhC,  the 
transmitted  probe  could  be  higher  than  the  reference  reading;  however,  it  is  shown 
in  Figure  4.30  (a)  that  the  wave  is  greatly  attenuated  as  it  transits  through  the  PhC. 
These  figures  indicate  that  this  is  a  good  design  for  a  PBG  at  632.8-nm  in  SU-8  photo¬ 
resist.  The  band  diagram  for  the  structure  in  the  SU-8  material  was  also  modeled 
in  the  Crystalwave  software  as  another  verification  of  the  software  being  evaluated. 
Figure  4.34  shows  the  band  diagram  of  the  SU-8  material  with  a  periodicity  of  320 
nm  and  a  radius  of  134  nm  using  a  triangular  structure.  Again,  this  figure  can  be 
compared  to  that  obtained  using  gap  maps  in  Figure  4.26  for  the  same  structure. 
This  band  diagram  is  also  shown  with  units  of  wavelength  in  Figure  4.34  (b). 
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Figure  4.30:  PhC  design  in  SU-8  with  a  periodicity 

of  320  nm  and  a  column  radius  of  134  nm  excited  by  a 
632.8-nm  Gaussian  mode  at  a  0°  incidence  showing  the 
(a)  Ex  field  and  (b)  Ez  field  and  (c)  intensity. 
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Figure  4.31:  PhC  design  in  SU-8  with  a  periodicity 

of  320  nm  and  a  column  radius  of  134  nm  excited  by  a 
632.8-nm  Gaussian  mode  at  a  0°  incidence  showing  (a) 
the  reference  Ex  amplitude  and  (b)  the  transmitted  Ex 
amplitude  and  (c)  the  reference  Ez  amplitude  and  (d)  the 
transmitted  Ez  amplitude. 
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Figure  4.32:  PhC  design  in  SU-8  with  a  periodicity 

of  320  nm  and  a  column  radius  of  134  nm  excited  by  a 
632.8-nm  Gaussian  mode  at  a  45°  incidence  showing  the 
(a)  Ex  field  and  (b)  Ez  field  and  (c)  intensity. 
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Figure  4.33:  PhC  design  in  SU-8  with  a  periodicity 

of  320  nm  and  a  column  radius  of  134  nm  excited  by  a 
632.8-nm  Gaussian  mode  at  a  45°  incidence  showing  (a) 
the  reference  Ex  amplitude  and  (b)  the  transmitted  Ex 
amplitude  and  (c)  the  reference  Ez  amplitude  and  (d)  the 
transmitted  Ez  amplitude. 
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Figure  4.34:  Band  diagram  of  SU-8  modeled  with 

Crystalwave  software  for  a  periodicity  of  320  nm  in  a 
triangular  pattern  showing  a  TE  band  gap  in  units  of  (a) 
frequency  (TE  gap  of  0.473504-0.515948  a/A)  and  (b) 
wavelength  (TE  gap  of  0.675813-0.620217  (/im)-1). 


4-4  Holographic  Lithography  Technique  of  Three-Dimensional  Face-  Centered- 
Cubic  Photonic  Crystals 

Hybrid  Technologies  of  Amherst,  New  York,  uses  a  holographic  lithography  tech¬ 
nique  to  create  three-dimensional  fee  structures  [3].  They  use  an  “umbrella”  technique 
with  their  setup,  incorporating  two  different  polarizations  of  the  beams.  Their  design 
is  adopted  from  Miklyaev  et  al.  [26]  and  Meisel  et  al.  [6].  The  basic  setup  of  the 
design  is  based  off  of  four  laser  beams.  Three  of  the  beams  are  set  up  in  much  the 
same  way  as  the  previous  holographic  technique  for  two-dimensional  PhC’s  (with  only 
two  different  polarizations).  The  fourth  beam  is  perpendicular  to  the  sample.  Figure 
4.35  shows  the  recording  geometry  used  by  Hybrid,  where  k\ ,  k2  and  k3  are  the  wave 
vectors  of  the  three  side  beams,  and  kQ  is  the  perpendicular  beam  [3]. 

These  beams  form  a  tetrahedron  with  a  base  parallel  to  the  x-y  plane.  The 
height  of  the  tetrahedron  is  governed  by  the  fourth  beam,  k$.  The  samples  they  cre¬ 
ated  are  PhC’s  designed  for  the  infrared  range.  Each  side  beam  of  the  tetrahedron 
forms  a  39°  angle  (in  the  air)  with  the  central  beam  [3].  With  regard  to  the  polar¬ 
ization,  two  different  setups  are  used.  The  first  setup  is  shown  in  Figure  4.35  (a); 
the  second  setup  is  shown  in  (b).  The  Erst  setup  illustrates  the  polarization  of  each 
beam  as  it  is  oriented  so  that  their  projections  on  the  x-y  plane  form  a  30°  angle  with 
the  x  direction.  The  second  setup  has  the  center  (perpendicular)  beam  polarized  in 
a  circular  manner,  and  the  three  side  beams  are  polarized  with  respect  to  their  plane 
of  incidence.  Figure  4.36  is  an  SEM  of  the  FCC  structures  created  by  the  “umbrella” 
technique  in  SU-8  negative  photo-resist.  The  expected  PBG  for  these  structures  is 
around  the  2.5-2.9-/mi  wavelength  range. 

4-5  Conclusion 

The  exploration  of  fabricating  a  PhC  with  a  PBG  in  the  visible  range  is  nec¬ 
essary  for  continued  advancement  the  held  of  PhC’s.  In  this  chapter,  the  design  of 
two-dimensional  triangular  PhC’s  for  fabrication  by  focused  ion  beam,  holographic 
polymerization/lithography,  and  electron  beam  lithography  techniques  were  shown. 
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Figure  4.35:  Hybrid  Technologies  laser  beam  recording  geometry  for 
three-dimensional  FCC  PhC’s.  (a)  Shows  the  setup  for  the  polarization 
of  all  the  laser  beams  projected  polarization  on  the  x-y  plane  forms  a 
30°  angle  with  the  x  axis;  (b)  shows  the  setup  for  the  polarization  of 
center  beam  to  be  circular  where  as  the  projection  of  the  side  beams  (1, 
2,  and  3)  on  the  x-y  plane  form  a  0°,  —60°  and  60°  angle  respectively, 
with  the  x  axis.  [3]. 


Figure  4.36:  SEM  images  of  FCC  structure  created  in  SU-8  by  Hy¬ 
brid  Technologies  using  holographic  process,  (a),  (b),  (c),  and  (d)  are 
different  views  and  magnifications  of  a  sample  PhC  they  fabricated  [3] . 


Depending  on  the  type  of  characterization  that  is  going  to  be  used  on  these  PhC’s, 
the  type  of  design  implemented  must  be  determined.  Therefore,  depending  on  the 
fabrication  process,  certain  preliminary  steps  must  be  accomplished  prior  to  the  cre¬ 
ation  of  the  PhC.  The  focused  ion  beam  fabrication  technique  required  the  fabrication 
of  a  mask  for  the  guided  beam,  i.e.,  the  reason  for  the  gold  strip  as  was  shown  in  this 
chapter.  The  main  point  to  grasp  from  this  chapter  is  that  to  obtain  a  PhC  with  a 
PBG  in  the  visible  wavelength  range  requires  the  size  of  the  PhC  features  to  be  on 
the  nanometer  scale.  The  possibility  for  these  designs  to  produce  a  PhC  with  a  PBG 
in  the  visible  wavelength  range  was  shown  to  be  possible  in  theory.  Next  chapter  will 
explore  the  results  of  fabricating  these  two-dimensional  triangular  structured  PhC’s 
designed  in  this  chapter  for  the  techniques  previously  listed.  Some  discussion  into  the 
setup  for  the  different  testing  techniques  will  be  discussed,  as  will  the  results  of  the 
PhC’s  that  were  able  to  be  fabricated. 
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V.  Fabrication  and  Testing  Results  and  Analysis  of  the 
Focused  Ion  Beam,  Electron  Lithography  and  Holographic 
Lithography  Fabrication  Techniques 


The  theoretical  possibility  for  PhC’s  to  be  fabricated  by  the  processes  listed  in 
Chapter  IV,  were  shown  to  be  viable.  Chapters  I  and  II  presented  the  history, 
theory  and  previous  experiments  of  PhC’s  and  were  used  as  a  basis  for  the  calculations 
of  Chapter  IV.  Chapter  III  introduced  a  necessary  design  tool,  FDTD  software,  used 
to  model  PhC  designs  to  prevent  wasted  time  and  money  during  the  fabrication  stage. 
To  actually  use  the  techniques  considered  in  this  thesis  to  fabricate  the  PhC’s  designed 
in  Chapter  IV  will  be  explored  in  this  chapter.  The  results  of  testing  these  PhC’s 
depended  on  the  ability  to  fabricate  these  structures.  The  rest  of  this  chapter  details 
the  results  of  each  fabrication  technique  outlined  in  the  previous  chapter  and  displays 
these  results. 

5.1  Focused  Ion  Beam  Fabrication  and  Testing  Results 

The  process  for  the  FIB  technique  was  outlined  in  the  previous  chapter.  The 
results  displayed  in  the  previous  chapter  are  of  the  more  successful  attempts;  however, 
the  process  involved  to  obtain  these  results  is  not  so  obvious.  The  first  gold  channel 
that  was  deposited  due  to  shadow  masking  was  extremely  wide  and  deep.  The  second 
attempt  required  making  a  much  finer  shadow  mask  and  decreasing  the  deposition 
thickness  to  lessen  the  channel  cutting  time  of  the  FIB.  This  better  sample  had 
a  gold  line  that  was  approximately  56-/iin  wide  and  1.5-pm  deep.  Figure  5.1  is  a 
magnification  of  Figure  4.1  which  shows  that  the  FIB  cut  down  into  the  substrate  on 
the  first  sample,  not  just  through  the  gold  strip.  This  unevenness  to  the  waveguide 
layer  could  cause  mode  loss  or  decoupling  properties. 

On  subsequent  trials,  more  care  was  taken  to  only  cut  down  to  the  substrate 
to  help  prevent  any  problem  this  may  cause.  Also,  from  Figure  5.1,  it  can  be  seen 
that  the  gold  appears  to  have  bubbled  away  from  the  indium  tin  oxide  (ITO)  layer. 
The  manufacturer  of  the  FIB  machine  advertises  that  there  is  no  heating  with  the 
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Figure  5.1:  SEM  image  of  gold  channel  used  in  the  FIB  process, 

showing  the  layers  penetrated  by  the  FIB.  It  is  shown  that  the  gold 
layer  actually  bubbled  and  separated  from  the  indium  tin  oxide  (ITO) 
layer  below  it. 
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FIB  process.  This  image  suggests  that  at  some  level  there  is  enough  heating  effect  to 
cause  the  problem  seen  in  Figure  5.1.  This  effect  may  be  lessened  or  negated  if  the 
ion  beam  power  is  decreased.  The  ion  beam  was  set  to  5,000  pA  for  the  cutting  of 
the  channel  and  the  channel  cutting  process  was  continually  monitored  to  determine 
when  the  gold  was  cut  through.  Figure  5.2  shows  the  initial  cut  of  the  channel;  notice 
the  irregularities  of  the  gold  surface.  The  deposition  of  the  gold  was  not  as  high  a 
quality  as  normal.  The  cutting  of  the  gold  channel  left  some  residual  gold  at  the  ends; 
this  was  cut  away  by  dropping  the  power  of  the  ion  beam  to  300  pA  and  isolating 
that  area  so  as  not  to  cut  down  into  the  substrate.  Figure  5.3  is  the  final  channel  cut 
into  the  56-yUin  wide  gold  stripe.  A  comparison  of  Figure  5.3  to  the  first  gold  channel 
cut,  Figure  4.1,  shows  the  quality  of  the  new  channel  is  obviously  better. 

The  previous  fabrications  were  designed  so  that  the  gold  would  act  as  the  mask 
for  the  waveguided  beam.  Another  attempt  was  to  replace  the  gold  stripe  with  air  for 
this  purpose.  In  other  words,  there  would  be  a  channel  of  material  (APC  or  silicon 
nitride)  between  large  expanses  of  air.  This  was  attempted  with  APC  and  is  shown 
in  Figure  5.4.  The  reason  for  the  mess  of  Figure  5.4  is  due  to  charging  in  the  FIB 
which  caused  drift  and  the  re-deposition  of  material.  The  average  drift  for  which  this 
FIB  is  rated  is  7  nm  per  minute.  The  drift  shown  in  Figure  5.4  was  much  greater 
than  that;  however,  even  7  nm  per  minute  is  enough  to  destroy  the  PhC  that  is  being 
created  due  to  its  small  size.  There  was  a  line  of  ITO  under  the  APC  area  that  was 
being  etched  by  the  FIB;  however,  the  charging  effect  due  to  the  APC  was  too  great. 
Applying  a  carbon  coating  in  order  to  control  the  charging  effect  was  not  an  option 
as  it  would  hinder  the  ability  to  prism  couple  a  beam  into  the  APC  layer.  Through 
subsequent  trials,  the  drift  was  used  to  aid  in  the  etching  process  instead  of  hindering 
it,  until  the  proper  channel  width  and  air  gap  was  achieved.  Figure  5.5  shows  this 
progression.  However,  being  able  to  manipulate  the  charging  effect  to  benefit  the 
experiment  may  be  possible  for  a  rectangular  object  of  this  large  size,  but  for  circular 
objects  on  a  small  scale  (300  nm),  it  is  a  different  matter.  For  this  reason,  it  was  not 
possible  to  fabricate  PhC  structures  of  a  visible  wavelength  range  in  this  material  by 
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Figure  5.2:  SEM  image  of  the  gold  channel  after  initial  cutting;  notice 
the  gold  residue  left  in  front  of  the  channel  that  must  be  cut  away  using 
300  pA  power. 
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Figure  5.3:  50X  microscopic  image  of  the  gold  channel  after  the 

residue  was  blasted  away. 


the  FIB. 

One  of  the  PhC  attempts  with  the  FIB  was  creating  a  metalloid  PhC  (a  PhC 
in  the  gold).  The  purpose  of  this  was  two  fold.  First,  since  creating  a  PhC  in  the 
waveguiding  material  directly  was  impossible,  this  would  demonstrate  the  ability  of 
the  FIB  to  create  a  PhC  on  that  scale.  The  second  was  to  create  a  PhC  of  the 
visible  wavelength  that,  instead  of  using  the  PBG  to  block  a  specific  wavelength 
region,  capitalizes  on  one  of  the  unique  properties  of  PhC’s  which  is  the  enhanced 
transmission  at  the  band  edges  of  the  PBG.  It  is  known  that  gold  would  block  this 
range  of  the  visible  spectrum;  however,  due  to  the  PhC,  there  can  be  transmission 
through  the  gold.  Figure  5.6  shows  an  attempt  that  was  made  in  the  first  gold  strip 
sample.  This  attempt  was  created  by  etching  the  channel  first  and  leaving  some  gold 
at  the  end,  then  etching  the  PhC  in  the  remaining  gold.  As  can  be  seen  from  Figure 
5.6,  not  only  did  the  PhC  get  etched  too  deeply,  but  the  gold  that  was  left  was  not 
a  straight  side  wall.  The  deposition  left  slanted  sides,  as  well  as  did  the  FIB  process. 
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Figure  5.4:  SEM  image  of  the  attempt  to  use  air  as  a  mask  for  the 
beam  leaving  a  10- /nn  wide  channel.  Due  to  bad  drift  and  re-deposition 
of  material,  it  can  be  seen  that  the  FIB  process  contains  too  high  of 
charging  for  this  material  and  the  effects  of  drift  cause  the  results  seen 
here,  he.,  a  non-usable  mess. 
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Figure  5.5:  SEM  image  of  APC  channel  in  air,  (a)  with  better  drift 
correction  and  (b)  with  best  drift  compensation  to  achieve  approxi¬ 
mately  15-/im  channel  of  APC  in  air.  (b)  Is  also  magnified  to  a  20-//m 
scale. 
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Figure  5.6  also  shows  a  top-down  view  of  this  occurance;  notice  the  stitching  errors 
and  the  random  mistakes  of  the  columns  created  by  the  FIB  process.  This  is  most 
likely  due  to  the  power  selected  for  the  ion  beam  during  etching. 

Figure  5.7  is  the  top-down  view  of  the  second  trial  attempt  with  the  1.5-/im 
deep  gold  line  sample.  The  stages  for  this  fabrication  started  with  the  fine  etching 
away  of  the  slanted  gold  side  left  from  the  deposition.  Next,  the  PhC  was  etched  into 
the  gold  (Figure  5.8),  which  shows  how  the  FIB  etching  process  unevenly  eats  away 
at  the  top  of  the  structure  as  the  etch  goes  deeper.  Finally,  the  etching  of  the  channel 
in  the  gold  was  accomplished  as  shown  in  Figure  5.8. 

The  first  test  sample  received  a  deposition  of  silicon  nitride,  and  the  second 
test  sample  (1.5  /irn  of  gold)  received  coatings  of  APC  as  their  waveguiding  media. 
Unfortunately,  the  second  test  sample  was  destroyed  in  the  coating  process  and  could 
not  be  put  back  into  the  FIB  for  the  final  etching  of  the  PhC  structures  in  the  APC; 
however,  it  could  still  be  evaluated  to  determine  the  waveguiding  characteristics  into 
the  channels  and  up  to  the  gold  PhC.  The  first  test  sample  was  never  taken  back 
for  its  final  etching  process;  due  to  charging  effects,  the  PhC’s  would  not  be  able 
to  be  produced.  However,  it  too  could  be  evaluated  by  waveguiding  into  the  silicon 
nitride  and  observing  the  effects  as  it  is  incident  on  the  structures.  The  process  of 
prism  coupling,  as  shown  in  Appendix  B,  portrays  the  necessity  to  have  the  prism  at 
a  higher  index  than  that  of  the  waveguide  material.  For  silicon  nitride,  this  means 
an  index  higher  than  2.0.  While  there  are  several  prisms  that  can  achieve  this,  the 
least  expensive  is  a  rutile  prism.  However,  due  to  the  birefringence  of  this  crystal, 
waveguiding  into  this  sample  was  unable  to  be  accomplished.  Waveguiding  into  the 
APC  (index  of  1.57)  is  easier,  as  a  non-birefringent  prism  could  be  used.  The  test  setup 
for  waveguiding  a  helium  neon  laser  using  a  prism  coupling  technique  is  illustrated  in 
Figure  5.9. 

The  helium  neon  laser  is  fed  through  a  pin  hole  and  a  collimator;  then,  it  is  fed 
through  an  adjustable  polarizer  that  is  set  to  only  pass  horizontal  waves  (for  this  set 
up,  creating  the  TE  polarization).  An  adjustable  aperture  is  used  to  decrease  the  size 


Figure  5.6:  SEM  image  of  metalloid  PhC  in  gold,  (a)  45°  look  angle 
at  the  metalloid  PhC  in  gold  and  (b)  top-down  view  of  metalloid  PhC 
in  gold. 
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Figure  5.7:  SEM  top-down  image  of  PhC  in  gold  where  the  gold 

strip  was  of  better  quality  and  the  etching  process  was  at  a  much  lower 
power  and  more  care  was  taken  to  try  and  not  etch  deeper  than  the 
gold  layer. 
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Figure  5.8:  SEM  image  of  a  couple  of  stages  of  the  FIB  process  for 
the  PhC  of  Figure  5.7,  where  (a)  shows  the  initial  etching  away  of  the 
gold  over  spray  and  the  etching  of  the  PhC,  and  (b)  shows  the  final 
stage  of  the  channel  being  etched  away  behind  the  PhC. 
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Figure  5.9:  Side  view  of  the  optical  setup  for  prism  coupled  waveg- 
uiding. 

of,  and  clean  up,  the  beam.  Figure  5.10  shows  a  Beamcode  image  of  the  beam  being 
used.  Next,  a  150-mm  focal  length  lens  is  used  to  focus  the  beam  to  a  point.  The 
focal  length  should  be  longer  than  the  distance  needed  to  couple  into  the  prism.  Two 
mirrors  in  a  penscope  configuration  are  used  to  raise  the  beam  to  obtain  the  proper 
angles  necessary  for  the  prism  coupling  effect.  The  top  mirror  is  on  a  micrometer 
adjustable  stage  for  fine  control  of  the  angle.  Figure  5.11  shows  the  channel  in  gold 
as  it  is  being  waveguided  to  by  a  632.8-nm  helium  neon  laser  with  an  incident  angle 
of  approximately  30°,  using  a  50X  objective  on  the  microscope.  Figure  5.11  (a)  shows 
the  channel  in  the  gold  with  the  light  on  and  (b)  shows  the  same  channel  with  the 
light  off.  Notice  the  trash  in  the  channel  from  Figure  5.11  (a);  this  should  give  a  good 
specular  reflection  if  the  beam  is  guiding  in  the  channel.  From  Figure  5.11  (b),  it 
can  be  seen  that  there  is  no  scatter  detected  from  the  channel.  This  could  be  due  to 
the  detectivity  of  the  camera  used  to  record  this  image.  A  near-held  scanning  optical 
microscopy  might  be  able  to  detect  the  presence  of  a  guided  mode  in  the  channel, 
or  light  may  not  be  guiding  down  the  channel.  Waveguiding  to  the  channel  with  a 
normal  incidence  angle  also  showed  no  scatter  down  the  channel,  indicating  no  mode 
was  successfully  propagated  down  the  channel. 

If  the  wave  is  not  propagating  down  the  channel,  it  could  be  due  to  the  coupling 
of  the  slab  waveguide  to  the  channel  waveguide.  This  situation  was  simulated  and 
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Figure  5.10:  Beamcode  image  of  the  HeNe  laser  beam  used  for  prism 
coupling. 


the  results  of  that  simulation  suggest  that  the  waveguided  mode  should  couple  into 
the  channel.  Part  of  the  evaluation  of  the  FDTD  software  included  the  modeling  of 
waveguiding.  The  modeling  results  of  the  coupling  the  slab  waveguided  beam  to  a 
channel,  in  the  gold,  is  illustrated  in  Figure  5.12.  Figure  5.12  (a)  is  the  Ex  held  (b) 
is  the  Ez  held  and  (c)  is  the  intensity,  all  at  a  0°  incidence  on  the  gold  channel.  The 
graphs  of  Figure  5.13  shows  the  amplitude  of  the  Ex  wave  at  0°  incidence  as  the  beam 
enters  the  channel,  as  shown  in  (a),  and  the  amplitude  of  the  beam  as  it  exits  the 
channel,  as  shown  in  (b).  Likewise,  Figures  5.13  (c)  and  (d)  show  the  amplitude  of 
the  Ez  wave  at  0°  incidence  as  it  enters  and  exits  the  channel,  respectively. 

Figure  5.14  (a)  is  the  Ex  held  (b)  is  the  Ez  held  and  (c)  is  the  intensity,  all  at  a 
10°  incidence  on  the  gold  channel.  The  graphs  of  Figure  5.15  shows  the  amplitude  of 
the  Ex  wave  at  10°  incidence  as  the  beam  enters  the  channel,  as  shown  in  (a),  and  the 
amplitude  of  the  beam  as  it  exits  the  channel,  as  shown  in  (b).  Likewise,  Figure  5.15 
(c)  and  (d)  show  the  amplitude  of  the  Ez  wave  at  10°  incidence  as  it  enters  and  exits 
the  channel,  respectively.  Figure  5.16  (a)  is  the  Ex  held  (b)  is  the  Ez  held  and  (c) 
is  the  intensity,  all  at  a  45°  incidence  on  the  gold  channel.  The  graphs  of  Figure  5.17 
shows  the  amplitude  of  the  Ex  wave  at  45°  incidence  as  the  beam  enters  the  channel, 
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Figure  5.11:  Gold  channel  and  APC  under  50X  magnification  with 
helium  neon  laser  beam  waveguided  to  the  channel  with  (a)  microscope 
light  on  bright  and  (b)  microscope  light  off. 
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Figure  5.12:  Crystalwave  modeling  results  of  the  cou¬ 
pling  of  the  slab  waveguide  to  the  gold  channel  where  (a) 
is  the  Ex  field  (b)  is  the  Ez  field  and  (c)  is  the  intensity 
at  a  0°  incidence  on  the  gold  channel. 


Figure  5.13:  The  amplitudes  of  (a)  the  reference  Ex 

held  (b)  the  transmitted  Ex  held  (c)  the  reference  Ez 
held  (d)  the  transmitted  Ez  held,  all  at  0°  incidence  on 
the  gold  channel. 
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Figure  5.14:  Crystalwave  modeling  results  of  the  cou¬ 
pling  of  the  slab  waveguide  to  the  gold  channel  where  (a) 
is  the  Ex  held  (b)  is  the  Ez  held  and  (c)  is  the  intensity 
at  a  10°  incidence  on  the  gold  channel. 
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Figure  5.15:  The  amplitudes  of  (a)  the  reference  Ex 

held  (b)  the  transmitted  Ex  held  (c)  the  reference  Ez 
held,  and  (d)  the  transmitted  Ez  held,  all  at  10°  incidence 
on  the  gold  channel. 
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Figure  5.16:  Crystalwave  modeling  results  of  the  cou¬ 
pling  of  the  slab  waveguide  to  the  gold  channel  where  (a) 
is  the  Ex  field  (b)  is  the  Ez  field  and  (c)  is  the  intensity, 
all  at  a  45°  incidence  on  the  gold  channel. 

as  shown  in  (a),  and  the  amplitude  of  the  beam  as  it  exits  the  channel,  as  shown 
in  (b).  Likewise,  Figure  5.17  (c)  and  (d)  show  the  amplitude  of  the  Ez  wave  at  45° 
incidence  as  it  enters  and  exits  the  channel,  respectively.  Figures  5.12  Through  5.17 
show  that  the  slab  waveguided  mode  should  clearly  couple  into  the  channel;  however, 
it  is  also  evident  that  the  magnitude  of  these  coupled  waves  is  dramatically  decreased 
and  might  possibly  require  a  better  detector  to  observe  the  scattering  effects  in  the 
gold  channel  that  was  tested.  The  modeling  depicts  the  beam  to  be  propagating  down 
the  channel  even  at  various  angles  to  the  opening.  Figure  5.18  is  a  three-dimensional 
image  of  the  intensity  of  the  beam  propagating  down  the  gold  channel.  The  source 
was  modeled  as  a  continuous  Gaussian  wave  at  632.8  nm. 
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Figure  5.17:  The  amplitudes  of  (a)  the  reference  Ex 

held  (b)  the  transmitted  Ex  held  (c)  the  reference  Ez 
held,  and  (d)  the  transmitted  Ez  held,  all  at  45°  incidence 
on  the  gold  channel. 


Figure  5.18:  A  three-dimensional  look  at  the  intensity 
of  the  wave  propagating  down  the  gold  channel. 
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Figure  5.19:  SEM  image  of  the  liquid  crystal /polymer 
sample  showing  the  periodicity  and  structure  of  the  holo¬ 
graphic  setup. 

5.2  Two-Dimensional  Holographic  Lithography  Fabrication  and  Testing 
Results 

The  holographic  technique,  as  outlined  in  the  previous  chapter,  is  a  two-step 
process  involving  a  liquid  crystal/polymer  substrate  for  initial  setup  justification,  then 
SU-8  photo-resist  for  final  structure  realization.  The  liquid  crystal /polymer  material 
is  suspended  between  two  plates  of  microscopie  slides  with  glass  cylinders  of  VS- pm. 
diameters  mixed  in  to  provide  the  proper  thickness.  The  sample  was  exposed  to  the 
UV  laser  beam  for  30  seconds.  The  top  glass  is  peeled  off  and  the  structure  left  on 
the  other  slide  is  then  imaged.  The  structure  was  coated  with  carbon  to  help  prevent 
charging;  Figure  5.19  shows  an  SEM  of  the  structure  obtained  from  the  holographic 
setup.  As  was  previously  mentioned  in  Chapter  IV  under  the  “Holographic  Lithog¬ 
raphy  Technique  of  Two-Dimensional  Photonic  Crystals”  section,  the  air  columns 
obtained  are  not  of  the  size  to  achieve  a  PhC  with  the  desired  PBG.  However,  this 
structure  is  used  to  determine  the  center-to-center  periodicity  and  geometric  structure 
of  the  holographic  setup.  Figure  5.19  shows  that  the  periodicity  is  around  320  nm  in 
a  triangular  pattern. 
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Figure  5.20:  SEM  image  of  the  SU-8  photo-resist  struc¬ 
ture  created  by  a  9.19-second  exposure  with  the  holo¬ 
graphic  setup. 


For  SU-8  photo-resist  samples,  the  required  spin-on  procedures  to  achieve  the 
15-/mi  thickness  are  listed  in  Appendix  F,  as  are  the  procedures  for  development  and 
the  exposure  time  chart.  The  actual  thickness  of  the  photo- resist  was  verified  using  an 
interferometer  technique.  The  average  thickness  was  determined  to  be  10  /mi.  This 
will  change  the  waveguided  modes  and  prism  coupling  angles  desired  to  waveguide 
into  this  material.  The  programs  in  Appendix  A  will  be  run  again  with  this  new 
thickness  to  determine  the  exact  angles  required.  The  exposure  chart  for  SU-8,  dis¬ 
played  in  Appendix  F,  is  based  on  SU-8  deposited  on  silicon;  silicon  is  too  reflective  at 
the  UV  wavelength  used  to  write  this  structure  and  destroys  the  interference  pattern; 
therefore,  a  glass  substrate  is  used.  Since  glass  is  transparent  at  these  wavelengths, 
the  exposure  times  from  the  chart  must  be  multiplied  by  a  minimum  factor  of  2  as  a 
basis,  based  on  the  advice  of  the  company  that  makes  the  SU-8.  The  exposure  times 
were  increased  from  5  seconds  to  10  seconds  in  steps  of  0.25  seconds.  It  was  found 
that  9.19  seconds  gave  the  best  structure  of  air  holes  to  material,  as  can  be  seen  in 
Figure  5.20. 

In  order  to  use  the  SU-8  as  a  waveguiding  medium  and  guide  to  the  PhC  created 
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by  the  exposure  of  the  SU-8,  two  masks  will  be  created  and  two  different  exposure 
periods  will  be  necessary.  The  first  mask  masked  off  the  entire  area  except  where  the 
three  beams  overlapped,  creating  a  circle  where  the  PhC  exists.  The  second  mask  is 
an  inverse  of  the  first  in  order  to  expose  the  rest  of  the  photo-resist  to  a  regular  beam 
(no  pattern)  so  that  it  would  not  develop  away  during  the  developing  period.  Thus, 
all  that  is  developed  away  are  the  air  columns  of  the  PhC.  The  SEM  picture  with 
an  exposure  time  of  9.19  seconds  is  shown  in  Figure  5.21  (a).  The  ruggedness  of  the 
features  in  this  image  may  be  indicative  of  an  imperfect  beam  quality.  The  periodicity 
and  triangular  structure  that  were  desired  are  evident  from  Figure  5.21  (a);  however, 
the  columns  of  air  do  not  penetrate  the  entire  thickness  of  the  photo-resist.  This 
might  be  attributed  to  the  process  used  to  scrape  the  sample  off  the  glass  substrate 
and  place  it  on  the  slide  for  the  SEM.  During  this  process,  due  to  the  small  size  of 
the  crystal  features,  the  entire  underside  of  the  crystal  could  be  caved  in  on  itself.  A 
different  exposure  time  may  be  required;  however,  Figure  5.21  (b)  shows  an  exposure 
of  9.49  seconds,  in  which  it  is  obvious  that  the  structure  is  now  being  exposed  away. 
The  thickness  of  the  sidewalls  indicate  that  a  shorter  exposure  time  is  required. 

5.3  Testing  Results  of  the  Three-Dimensional  Holographic  Lithography 
Technique 

The  three-dimensional  samples  provided  by  Hybrid  Technologies  were  charac¬ 
terized  by  performing  FTIR  measurement  set  up  for  directional  hemispheric  reflection 
measurement.  The  instrument  used  was  a  Bomen  DA3  FTIR  [5].  The  sample  was  set 
at  a  20°  angle  off  of  the  incident  beam  to  minimize  any  loss  reflected  back  down  the 
path  of  the  beam.  The  results  of  the  test  are  shown  in  Figure  5.22.  It  is  shown  that 
there  is  a  37  percent  reflection  around  3.1  pm.  This  is  not  a  band  gap;  however,  if 
it  were  a  band  gap,  there  would  be  100  percent  reflection  instead  of  only  37  percent. 
This  can  be  accounted  for  due  to  the  imperfect  fee  structure  that  Hybrid  Technologies 
created.  Since  Snell’s  law  was  not  truly  accounted  with  their  setup,  the  fee  structure 
is  slightly  ascued.  The  other  contributer  to  this  less-than-perfect  reflection  coefficient 
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Figure  5.21:  SEM  of  two-dimensional  PhC  in  SU-8 

fabricated  by  holographic  lithography  technique  with  an 
exposure  time  of  (a)  9.19  seconds  and  (b)  9.49  seconds. 
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Figure  5.22:  FTIR  directional  hemispheric  reflection 

results  of  the  three-dimensional  sample  from  Hybrid 
Technologies. 

could  be  due  to  the  number  of  periods  of  the  PhC.  After  approximately  ten  periods, 
there  is  enough  attenuation  of  the  wave  to  produce  a  PBG,  less  than  that  and  it  may 
not  have  enough  periods  to  produce  a  PBG. 

5-4  Conclusion 

The  results  of  the  FIB  fabrication  for  PhC’s  with  a  PBG  on  the  order  of  visible 
wavelength  was  shown  to  not  have  the  fidelity  necessary  to  handle  the  charging  effects 
in  any  suitable  material  for  propagating  a  visible  beam.  However,  it  was  shown  that 
the  FIB  had  the  technical  accuracy  to  create  a  PhC  on  that  scale  in  materials  in  which 
charging  was  not  a  problem.  The  two-dimensional  holographic  technique  was  shown 
to  be  able  to  produce  a  structure  of  the  correct  dimensions  necessary  to  produce  a 
PBG  in  the  visible  as  shown  by  the  modeling  in  the  previous  chapter.  However,  the 
two-dimensional  PhC  that  was  created  with  the  holographic  technique  lacked  columns 
of  air  that  extended  all  the  way  through  the  material.  This  can  be  attributed  to  the 
method  that  was  used  to  transfer  the  sample  from  the  glass  slide  to  the  post  used  in 
the  SEM.  The  three-dimensional  sample  provided  by  Hybrid  Technologies  has  some 
reflection  properties  indicating  that  it  contains  some  properties  of  a  PhC  with  a  PBG; 
however,  it  does  not  contain  a  PBG  or  there  would  be  total  reflection  for  the  PBG 
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region.  In  the  next  chapter,  the  results  of  this  thesis  will  be  summarized,  providing 
the  most  important  and  relevant  lessons  learned.  Future  work  and  recommendations 
will  also  be  discussed  to  continue  this  research  and  further  the  knowledge  in  this  field. 
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VI.  Conclusions  and  Recommendations 


This  chapter  summarizes  the  results  obtained  in  the  previous  chapters  and  out¬ 
lines  follow-on  work  in  this  field.  The  goal  was  to  explore  certain  fabrication 
techniques  for  their  viability  to  produce  two-dimensional  triangular  PhC’s  with  PBG’s 
in  the  visible  range  and  to  further  the  knowledge  base  in  this  field  regarding  PhC’s 
with  PBG’s  in  the  visible  range.  As  detailed  in  this  thesis,  the  two-dimensional  fab¬ 
rication  processes  pursued  were  FIB,  electron  lithography,  and  holographic  photo¬ 
polymerization/lithography  techniques.  The  fabrication  techniques  and  materials 
used  to  create  the  PhC  in  part  determined  the  characterization  technique  required  to 
investigate  the  PBG.  Characterization  techniques  examined  included:  the  coupling  of 
a  beam  by  means  of  a  prism  into  a  wave-guiding  medium  in  which  the  PhC  has  been 
fabricated  and  Fourier  transform  infrared  spectrometer.  The  edge  firing  technique 
was  considered  for  the  electron  beam  lithography  process.  The  FIB  and  holographic 
lithography  techniques  required  a  prism  guiding  technique  to  evaluate  the  PhC  fabri¬ 
cated,  however,  as  shown,  the  FIB  lacked  the  fidelity  to  create  PhC’s  in  the  required 
material  needed  to  accommodate  this  characterization.  The  holographic  lithography 
fabrication  process  proved  to  be  capable  to  create  a  PhC  in  the  proper  material,  how¬ 
ever,  time  constraints  prevented  further  characterization  of  this  PhC. 

The  design  of  the  PhC  was  based  on  a  program  created  to  display  gap  maps 
for  triangular  structures  and  then  the  PBG  of  the  structure,  designed  from  the  gap 
maps,  was  modeled  using  another  program  created  to  display  band  diagrams  for  trian¬ 
gular  structures,  which  was  compared  to  the  commercial  FDTD  softwares  evaluated 
in  this  thesis.  FDTD  softwares,  Crystalwave  and  OptiFDTD,  were  used  to  compare 
the  modeling  and  simulation  to  the  design.  Finally,  a  three-dimensional  PhC  was 
tested  with  the  FTIR  technique,  however,  a  complete  PBG  was  not  discovered  for 
this  sample  either.  There  was  some  merit  to  the  results  obtained. 
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6.1  Fabrication  using  the  Focused  Ion  Beam 

The  fabrication  of  the  10- pm  channel  in  the  gold  strip  was  successful;  however, 
the  waveguide  testing  conducted  showed  no  transmission  of  the  coupled  wave  into  the 
channel.  This  is  contrary  to  the  modeled  results  and  is  believed  to  be  a  detectivity 
problem  of  the  camera  being  used.  The  data  from  the  modeling  showed  a  drop  of 
about  77  percent  for  the  Ex  field  at  a  0°  incident  angle  and  a  82  percent  decrease  for 
a  10°  angle  of  incidence.  The  PhC’s  that  were  attempted  in  the  required  material 
for  prism  waveguiding  were  unable  to  be  fabricated  due  to  drift  problems  caused  by 
charging  of  the  material.  This  charging  effect  was  too  great  even  at  a  low  FIB  power 
setting  of  300  pA.  It  was  shown  that,  technically,  the  FIB  process  can  create  PhC’s 
of  the  dimensions  required  for  a  PBG  in  the  visible  range;  however,  a  material  that 
is  not  transmittable  to  the  visible  wavelength  range  was  required.  Overall,  the  FIB 
process  is  incapable  of  producing  a  PhC  of  this  size  in  a  material  in  which  a  visible 
wave  can  be  guided,  due  to  charging  effects  that  lead  to  drift  problems. 

6.2  Electron  Beam  Lithography  Fabrication 

Although,  at  this  time,  a  PhC  of  this  design  has  not  been  fabricated,  the  mod¬ 
eling  data  shows  that  the  design  will  be  PhC  with  a  37  percent  PBG.  The  modeling 
demonstrated  around  a  95  percent  attenuation  for  the  incident  angle  of  45°  for  the 
Ex  held  of  the  incident  wave.  The  slight  amount  of  wave  that  was  transmitted  could 
actually  be  an  anomaly  of  the  modeling  software.  The  defined  boundaries  allowed  for 
reflection  bounces;  this  could  account  for  the  slight  amplitude  shown  at  the  transmit¬ 
ted  boundary  of  the  PhC. 

6.3  Holographic  Two-Dimensional  Fabrication 

This  technique  is  a  viable  process  to  produce  a  PhC  with  feature  sizes  capable  to 
produce  a  PBG  in  the  visible  wavelength  range.  It  was  shown  that  a  two-dimensional 
triangular  structure  of  the  correct  periodicity  and  column  radii  was  fabricated  in  a 
negative  photo- resist  material  (SU-8).  The  columns  may  have  been  collapsed  due  to 
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the  developing  process,  where  by,  a  razor  blade  was  used  to  lift-off  the  PhC  while 
still  emerged  in  isopropyl  alcohol.  The  image  of  the  final  result  shows  the  proper 
periodicity  and  approximately  the  correct  volume  fraction  with  an  exposure  time  of 
9.19  seconds.  The  actual  spin-on  techniques  for  a  15-/im  thickness  were  shown  to 
produce  approximately  a  10- /mi  layer.  The  exposure  times  had  to  be  increased  from 
the  charts  in  Appendix  F  due  to  the  substrate  that  the  chart  was  based  on.  The 
silicon  substrate  is  what  the  photo-resist  depicts  for  its  exposure  times,  which  for 
the  UV  wavelength  used  to  create  the  PhC,  is  very  reflective;  this  would  negate  the 
interference  pattern  desired  for  the  holographic  technique.  Thus,  glass  was  used  as 
the  substrate.  Glass  will  transmit  the  UV  beam,  and  thus,  some  power  will  be  lost 
that  the  chart  assumed  would  not  be  lost.  This  accounts  for  the  increased  exposure 
times  required  to  create  this  design  on  a  glass  substrate.  The  PhC  design  which  was 
created  in  the  photo-resist  has  been  shown  by  the  modeling  software  to  produce  an 
8.5  percent  PBG  for  the  TE  wave,  with  an  82  percent  attenuation  of  the  Ex  field  at 
a  0°  incident  angle. 

6-4  Testing  of  the  Three-Dimensional  PhC 

The  testing  conducted  using  the  FTIR  shows  that  there  is  some  reflection  prop¬ 
erties  of  this  crystal  sample.  The  results  show  approximately  a  37  percent  reflection 
at  the  wavelength  of  3.1  /mi.  The  holographic  lithography  fabrication  technique  can 
produce  a  three-dimensional  PhC  as  shown,  however,  the  sample  tested  was  not  a  true 
fee  structure,  which  is  obvious  due  to  Hybrid  Technologies  fabrication  process;  thus, 
there  is  no  true  PBG.  Due  to  the  feature  size  of  this  PhC  structure,  the  existence  of  a 
PBG,  if  one  existed,  would  be  in  the  near  infrared,  which  is  in  line  with  the  reflection 
properties  shown  in  the  test  result.  While  there  is  merit  to  the  three-dimensional 
PhC  designed  by  Hybrid  Technologies,  more  work  in  this  area  is  required  to  be  able 
to  produce  a  true  three-dimensional  PBG  with  this  process.  Eventually,  the  desire 
will  be  to  use  this  process  to  create  a  three-dimensional  PhC  with  a  PBG  in  the  visible 
wavelength  range. 
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6.5  Evaluation  of  the  FDTD  Modeling  Software 

Both  softwares  evaluated  were  more  than  adequate  for  the  PhC  design  of  this 
thesis,  excluding  the  three-dimensional  PhC  design;  however,  great  care  should  be 
taken  before  purchasing  a  particular  software.  Both  gave  similar  results  for  the  designs 
attempted;  however,  the  presentation,  user-friendliness  and  analysis  features  were 
very  different.  Crystalwave  was  very  user-friendly,  presented  colorful  presentations, 
and  presented  some  unique  special  feature  applications  (he.,  thin  oxidation  walls 
around  the  PhC  columns,  randomization  feature  to  simulate  real  world  fabrication 
process,  optimizer  to  select  the  proper  size  and  placement  of  columns  based  on  the 
desired  design  goals,  etc.).  OptiFDTD  seems  to  be  more  flexible  on  the  input  of  the 
structures,  yet  not  as  user  friendly,  and  the  special  features  are  not  as  abundant  at  this 
time.  Both  companies  claim  to  be  constantly  upgrading  and  improving  their  software 
capabilities.  Support  for  both  softwares  throughout  the  trial  period  was  excellent. 

6.6  Future  Recommendations 

For  follow-on  work,  finishing  up  the  characterization  of  the  holographic  and  elec¬ 
tron  beam  techniques  should  be  accomplished.  The  holographic  lithography  process 
needs  to  be  set  up  as  a  dedicated  table  with  quality  equipment  that  is  kept  clean. 
The  laser  needs  to  have  the  best  beam  quality  possible.  Once  these  parameters  are 
obtained,  using  the  setup  outlined  in  this  thesis  along  with  the  power  settings  and  the 
9.19-second  exposure  time,  a  new  sample  should  be  fabricated.  In  order  to  prevent 
damaging  the  sample  in  developing,  instead  of  lifting-off  the  sample  to  place  it  on 
the  slide  for  the  SEM,  the  sample  should  be  left  on  the  glass  substrate  and  the  glass 
substrate  should  be  cut  down  to  the  proper  size  to  place  the  glass  with  the  PhC  on 
the  SEM  slide.  After  verification  of  the  PhC,  the  next  sample  should  be  fabricated 
using  the  same  circular  mask  as  before,  however,  this  time  after  exposing  the  PhC 
area,  the  PhC  area  should  be  masked  off  and  the  rest  of  the  photo-resist  exposed  to 
regular  UV  intensity  light.  This  will  allow  only  the  columns  of  air  to  be  developed 
away  during  the  developing  process  and  then  the  PhC  can  be  tested  by  use  of  the 
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prism  coupling  described  in  this  thesis. 

In  regards  to  the  electron  beam  lithography  process,  the  samples  of  Al0  6Ga0AAs 
that  have  been  sent  to  the  University  of  Illinois  should  have  a  channel  waveguide  of 
5-15-/im  width  etched  into  it  and  then  have  the  PhC  as  designed  in  this  thesis  etched 
into  the  waveguide.  The  sample  should  then  be  oxidized  and  selectively  etched  to 
produce  the  structure  suspended  in  air  as  described  in  this  thesis.  Edge  firing  tech¬ 
niques  should  be  used  to  determine  the  PBG  of  the  PhC  designed. 

A  software  package  like  the  ones  evaluated  in  this  thesis  should  be  purchased 
and  then  can  be  used  to  compare  the  test  results  of  the  previously  mentioned  testing 
to  validate  the  software.  Once  this  is  accomplished,  great  work  can  be  accomplished 
using  the  software  to  create  other  PhC  designs  to  guide  light,  stop  light,  block  light, 
determine  new  PhC  designs  that  utilize  new  fabrication  techniques  and  ultimately 
discover  a  new  process  for  producing  large  quantities  of  PhC’s.  Along  these  lines,  the 
software  could  be  used  to  further  the  work  started  in  designing  a  VSCEL  structure 
that  couples  longitudinally.  The  idea  is  that  the  highly  efficient  VSCEL  laser  can 
be  coupled  in  a  wave  guide  channel  with  other  VSCEL  beams  of  sightly  different 
wavelengths.  Also,  designing  PhC’s  to  have  a  broad  PBG  would  be  of  significant 
importance  to  future  work. 
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Appendix  A.  Prism  Coupling /Wave- Guiding  Calculation  Program 

This  is  the  MATHEMATICA®  program  that  calculates  the  number  of  modes  al¬ 
lowed  to  propagate  in  a  given  medium  depending  on  its  thickness  and  the  indices 
of  the  core  and  cladding  layers.  This  program  then  determines  the  roots  and  applies 
them  to  the  prism  coupling  calculations  required  and  detailed  in  appendix  B,  to  de¬ 
termine  the  incident  angle  necessary  to  accomplish  total  internal  reflections  (TIR)  in 
the  prism  allowing  wave  coupling  into  the  core  of  the  waveguide. 

A.l  Program  for  Prism  Coupling 

This  program  is  courtesy  of  Emily  Heckman.  The  inputs  are  ’nc’  the  index  of 
the  core,  ’nf’  the  index  of  one  of  the  claddings,  ’ns’  the  index  of  the  other  cladding, 
’h’  the  thickness  of  the  core.  After  the  first  part  of  the  program  is  executed,  then 
execute  the  second  part  (finding  the  root),  note:  the  value  which  is  assigned  to  look 
for  the  root  must  be  close  to  the  crossing  point  of  interest  based  on  the  graph  above 
it.  The  result  of  root  is  input  into  the  square  root  of  the  next  part.  The  equation 
following  that  result  has  the  prism  index  as  an  input  (the  denominator  of  the  Arcsin, 
and  the  result  of  the  previous  set  of  equations  is  placed  in  the  numerator.  The  answer 
that  results  is  the  6rn  of  the  prism  wave-guiding  computations. 

nc  =  1.457; 
nf  =  1.71; 
ns  =  1; 
h  =  2?r  • 

^  0.6328’ 

/3 [kappa _]  :=  y/k2nf2  —  kappa2] 
gammac[kappa_]  :=  sj /3[kappa]2  —  k2nc2\ 
gammas[kappa_]  :=  ■sj /3[kappa}2  —  k2ns2\ 
h  =  2.0; 

yl[kappaA  :=  Tan[h  kappa]] 
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v2\kmmn  1  •=  {garnmac[kappa]  +gammas  [kappa] ) 

y  yrv  Li  \Aj — J  .  /l  —  gammac[kappa]gammas[kappa]\\1 

'  ” ”  '  kappa?  ' ' 

PZof[{?/l[fcappa],  y2[kappa]},  {kappa,  0,  20},  PlotRange  —>  {  —  10,  10}] 
Run  the  program  to  get  the  plot.  Then  execute  the  rest  of  the  program. 
FindRoot[yl[kappa]  ==  y2[kappa],  {kappa,  1}] 

Execute  this  command  to  get  desired  result  of: 

{kappa  — ■>  1.43835} 

Now  plug  this  number  into  equation  below. 

N Solve yk2nf2  -  1.438352] 

Run  this  segment  to  get  the  following  answer. 

IV  Solve  [16.9179] 


Put  this  value  into  the  equation  below  along  with  the  prism  index  of  1.77  in 
this  case. 

N Solve [. ArcSin [ fyt9 ]  ^p] 

Execute  the  above  command  and  the  answer  is  the  6m  shown  below. 


NSolve[ 74.2871] 
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Appendix  B.  Prism  Coupling  Theory 


ry  his  appendix  explains  the  theory  behind  prism  coupling. 

B.  1  Theory 

Prism  coupling  is  accomplished  in  order  to  couple  a  wave  into  a  waveguiding 
medium  in  the  transverse  direction;  this  is  usually  accomplished  when  the  ability  to 
couple  into  the  waveguide  via.  edge  fire  or  butt  coupling  is  not  applicable.  In  order 
for  the  proper  coupling  to  occur,  the  phase  velocities  of  the  waves  in  the  prism  and 
the  waveguiding  medium  must  have  the  same  velocities  in  the  z-direction  [10].  This 
phase  matching  requirement  means  that: 

hm  =  kn i  sin(0m)  =  27r/A0ni  sin(6>m)  (B.l) 

The  propagation  constant  /3m  (z-component  of  the  waveguided  mode  in  the 
waveguide  medium)  must  satisfy  the  expression  f3m  >  kn\  for  waveguided  modes. 
Where  k  =  27t/A,  n\  is  the  index  of  the  cladding  that  the  coupling  is  occurring 
through  and  9m  is  the  angle  between  the  normal  to  the  waveguide  modes  propagation 
and  the  input  wave.  Figure  B.l  shows  a  depiction  of  this  situation  with  no  prism 
involved. 

Combining  the  above  two  equations;  shows  that  sin(6lrn)  >  1,  which  is  a  math¬ 
ematical  impossibility.  This  is  the  reason  a  prism  must  be  used  to  accomplish  the 
coupling  procedure.  The  prism  must  have  an  index  that  is  higher  than  that  of  the 
adjacent  layer  to  the  prism  (usually  the  cladding  just  above  the  waveguide  layer) 
and  the  waveguiding  layer  as  well.  If  a  beam  is  directed  into  the  prism  hypotonus  as 
shown  in  Figure  B.2  [10],  with  an  input  angle  such  that  the  angle  between  the  incident 
beam  on  the  inner  surface  of  the  prism  adjacent  to  the  cladding  layer  achieves  total 
internal  reflection  (TIR),  a  standing  wave  is  set  up  in  the  prism.  This  mode  is  moving 
in  the  z-direction  with  a  constant  phase  of  f3p.  Assuming  the  waveguide  region  has 
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Figure  B.l:  Diagram  of  coupling  without  a  prism. 

several  modes  that  can  exist,  then  there  are  several  /3m  (constant  phase  terms  in  the 
z-direction)  that  can  exist. 

The  key  is  to  get  one  of  these  modes  to  match  np  with  the  mode  in  the  prism 
(3p.  Both  of  these  possible  modes  (/3m,  j3p)  have  tails  that  extend  beyond  the  region  of 
the  standing  wave.  The  standing  wave  in  the  prism  has  a  tail  that  extends  into  the 
ni  layer  and  the  wave  mode  in  the  waveguide  has  a  tail  that  extends  into  the  ri\  layer 
also  (reference  Figure  B.2).  When  these  tails  overlap  in  a  constructive  manner,  an 
exchange  of  energy  takes  place  and  the  mode  is  successfully  coupled.  The  equation 
that  explains  the  matching  of  these  terms  (/3m  =  /3P)  is  shown  below. 

27mp/Aosin(0m)  =  (3m  (B.2) 
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Figure  B.2:  Diagram  of  coupling  with  a  prism. 
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Realize  of  course  that  the  equation: 


9m  >  @c  =  arcsinfYi]  /np)  (B.3) 

must  also  be  satisfied  for  TIR  to  occur,  where  9C  is  the  critical  angle  [10].  The 
interaction  between  the  prism  standing  wave  mode  and  the  waveguide  mode  only 
occurs  over  length  L  (see  Figure  B.2).  For  a  complete  interchange  of  energy  between 
the  tails  of  these  modes,  the  interaction  of  these  modes  in  the  z-direction  must  satisfy 
kL  =  7t/2,  k  is  the  coupling  coefficient.  This  coefficient  depends  on  the  indices  of  the 
prism,  cladding  and  core  and  the  prism  spacing  between  the  prism  and  the  cladding. 
Therefore  the  required  length  for  coupling  is  based  on: 

L  =  W/cos(6m)  =  tt/2  k  (B.4) 

This  can  be  solved  for  k  by  applying  a  known  L,  giving: 

k  —  n  cos(9m)/2W  (B.5) 

In  general  the  better  the  beam,  i.  e.  the  more  uniform  the  electric  field  is  across 
the  the  width  W,  resulting  in  a  better  coupled  wave.  Also  to  get  almost  100  percent 
of  the  beam  coupled,  the  trailing  edge  of  the  beam  must  exactly  intersect  the  right- 
angle  corner  of  the  prism.  If  it  is  not  perfectly  in  the  corner,  then  some  will  be  lost 
to  reflection  or  transmission. 

The  equation  that  shows  that  the  prism  must  be  of  an  higher  index  than  that 
of  the  cladding  and  the  waveguide  medium  is  shown  below: 

f3m  =  kn2  =  2n  /  \0n2  (B.6) 


It  can  be  seen  that  np  >  n2  since  the  sin(#m)  <  1. 
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Appendix  C.  Band  Diagram  and  Gap  map  Programs 


This  appendix  displays  the  band  diagram  program  in  Matlab®  created  by  Dr. 
Haus  and  the  gap  map  program  in  Matlab®  created  by  Jeff  Clark. 


C.l  Band  Diagram  Created  by  Dr.  Haus 

The  following  Matlab®  programs  are  the  E  and  H  field  band  diagrams  for 
triangular  structured  two-dimensional  PhC’s.  The  inputs  are  the  periodicity  ’d’  in 
/mi  and  the  radius  ’a’  in  //m,  along  with  the  dielectric  of  the  column  material  ’EPSA’ 
and  the  dielectric  of  the  material  ’EPSB’.  Another  parameter  that  can  be  adjusted  is 
the  number  of  plane  waves  used  to  evaluate  the  design  ’Nl’. 


Listing  C.l:  E  field  band  diagram  program. 

pendix3  /  twoDEprog2gapmapver3.  m) 


°/0  These  programs  are  for  air  columns  in  a 
7.  program  for  calculating  the  E-polarized 
triangular 

7.  lattice.  A  parameter  determining  the 

7,  waves 

clear  ;  clc ;  elf  ; 

Nl  =  9;  7.  number  of  plane  waves  (i.e.  9) 

pi  =  acos (-1.0)  ; 

7othe  lattice  constant  is  scaled  to  unity, 
column  . 


dielectric 

band  structure  of  the 
number  of  plane 


INPUT 

The  raduis  of  the 


d=.320;  7.  I  think  this  is  the  spacing  of  columns  (microns)  INPUT 
a=.134;  7»  I  think  this  is  the  radius  of  the  columns  (microns)  . 

INPUT 


7.  Volume  fraction  of  column  in  the  unit  cell, 
bbeta  =  3 . 14159* a*a/d/d*2 . 0/ sqrt (3 . 0) ; 

EPSA  =  1  ;  7.  column  DIELECTRIC  INPUT 
MatNaml  =  ’air1;  7»name  of  material 
EPSB  =2.957;  7«  material  DIELECTRIC  INPUT 
MatNam2  =  ’  SU-8’;  7»name  of  material 

7.  The  major  symmetry  points  in  the  Brillouin  zone  are  defined  (for 
7.  triangular  geometry) 


gxg  =  0.0;  gyg  =  0.0; 

gxx  =  0.0;  gyx  =  2 . 0*pi/sqrt (3 . 0) ; 

gxj  =  2.0*pi/3.0;  gyj  =  2 . 0*pi/sqrt (3 . 0) ; 

7.  The  itinerary  around  the  brillouin  zone  is  laid  out. 
gxi(l)  =  gxx; 
gyi(l)  =  gyx; 
gxi  (2 )  =  gxg; 


(ap- 
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gyi (2) 

=  gyg 

gxi (3) 

=  gXJ 

30  gyi (3) 

=  gyj 

gxi (4) 

=  gxx 

gyi (4) 

=  gyx 

35 

•/.  Construct  the  wave  vectors  used  in  the  plane  wave  calculation 
L  =  0  ; 

gxl  =  2.0*pi; 

40  gyl  =  -2 . 0*pi/sqrt (3 . 0) ; 

gx2  =  0.0; 

gy2  =  2.0*pi*2.0/sqrt(3.0) ; 

7„  generate  all  the  wavevectors  within  a  hexagon  of 

°/0  diameter  2*nl. 

45 

for  ii=-Nl : N1 ; 
for  j  =-Nl : N1 ; 

if  (  (-(Nl  +  1) <(ii  +  j ) )  &  ((ii  +  j)<(Nl  +  l) ))  , 

L  =  L  + 1 ; 

50  gx ( L )  =  gxl * ii  +  gx2*j ; 

gy (D  =  gyi*ii  +  gy2*j ; 
end  ; 
end  ; 
end  ; 
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7.  Construct  the  dielectric  matrix 

7.  the  volume  of  the  unit  cell  is  defined 

60  vcel  =  sqrt (3 . 0) /2 . 0*d*d ; 

for  i i  = 1 : L ; 

EPS (ii ,  ii)  =  bbeta/EPSA+ (1 . 0-bbeta) /EPSB ; 
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i  1  =  i  i  +  1  ; 
for  j  = i 1 : L ; 

xl  =  gx(ii)  -  gx ( j ) ; 

70  x2  =  gy(ii)  -  gy(j) ; 

x  =  sqrt (xl *xl+x2*x2) *a/d ; 

EPS (ii , j ) =2 . 0*pi * (1/EPSA -1/EPSB) *a*a/ ( vcel) *bessel j  (1  ,  . . . 
x)  /x  ; 

7.  construct  the  inverse -dielectric  matrix 

7. 

75  /  The  matrix  is  symmetric 

EPS  (  j  ,  i i  )  =  EPS  (  ii  ,  j  )  ; 

end  ; 

end 
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7o  Use  the  itinerary  constructed  above  to  move  around  the  ... 
Brillouin  zone 
n  =  L  ; 

7.  Counter  for  the  number  of  points  in  the  itinerary 
nknt  =  1  ; 
for  nk  =  1:3; 

nkl  =  nk+1 
nmax=  20; 
for  1=1 : nmax ; 

dkx  =( 1 - 1 )*( gxi ( nkl ) -gxi ( nk) )/( nmax )  +  gxi(nk); 

dky  =( 1 -1 )*( gyi (nkl ) -gyi (nk) )/( nmax)  +  gyi(nk); 
for  ii=l : n ; 

a(ii ,ii)=EPS(ii , ii ) * ( ( dkx+gx (ii))~2+( dky+gy (ii))~2)/4.0/... 
pi/pi ; 

for  j  =  ii  +  1 : n ; 

a ( ii , j ) =  EPS ( ii , j ) * ( ( dkx  +  gx ( ii ) ) * ( dkx  +  gx ( j ) ) + ( dky  +  gy . .  . 
( ii ) ) *  .  .  . 

(dky+gy (j)))/4. 0/pi/pi; 
a ( j , ii ) =a ( ii ,j) ; 

end 

end 

7. 

7.  Find  the  eignevalues 
v=eig ( a) ; 

7. 

7.  Store  the  lowest  10  eigenvalues  for  plotting 
for  li =1 : 10 ; 

rr (nknt  ,li)  =  sqrt(v(li))  ; 
kxx (nknt ) =nknt -1 ; 
end 

nknt  =  nknt  +  1  ; 

end 


end 


7.  Plot  lowest  six  eigenvalues 


f igur e  (  1 ) 
hold  on  ; 


plot  (kxx  ,  rr  ( 

,D  , 

’  r 

)  ; 

plot  (kxx  ,  rr  ( 

,2)  , 

’  b 

)  ; 

plot  (kxx  ,  rr  ( 

,3)  , 

’g 

)  ; 

plot  (kxx  ,  rr  ( 

,4)  , 

’y 

)  ; 
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plot (kxx 

rr  ( 

,5)  , 

’  r 

) 

plot (kxx 

rr  ( 

,6)  , 

’  b 

) 

130  plot (kxx 

rr  ( 

,7)  , 

’g 

) 

plot (kxx 

rr  ( 

,8)  , 

’y 

) 

plot (kxx 

rr  ( 

,9)  , 

’  r 

) 

plot (kxx 

rr  ( 

,  10) 

,  ’1 

’ 

hold  off 

135  xlabel(’r/a’,  ’font weight'  ,  ’ bold1,  ’ f ontsize ’  ,16) 

ylabel ( ’ \omegaa/2\pic 1  ,  ’ fontwe ight ’  ,  ’bold1  ,  ’fontsize 1  ,16) 

7,  title  ({’Gap  Map  for  triangular  structure  columns  of  \epsilon... 

=  ’  , int2str(EPSA)  ,  .  .  . 

7,  ’  and  material  \epsilon=  ’  ,  int2str  (EPSB)  ;  ’  E  polarization  red  .. 

and  blue,  H  polarization  green  and  black ’},’ f ontwe ight  bold ..  . 

’  ,  ’ f ontsize  ’  ,  16) 

140  title ({[’ Band  Diagram  of  E  Field  Polarization  for  triangular  ... 

structure’]  ;[’ columns  of  \ eps i lon=  ’  , num2st r ( EPSA )  ,  ’  ’,  MatNaml  .  . 

,  ’  and  material  \epsilon=’  , num2str ( EPSB )  ,  ’  ’  ,  MatNam2]}) 

Listing  C.2:  H  field  band  diagram  program. 

pendix3  /  twoDHprog2gapmapver3  .m) 

7.  program  for  calculating  the  H-polarized  band  structure  of  the  ... 
triangular 

7.  lattice.  Lines  are  inserted  to  convert  it  to  band  structure  ... 
computation 

7.  for  the  E-polarized  waves.  A  parameter  determining  the  number  ... 

of  plane 
7.  waves 

5  clear  ;  clc ;  elf  ; 

N1  =  9;  7.  number  of  plane  waves  (i.e.  9)  INPUT 

pi  =  acos (-1.0)  ; 

7«the  lattice  constant  is  scaled  to  unity.  The  raduis  of  the  ... 
column  . 

d=.320;  7.  I  think  this  is  the  spacing  of  columns  INPUT 
10  a=0.134;  %  I  think  this  is  the  radius  of  the  columns  INPUT 

7.  Volume  fraction  of  column  in  the  unit  cell, 
bbeta  =  3 . 14159* a*a/d/d*2 . 0/ sqrt (3 . 0) ; 

EPSA  =  1  ;  7.  column  DIELECTRIC  INPUT 
MatNaml  =  ’Air’;  7»name  of  material 
15  EPSB  =2.957;  '/,  material  DIELECTRIC  INPUT 
MatNam2  =  ’SU-8’;  7,name  of  material 

7.  The  major  symmetry  points  in  the  Brillouin  zone  are  defined  (for 
7.  triangular  geometry) 


20 

gxg  =  0.0;  gyg  =  0.0; 

gxx  =  0.0;  gyx  =  2 . 0*pi/sqrt (3 . 0) ; 

gxj  =  2.0*pi/3.0;  gyj  =  2 . 0*pi/sqrt (3 . 0) ; 

7.  The  itinerary  around  the  brillouin  zone  is  laid  out. 
25  gxi  ( 1 )  =  gxx  ; 


(ap- 
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gyi  (i) 

=  gyx 

gxi (2) 

=  gXg 

gyi (2) 

=  gyg 

gxi (3) 

=  gXJ 

30  gyi (3) 

=  gyj 

gxi (4) 

=  gxx 

gyi (4) 

=  gyx 

35 

•/„  Construct  the  wave  vectors  used  in  the  plane  wave  calculation 
L  =  0  ; 

gxl  =  2.0*pi; 

40  gyl  =  -2 . 0*pi/sqrt (3 . 0) ; 

gx2  =  0.0; 

gy2  =  2 . 0*pi *2 . 0/ sqrt (3 . 0) ; 

°/0  generate  all  the  wavevectors  within  a  hexagon  of 

7«  diameter  2*nl. 

45 

for  ii=-Nl : N1 ; 
for  j  =-Nl : N1 ; 

if  (  (-(Nl  +  1) <(ii  +  j ) )  &  ((ii  +  j)<(Nl  +  l) ))  , 

L  =  L  + 1 ; 

50  gx ( L )  =  gxl * ii  +  gx2*j ; 

gy  (D  =  gyi*ii  +  gy2*j ; 

end  ; 
end  ; 
end  ; 

55 

7.  Construct  the  dielectric  matrix 

7.  the  volume  of  the  unit  cell  is  defined 

60  vcel  =  sqrt (3 . 0) /2 . d0*d*d ; 

for  i i  = 1 : L ; 

EPS(ii,ii)  =  bbeta*EPSA  +  ( 1 . 0 - bbet a) * EPSB ; 

7,  H-method  insert  the  following  line  to  replace  the  ... 
preceding  line. 

7»  Construct  the  inverse  dielectric  matrix 
65  7.  EPS  (  ii  ,  ii  )  =  bbet  a/ EPSA  +  (  1 . 0  -  bbet  a) /EPSB  ; 


i  1  =  i  i  +  1  ; 

70  f  or  j  =  i  1  :  L  ; 

xl  =  gx(ii)  -  gx ( j ) ; 
x2  =  gy(ii)  -  gy ( j ) ; 
x  =  sqrt (xl*xl+x2*x2)*a/d; 

EPS ( ii  ,  j )  =  2.0*pi*(EPSA-EPSB)*a*a/(vcel)  *  bessel j  (1  ,  . . . 
x)  /x  ; 


132 


75 

80 

85 

90 

95 

100 

105 

110 

115 


"/,  H-method  insert  the  following  line  to  replace  the  .. 
preceding 

7,  line,  construct  the  inverse -dielectric  matrix 
7,  EPS(ii  ,j)=2.  dO*pi*(l/EPSA-EPSB)*a*a  /  (vcel)*besselj 

(l,x)/x; 

7.  The  matrix  is  symmetric 

EPS  (  j  ,  i i  )  =  EPS  (  ii  ,  j  )  ; 

end  ; 

end 


7,  Use  the  itinerary  constructed  above  to  move  around  the 
Brillouin  zone 


n  =  L  ; 

7.  Counter  for  the  number  of  points  in  the  itinerary 
nknt  =  1  ; 
for  nk  =  1:3; 


nkl  =  nk+1 
nmax=  20; 
for  1=1 : nmax ; 

dkx  =( 1 - 1 )*( gxi ( nkl ) -gxi ( nk) )/( nmax )  +  gxi(nk); 

dky  =( 1 -1 )*( gyi (nkl ) -gyi (nk) )/( nmax )  +  gyi(nk); 

for  ii=l : n ; 


a ( ii  ,  ii )  =  ( ( dkx  +  gx ( ii ) ) ~2  +  ( dky  +  gy ( i i ) ) ~ 2 ) /4 . 0/ pi /pi ; 
b(ii,ii)  =  EPS(ii,ii); 
for  j=ii+l:n; 

a  (  ii  ,  j  )  =  0.0; 

a( j , ii)  =  0.0; 
b ( ii , j )  =  EPS  (  ii  ,  j  )  ; 
b  ( j  ,  ii )  =  EPS ( ii  ,  j  )  ; 


end 


end 

7,  H-method  insert  the  following  lines  and  replace  the  ... 

preceding  10 
7,  lines 

7,  for  i  i  =  1  :  n  ; 

7.  a(ii  ,ii)=EPS(ii  ,  ii  )  *  (  (  dkx  +  gx  (ii))~2+(  dky  +  gy  (ii))~2)/4.0/. 

pi/pi ; 

7,  for  j  =  i  i  + 1  :  n  ; 

7.  a  (  ii  ,  j  )  =  EPS  (  ii  ,  j  )*((  dkx  +  gx  (  ii  ))*(  dkx  +  gx  (j  ))  +  (  dky  +  gy  (  . 

ii ) ) *  .  .  . 

7.  (dky  +  gy  (j  ))) /4 . 0/ pi /pi  ; 

7.  a  (  j  ,  ii  )  =a  (  ii  ,  j  )  ; 

7,  end 

7,  end 


7.  Find  the  eignevalues 
v  =  eig (a , b)  ; 

7.  Find  the  eignevalues  for  the  H-method  replace  the 
preceding  line 
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7.  v=eig  (a)  ; 

7,  Store  the  lowest  10  eigenvalues  for  plotting 
for  li=l : 10 ; 

rr (nknt ,li)=sqrt(v(li)) ; 
kxx (nknt ) =nknt -1 ; 
end 

nknt  =  nknt  +  1  ; 

end 

end 


7,  Plot  lowest  six  eigenvalues 


f igur e  (  1 ) 
hold  on  ; 
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xlabel(’r/a’  ,  'font weight  '  ,  'bold'  ,  ’ f ontsize  '  ,16) 
ylabel('\omegaa/2\pic’  ,  'f  ontwe ight ’  ,  'bold'  ,  'fontsize '  ,16) 

7.  title  ({’Gap  Map  for  triangular  structure  ';'  columns  of  \epsilon... 

=  '  , int2str(EPSA)  ,  .  .  . 

I  '  and  material  \epsilon= ' , int2str (EPSB) ; ’ E  polarization  red  ... 

and  blue,  H  polarization  green  and  black f ontwe ight  bold ..  . 

’  ,  ’ fontsize  ’  ,  16) 

title ({[’ Band  Diagram  of  H  Field  polarization  for  triangular  ... 

structure']  ;[’ columns  of  \ eps i lon=  ’  , num2st r ( EPSA )  ,  ’  ',  MatNaml  .  . . 

,  ’  and  material  \epsilon=’  , num2str ( EPSB )  ,  ’  ’  ,  MatNam2]}) 


C.2  Gap  Map  Program  Created  by  Jeff  Clark 

The  following  program  was  designed  to  allow  gap  maps  of  two-dimensional  tri¬ 
angular  structured  PhC’s  to  be  created  based  off  of  the  previous  two  programs  by  Dr. 
Haus  (being  used  as  subroutines).  Input  parameters  are  ’Nl’  number  of  plane  waves 
used  to  evaluate  the  desired  structures,  ’EPSA’  the  dielectric  of  columns,  ’EPSB’  the 
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dielectric  of  the  material,  ’NDC’  number  of  dispersion  curves  calculated,  ’Prec’  the 
precision  of  the  PBG  desired,  ’a  range’  the  radius  range  for  the  calculations  of  the 
gap  map,  ’a  range  increments’  the  value  that  the  ’a  range’  is  incremented  from  one 
end  of  the  range  to  the  other. 

Listing  C.3:  E  field  band  diagram  program.  ( app endix3 /main,. outinev er 5 . m) 

“/.main  routine 
pack 

clear ; clc ; elf  ; 


”/o  Input  s 

d  =  l ;  '/.lattice  spacing 
0/Oa=.48;  '/.radius 

Nl=9;  l  number  of  plane  waves 
EPSA  =  1;  ’/.epselon  of  columns 
MatNaml  =  ’Air’;  '/.name  of  material 
EPSB  =  11.94;  “/.  epselon  of  material 

MatNam2  =  ’AlGaAs’;  '/.name  of  material 
NDC  =  10;  '/.number  of  dispersion  curves 
Prec  =  2;  %  precision  of  the  PBG 

a.range  =  [.34  .38]  ;  %  range  of  radius  of  columns 

a_range_increments  =  .01;  "/.  increment  range  radius  of  columns 

7,  end  input  s 


°/.  determine  radius  to  send  in  subroutine 
count  =  0 ; 

for  t  =  a_r ange ( 1 ): a_range_ increment s : a.range (2) 
a  =  t  ; 

count  =  count  +1; 


“/call  sub  routine  for  E  and  H  field 

rr  =  Ef  ie  lddi  sper  s  ion  (N1  ,  d  ,  a  ,  EPSA  ,  EPSB  ,  NDC  )  ; '/.  this  is  actually... 
the  H  field  subroutine 

ss  =  Hf  ie  lddi  sper  s  ion  (N1  ,  d  ,  a  ,  EPSA  ,  EPSB  ,  NDC  );  °/n  this  is  ... 
actually  the  E  field  subroutine 


“/„  precision  of  band 
factor  =  10~Prec; 

gap 

gg  =  rr.*f actor;  % 

this  is  s 

curve 

gg  =  round (gg);  % 

rr  =  gg . /factor ; 

numbers  t 

factor  =  10'Prec; 
gg2  =  ss . *  f  actor ;  ) 

{  this  is 

curve 

gg2  =  round (gg2) ; 
ss  =  gg2. /factor; 

•/.  numbers 

o  we  can  round  the  dispersion 
o  the  percsion  requested 

so  we  can  round  the  dispersion 
to  the  percsion  requested 


135 


45 

50 

55 

60 

65 

70 

75 

80 


%  calculations  for  PBG 
for  z  =  1  :  NDC  - 1 

bandH  =  max(rr(:  ,  z))  ;  °/0  finds  the  max  number  along  this  . 

dispersion  curve 

bandL  =  min  ( rr  (  :  ,  z  + 1 )  )  ;  °/0  finds  the  min  number  along  the  . 
disperision  curve 

if  bandL-bandH  >  0  %  this  determins  if  there  is  a  gap  or 

not 

PBGLow ( count , z ) =  bandH;  7«  now  we  store  the  results 
PBGHigh ( count , z) = bandL ; 

else 

PBGLow(count ,z)=0; 

PBGHigh ( count ,z)=0; 


end 

z  =  z  + 1  ; 

end 

1  now  for  the  H  field 
for  z  =  1  :  NDC  - 1 

HbandH  =  max(ss(:  ,  z))  ;  °i  finds  the  max  number  along  this 
dispersion  curve 

HbandL  =  min  ( ss  (  :  ,  z  + 1 )  )  ;  °i  finds  the  min  number  along  the 
disperision  curve 

if  HbandL  -  HbandH  >  0  °l  this  determins  if  there  is  a  gap 
or  not 

HPBGLow  ( count  ,  z )  =  HbandH;  ”/0  now  we  store  the  results 
HPBGHigh ( count , z) = HbandL ; 

else 

HPBGLow  ( count  ,  z )  =  0;  °/0  now  we  store  the  results 

HPBGHigh ( count ,z)=0; 

end 

z  =  z  + 1  ; 

end 

x(count)  =  (  r  ound  (  (  a/d)  .  *  f  act  or  )  )  .  /  f  act  or  ;  °i  this  records  the 
r/a  factor  for  the  x  axis 

end 


•/.  this  gets  rid  of  the  zeros  and  records  only  integer  numbers  for 
the  most 
7.  part 
u  =  1  ; 
ul  =  1  ; 

for  c  =  1 : count 
u  =  1  ; 
ul  =  1  ; 

for  v  =  1  :  s ize  ( PBGLow  ;  ) 
if  PBGLow(c,v)  >  0 

y(c,u)  =  PBGLow(c,v); 
u  =  u  +  1  ; 
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end 

if  PBGHigh(c,v)  >  0 

y2(c,ul)  =  PBGHigh ( c , v) ; 
ul  =  ul  +  1 ; 

end 

end 

end 

'/.This  is  the  same  as  above  but  for  the  H  field 
u  =  1  ; 
ul  =  1  ; 

for  c  =  1 : count 
u  =  1  ; 
ul  =  1  ; 

for  v  =  1  :  s ize  ( HPBGLow  ’  ) 
if  HPBGLow(c,v)  >  0 

Hy(c,u)  =  HPBGLow (c , v) ; 
u  =  u  +  1  ; 

end 

if  HPBGHigh (c , v)  >  0 

Hy2 ( c , ul )  =  HPBGHigh (c ,v) ; 
ul  =  ul  +  1  ; 

end 

end 

end 

7o  determine  what  to  plot 

[i,j]  =  f ind ( PBGLow )  ;  %  finds  the  indicies  of  the  nonzero  integers 

j(max(size(j))+l)=  0;  %  add  a  0  to  end  of  string  to  use  as  stop  .. 

bit  in  next  code 

i  (max  ( s  ize  ( i )  )  + 1 )  =  0  ;  °/0  add  a  0  to  end  of  string  to  use  as  stop  bit 
in  next  code 
counter  =  1 ; 

7o  now  same  as  above  for  the  H  field 

[Hi,Hj]  =  f ind ( HPBGLow ) ;  7  finds  the  indicies  of  the  nonzero  ... 
integers 

Hj  (max  ( size  (Hj  ))  +1)  =  0;  7  add  a  0  to  end  of  string  to  use  as  stop 

bit  in  next  code 

Hi(max(size(Hi))+l)=  0 ; 7  add  a  0  to  end  of  string  to  use  as  stop  . 

bit  in  next  code 
Hcounter  =  1 ; 


7  plot 

7plots  the  dispersion  curves  on  the  gap  map 
f igure (2) 

while  counter  <  max ( s ize ( j ) ) - 1  7  make  sure  we  don’t  violate  the 

max  length  of  the  matrix 
newcounter  =  1 ; 
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while  isequal (j ( counter ), j ( counter +1) )  %  finds  all  column  ... 
numbers  so  they  can  be  ploted  with  corresponding  x  values 
if  counter  <  max ( s ize ( j ) ) - 1 

inew ( newcount er )  =  i(counter); 
j new ( newcount er )  =  j (counter); 
counter  =  counter  +1; 
newcounter  =  newcounter +1 ; 

end 

end 

inew ( newcount er )  =  i(counter); 
j new ( newcount er )  =  j (counter); 

plot(x(inew(l:end))  ,PBGLow(inew(l:end)  , j  new ( 1 : end ) )  ,  ’  r  ’  ) 
hold  on 

plot(x(inew(l:end)) ,PBGHigh(inew(l:end) ,jnew(l:end)) , ’  b  ’  ) 
hold  on 
clear  inew 
clear  jnew 

counter  =  counter  +1; 

end 

7„  now  print  the  H  field  gap  map 
Hold  on 

while  Hcounter  <  max  ( size  (Hj  )) -1  °/n  make  sure  we  don’t  violate  the... 

max  length  of  the  matrix 
Hnewcounter  =  1; 

while  isequal  (Hj  (Hcounter ),  Hj  (Hcounter +  1)  )  °/0  finds  all  column  ... 
numbers  so  they  can  be  ploted  with  corresponding  x  values 
if  Hcounter  <  max ( s ize ( H j ) ) -1 

Hinew (Hnewcounter )  =  Hi ( Hcount er ) ; 

Hjnew (Hnewcounter )  =  H j ( Hcount er ) ; 

Hcounter  =  Hcounter  +1; 

Hnewcounter  =  Hnewcounter +1 ; 

end 

end 

Hinew (Hnewcounter )  =  Hi ( Hcounter ) ; 

H j new ( Hnewcount er )  =  Hj ( Hcounter ) ; 

plot (x (Hinew (1 : end) )  , HPBGLow (Hinew (1 : end)  , H j  new ( 1 : end))  ,  ’  g  ’  ) 
hold  on 

plot (x (Hinew (1 : end) )  , HPBGHigh (Hinew ( 1 : end )  ,Hjnew(l:end))  ,  ’  k  ’  ) 
hold  on 
clear  Hinew 
clear  Hjnew 

Hcounter  =  Hcounter  +1; 

end 

xlabel( ’  r/a  ’  ,  ’font weight ’  ,  ’bold’  ,  ’ f ontsize ’  ,16) 
ylabel ( ’\omegaa/2\pic ’  ,  ’ f  ontwe ight ’  ,  ’bold’  ,  ’fontsize ’  ,16) 
title({[’Gap  Map  for  triangular  structure’]  ; [’columns  of  \epsilon... 
=  ’  , num2str ( EPSA )  ,  ’  ’,  MatNaml  ,  ’  and  material  \ epsilon= ’  , num2str . .  . 

(EPSB),’  ’,  MatNam2];[’H  polarization  red  and  blue,  E  ... 
polarization  green  and  black’]}) 
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Listing  C.4:  H  field  band  diagram  program.  (appendix3/Efielddispersion.m) 

function  rr  =  Ef ielddispersion (N1 , d , a , EPSA , EPSB , NDC) 

7,  program  for  calculating  the  E-polarized  band  structure  of  the  ... 
triangular 

7,  lattice.  Lines  are  inserted  to  convert  it  to  band  structure  ... 
computation 

7.  for  the  H-polarized  waves.  A  parameter  determining  the  number  .. 
of  plane 
5  7.  waves 

7.N1  =9;  7.  number  of  plane  waves  (i.e.  9)  INPUT 

pi  =  acos (-1.0)  ; 

7othe  lattice  constant  is  scaled  to  unity.  The  raduis  of  the  ... 
column  . 

10  % d  =  1  ;  7,  I  think  this  is  the  spacing  of  columns  INPUT 

7«a  =  0.48;  7»  I  think  this  is  the  radius  of  the  columns  INPUT 

7.  Volume  fraction  of  column  in  the  unit  cell, 
bbeta  =  3 . 14159* a*a/d/d*2 . 0/ sqrt (3 . 0) ; 

7. EPSA  =  1  .  ;  1  column  DIELECTRIC  INPUT 
15  7.EPSB  =  13;7.  material  DIELECTRIC  INPUT 

7.  The  major  symmetry  points  in  the  Brillouin  zone  are  defined  (for 
7.  triangular  geometry) 


20  gxg  =  0.0;  gyg  =  0.0; 

gxx  =  0.0;  gyx  =  2 . 0*pi/sqrt (3 . 0) ; 

gxj  =  2.0*pi/3.0;  gyj  =  2 . 0* pi / sqrt  (3 . 0)  ; 

7.  The  itinerary  around  the  brillouin  zone  is  laid  out. 
gxi(l)  =  gxx; 

25  gyi  (1)  =  gyx  ; 
gxi  ( 2 )  =  gxg; 

gyi  (2)  =  gyg; 

gxi  (3)  =  gxj  ; 
gyi  (3)  =  gyj  ; 

30  gxi (4)  =  gxx  ; 
gyi  (4)  =  gyx; 


35  7.  Construct  the  wave  vectors  used  in  the  plane  wave  calculation 

L  =  0  ; 

gxi  =  2 . 0*pi ; 

gyi  =  -2 . 0* pi / sqrt (3 . 0 ) ; 

40  gx2  =  0.0; 

gy2  =  2 . 0*pi *2 . 0/ sqrt (3 . 0) ; 

7.  generate  all  the  wavevectors  within  a  hexagon  of 

/  diameter  2*nl. 

45  for  ii=-Nl:Nl; 

for  j=-Nl:Nl; 

if  (  (-(Nl  +  1) <(ii  +  j))  k  ((ii  +  j)<(Nl  +  l) ))  , 
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L  =  L  + 1 ; 

gx  (L)  =  gxl * ii  +  gx2 *  j  ; 

gy  (D  =  gyi*ii  +  gy2*j ; 

end  ; 
end  ; 
end  ; 


7„  Construct  the  dielectric  matrix 

*/„  the  volume  of  the  unit  cell  is  defined 

vcel  =  sqrt (3 . 0) /2 . d0*d*d ; 
for  i i  = 1 : L ; 

EPS(ii,ii)  =  bbeta*EPSA  +  ( 1 . 0 - bbeta) *EPSB ; 

°/0  H -method  insert  the  following  line  to  replace  the 
preceding  line . 

1  Construct  the  inverse  dielectric  matrix 
1  EPS ( ii , ii ) =  bbeta/EPSA+(l .0-bbeta)/EPSB ; 


i  1  =  i  i  +  1  ; 
for  j  = i 1 : L ; 

xl  =  gx(ii)  -  gx ( j ) ; 
x2  =  gy(ii)  -  gy (j) ; 
x  =  sqrt (xl*xl+x2*x2)*a/d; 

EPS ( ii  ,  j )  =  2.0*pi*(EPSA-EPSB)*a*a/(vcel)  *  besselj(l, 
x)  /x  ; 

/  H-method  insert  the  following  line  to  replace  the  .. 
preceding 

7,  line,  construct  the  inverse -dielectric  matrix 
V.  EPS(ii  ,j)=2.  dO*pi*(l/EPSA-EPSB)*a*a/  (vcel)*besselj 
(1 , x) /x  ; 

“/,  The  matrix  is  symmetric 

EPS  (  j  ,  i i  )  =  EPS  (  ii  ,  j  )  ; 

end  ; 

end 


°/0  Use  the  itinerary  constructed  above  to  move  around  the 
Brillouin  zone 
n  =  L  ; 

°/0  Counter  for  the  number  of  points  in  the  itinerary 
nknt  =  1  ; 
for  nk  =  1:3; 

nkl  =  nk+1; 
nmax=  20; 
for  1=1 : nmax ; 

dkx  =  (  1 - 1 )*( gxi ( nkl ) -gxi ( nk) )/( nmax )  +  gxi(nk); 

dky  =( 1 -1 )*( gyi (nkl ) -gyi (nk) )/( nmax )  +  gyi(nk); 
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for  ii=l : n ; 


a ( ii  ,  ii )  =  ( ( dkx  +  gx ( ii ) ) 
b(ii,ii)  =  EPS ( ii , ii ) ; 
for  j  =ii  +  l : n ; 
a ( ii  ,  j  ) 


a(j , ii) 
b ( ii  ,  j  ) 
b ( j , ii ) 


end 


'2  + 


0.0; 

0.0; 

EPS ( ii  ,  j  )  ; 
EPS ( ii  ,  j  )  ; 


( dky+gy (ii))"2)/4. 0/pi/pi; 


end 

7,  H-method  insert  the  following  lines  and  replace  the  ... 

preceding  10 
*/.  lines 

7,  for  i  i  =  1  :  n  ; 

7.  a(ii  ,ii)=EPS(ii  ,  ii  )  *  (  (  dkx  +  gx  (ii))~2+(  dky  +  gy  (ii))~2)/4.0/ 

pi/pi ; 

7.  for  j  =  i  i  + 1  :  n  ; 

7.  a  (  ii  ,  j  )  =  EPS  (  ii  ,  j  )*((  dkx  +  gx  (  ii  ))*(  dkx  +  gx  (j  ))  +  (  dky  +  gy  ( 

ii ) ) * .  .  . 

7.  (dky  +  gy  (j  )))  /4 . 0/  pi  /pi  ; 

7.  a ( j  ,  ii )  =a ( ii  ,  j  )  ; 

7,  end 

7,  end 


7.  Find  the  eignevalues 
v  =  eig (a , b)  ; 

7,  Find  the  eignevalues  for  the  H-method  replace  the 
preceding  line 
7.  v=eig  (a)  ; 

7.  Store  the  lowest  10  eigenvalues  for  plotting 
for  li=l : NDC ; 
rr (nknt ,li)=sqrt(v(li)) ; 
kxx (nknt ) =nknt -1 ; 
end 


nknt  =  nknt  +  1  ; 

end 


end 


7.  Plot  lowest  six  eigenvalues 


f igur e  (  1 ) 
hold  on  ; 


plot  (kxx  ,  rr  ( 

,D  , 

’r 

)  ; 

plot  (kxx  ,  rr  ( 

,2)  , 

’  b 

)  ; 

plot  (kxx  ,  rr  ( 

,3)  , 

’g 

)  ; 

plot  (kxx  ,  rr  ( 

,4)  , 

’y 

)  ; 

plot  (kxx  ,  rr  ( 

,5)  , 

’  r 

)  ; 

141 


plot (kxx ,rr(:  ,6)  ,  ’  r  1  )  ; 
hold  of  f  ; 

Listing  C.5:  H  field  band  diagram  program.  (appendix3/Hfielddispersion.m) 

function  ss  =  Hf ielddispersion (N1 , d , a , EPSA , EPSB , NDC) 

7,  program  for  calculating  the  H-polarized  band  structure  of  the  ... 
triangular 

°l  lattice.  A  parameter  determining  the  number  of  plane 

7,  waves 
5 

7.N1  =9;  7.  number  of  plane  waves  (i.e.  9)  INPUT 

pi  =  acos (-1.0)  ; 

7«the  lattice  constant  is  scaled  to  unity.  The  raduis  of  the  ... 
column  . 

7od  =  l;  7.  I  think  this  is  the  spacing  of  columns  INPUT 
10  %a=0.48;  7»  I  think  this  is  the  radius  of  the  columns  INPUT 
7.  Volume  fraction  of  column  in  the  unit  cell, 
bbeta  =  3 . 14159* a*a/d/d*2 . 0/ sqrt (3 . 0) ; 

7«EPSA  =  1  .  ;  7.  column  DIELECTRIC  INPUT 
7«EPSB  =  13  ;  7o  material  DIELECTRIC  INPUT 
15  /  The  major  symmetry  points  in  the  Brillouin  zone  are  defined  (for 
/  triangular  geometry) 


gxg  =  0.0;  gyg  =  0.0; 

20  gxx  =  0.0;  gyx  =  2 . 0*  pi  /  sqrt  (3 . 0)  ; 

gxj  =  2.0*pi/3.0;  gyj  =  2 . 0* pi / sqrt  (3 . 0)  ; 

°l  The  itinerary  around  the  brillouin  zone  is  laid  out. 
gxi(l)  =  gxx; 
gyi(l)  =  gyx; 

25  gxi  (2)  =  gxg ; 

gyi  (2)  =  gyg; 

gxi  (3)  =  gxj  ; 
gyi  (3)  =  gyj  ; 
gxi  (4)  =  gxx  ; 

30  gyi  (4)  =  gyx ; 


7.  Construct  the  wave  vectors  used  in  the  plane  wave  calculation 
35 

L  =  0  ; 

gxi  =  2 . 0*pi ; 

gyi  =  -2 . 0* pi / sqrt (3 . 0 )  ; 

gx2  =  0.0; 

40  gy2  =  2.0*pi*2.0/sqrt(3.0) ; 

7.  generate  all  the  wavevectors  within  a  hexagon  of 

7.  diameter  2*nl. 

for  ii=-Nl : N 1 ; 

45  for  j=-Nl:Nl; 

if  (  (-(Nl+1) <(ii+j ) )  &  ((ii+j)<(Nl+l) )) , 
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L  =  L  + 1 ; 

gx  (L)  =  gxl * ii  +  gx2 *  j  ; 
gy (D  =  gyi*ii  +  gy2*j ; 
end  ; 
end  ; 
end  ; 


7„  Construct  the  dielectric  matrix 

*/„  the  volume  of  the  unit  cell  is  defined 

vcel  =  sqrt (3 . 0) /2 . 0*d*d ; 
for  i i  = 1 : L ; 

EPS (ii , ii) =  bbeta/EPSA+ (1 . 0-bbeta) /EPSB ; 


i  1  =  i  i  +  1  ; 
for  j  = i 1 : L ; 

xl  =  gx(ii)  -  gx ( j ) ; 
x2  =  gy(ii)  -  gy(j) ; 
x  =  sqrt (xl*xl+x2*x2)*a/d; 

EPS (ii , j ) =2 . 0*pi*(l/EPSA-l/EPSB)*a*a/(vcel)*besselj  (1  ,  . . . 
x)  /x  ; 

/  construct  the  inverse -dielectric  matrix 

% 

°/o  The  matrix  is  symmetric 

EPS  (  j  ,  i i  )  =  EPS  (  ii  ,  j  )  ; 

end  ; 

end 


*/„  Use  the  itinerary  constructed  above  to  move  around  the  ... 
Brillouin  zone 
n  =  L  ; 

”/o  Counter  for  the  number  of  points  in  the  itinerary 
nknt  =  1  ; 
for  nk  =  1:3; 

nkl  =  nk+1; 
nmax=  20; 
for  1=1 : nmax ; 

dkx  =(  1 - 1 )*( gxi ( nkl ) -gxi ( nk) )/( nmax )  +  gxi(nk); 

dky  =( 1 -1 )*( gyi (nkl ) -gyi (nk) )/( nmax )  +  gyi(nk); 

for  ii=l : n; 

a(ii ,ii)=EPS(ii , ii ) * ( ( dkx+gx (ii))~2+( dky+gy (ii))~2)/4.0/... 
pi/pf ; 

for  j  =  ii  +  1 : n ; 

a ( ii , j ) =  EPS ( ii , j ) * ( ( dkx  +  gx ( ii ) ) * ( dkx  +  gx ( j ) ) + ( dky  +  gy . .  . 
(ii))*.  .  . 
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end 


(dky+gy (j)))/4. 0/pi/pi; 
a ( j  ,  ii)=a(ii , j) ; 


end 

7, 


7o  Find  the  eignevalues 
v=eig ( a) ; 

7, 

7.  Store  the  lowest  10  eigenvalues  for  plotting 
for  li =1 : NDC ; 
ss (nknt ,li)=sqrt(v(li)) ; 
kxx (nknt ) =nknt -1 ; 
end 

nknt  =  nknt  +  1  ; 

end 


end 


7.  Plot  lowest  six  eigenvalues 


f igur e  (  1 ) 
hold  on  ; 


plot  (kxx  ,  ss  ( 

,D  , 

’r 

) 

plot  (kxx  ,  ss  ( 

,2)  , 

’  b 

) 

plot  (kxx  ,  ss  ( 

,3)  , 

’g 

) 

plot  (kxx  ,  ss  ( 

,4)  , 

’y 

) 

plot  (kxx  ,  ss  ( 

,5)  , 

’  r 

) 

plot  (kxx  ,  ss  ( 

,6)  , 

’  r 

) 

hold  off  ; 
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Appendix  D.  Index  Program  for  AlGaAs 


This  appendix  shows  the  Matlab®  program  to  calculate  the  indices  of  AlGaAs 
courtesy  of  Lt.Col.  Lott. 


D.l  The  MATLAB®  Index  Program 

The  following  code  returns  the  real  and  imaginary  (if  there  is  any  absorption  at 
that  wavelength  and  mole  fraction)  index  of  AlGaAs.  This  code  was  based  off  of  the 
C-code  created  by  Lt.Col.  Lott  and  was  converted  to  Matlab®  by  Capt.  Ochoa. 
This  code  is  a  function  and  can  easily  be  used  from  the  Matlab®  command  line. 
The  inputs  are  the  mole  fraction  of  Aluminum  and  the  wavelength  in  angstroms. 

Listing  D.l:  AlGaAs  index  calculation  program.  (appendix4/AlGaAs2Indxr.m) 

function  IndexOf Ref ract ion  =  AlGaAs2Indxr (x ,  lambda) 

7.  AlGaAs2Indxr  Index  of  Refraction  of  AlGaAs  between  1.5  to  2.2  eV 
°/0  IndexOf  Ref  raction  =  AlGaAs2Indxr  (x  ,  lambda); 

7.  x  -  mole  fraction  of  A1 
%  lambda  -  Angstroms 

7. 

7.  Calculates  the  real  and  complex  index  of  refraction  for  AlxGal-... 
xAs 

7.  between  1.5  to  2.2  eV  .  Ref:  Aspnes  et  al  JAP  60(2)  15  July  1986 

7.  8265.3  Ang  to  5635  Ang 

/  based  on  C-code  by  Lt  Col  James  A.  Lott 

/  note:  does  not  call  AlGaAslIndxr  for  x>=.8...  may  need  to  modify... 
later . . . 

7.  Converted  to  Matlab  by:  Lt  Bernhard,  Capt  Ochoa  1/21/00 

7o  initialize  B 

B  =  0; 

if  x  <  0 
x  =  0 ; 

end 

if  x  >  1 
x  =  1  ; 

end 

dumvar  =  x*10; 

dumint  =  f loor ( dumvar ) ; 

7.  convert  lambda  to  energy  in  eV 
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hbo  =  1 . 2398/ ( lambda*10~ -4) ; 

for  n  =  1:2 

if  dumint  ==  0 

if  lambda  >  6440 

A  =  5 . 4323- . 00036682* lambda+1 . 85 1 5E-8* lambda ~2 ; 
else 

A  =  12. 117-. 0025398* lambda  + 1 . 945E-7* lambda  ~ 2 ; 

end 

B  =  -.40018+. 31110* hbo-. 0023811* hbo~2; 

end 

if  dumint  ==  1 

A  =  7. 1961-. 00087025* lambda +  5 . 279E -8* lambda  ~ 2 ; 
if  hbo  <  1.5 
B  =  0  ; 
else 

B  =  -  .  4627+ . 31595*hbo ; 

end 

end 

if  dumint  ==  2 

if  lambda  >  7300 

A  =  8 . 954- . 00121 13* lambda+6 . 6093E-8* lambda'2 ; 
else 

A  =  7. 9718-. 0011754* lambda  +  7 . 97 35 E -8* lambda  ~ 2 ; 

end 

if  hbo  <  1.65 
B  =  0; 
else 

B  = -  .  59874+ .36371* hbo; 

end 

end 

if  dumint  ==  3 

A  =  6 . 2665- . 00061784* lambda +  3 . 3085E-8* lambda  ~ 2  ; 
if  hbo  <  1.71 
B  =  0; 
else 

B  =  -  .723  24+ .42  533* hbo ; 

end 

end 

if  dumint  ==  4 

A  =  6. 3718-. 0006667* lambda +  3 . 6448E-8* lambda  ~2 ; 
if  hbo  <  1.82 
B  =  0; 
else 

B=-.67961+. 37484* hbo ; 

end 

end 
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if  dumint  ==  5 

A  =  6. 6566-. 000777* lamb da +  4 . 45688E-8* lambda  ~2 ; 
if  hbo  <  1.99 
B  =  0; 
else 

B  =  -  1 . 3012+ . 655*  hbo ; 

end 

end 


if  dumint  ==  6 

A  =  6 . 82823- . 0008582* lambda  +  5 . 1 38948E -8* lambda  "  2  ; 
if  hbo  <  2.1 
B  =  0; 
else 

B  =  -  1 .213+ . 58*hbo ; 

end 

end 


if  dumint  ==  7 

A  =  7. 1055-. 0009804* lambda +  6 . 10584E-8* lambda"  2 ; 
B  =  0; 

end 


if  dumint  ==  8 

A  =  6 . 52743- . 00086692* lambda  +  5 . 5396E-8* lambda  “2  ; 

B  =  0; 

end 

if  dumint  ==  9 

A  =  5. 1113-0.00047623*  lambda  +  2 . 7779E-08*  lambda  ~  2  ;  7,5635  to... 

8265 
B  =  0; 

end 


if  dumint  ==  10 

A  =  4.7921-0.00041411*  lambda  +  2 . 405E-08*  lambda  ~  2  ;  7,5635  to... 

8265 
B  =  0; 

end 

if  B  <  0 
B  =  0; 

end 

A1  (n)  =  A  ; 

B1  (n)  =  B  ; 


147 


dumint  =  dumint  +  1 ; 

135 

end 

rindex  =  A 1 ( 1 )  -  ( A 1 ( 1 )  -  Al(2))*(10*x  - 

cindex  =  ( B 1  ( 1 )  -  ( B 1  (  1 )  -  Bl(2))*(10*x 

140 

IndexOf Ref raction  =  rindex  -  i*cindex; 


floor  (10*x) )  ; 

-  f loor  (10*x) ) ) 
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Appendix  E.  Holographic  Interference  Program 


This  appendix  explains  the  equations  used  in  the  program  that  predicts  the  holo¬ 
graphic  interference  pattern  used  in  the  holographic  polymerization/lithography 
process. 

E.l  Geometry  and  Equations  of  the  Interference  Pattern  of  Two-Dimensional 
PhC  setup 

The  following  text  is  courtesy  of  Dr.  Sutherland.  The  text  explains  the  equa¬ 
tions  that  were  typed  into  Mathcad®  which  is  displayed  after  this  text  explanation. 

The  following  sketch  illustrates  the  geometry  and  equations  that  follow. 


Z 


Figure  E.l:  The  geometry  that  explains  the  equations  used  in  the 

holographic  polymerization/lithography  process. 


The  sample  is  in  the  x-y  plane,  and  the  z-axis  is  perpendicular  to  the  sample. 
The  wave  vectors  of  figure  E.l  are  given  by: 
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ki  =  —k(sm(9)x  +  cos  (9)z) 

k2  =  —k((sm(9)  cos((J)2))x  +  (sin(0)  sm((j)2))y  +  cos (9)z) 
k3  =  —  k((sm(9)  cos(03))£  +  (sin(0)  sin(03))y  +  cos (9)z) 


Where  k  is  the  wave  vector  (2n)/\  and  A  is  the  wavelength  of  the  laser  that  is 
writing  the  structure.  The  fields  are  given  by: 


Ei=y 

E2  =  —  cos (9)  cos {4>2)x  —  cos (9)  sin (cj>2)y  +  sin (9)z 
E3  =  V2(  —  (sin(a)  sin(03)  +  cos(a)  cos (9)  cos(cj)3))x  +  (sin(a)  cos(03) 
—  cos(a)  008(0)  sin (<f)3))y  +  cos(a)  sin(0)i) 


The  angel  a  is  the  angle  of  polarization  in  the  local  plane  of  polarization  of  beam  3; 
in  other  words  this  is  the  the  polarization  of  beam  1  that  has  an  added  angle  a  with 
respect  to  the  plane  of  incidence  in  the  z-plane.  The  intensity  in  the  sample  is  given 
by: 

I(r)  =  4+  Smnexp(iGmn)  (E.3) 

m^n 

where  Gmn  =  km  —  kn  are  the  reciprocal  lattice  vectors,  and  Smn  =  Em  ■  En  are  the 
structure  factors  and  m,n  =  1...3. 


It  is  worth  noting  that  the  9  values  govern  the  overall  pattern,  (f)  governs  the 
overall  periodicity  and  a  governs  the  intensity  of  each  pattern  or  the  volume  fraction. 

E.2  Holographic  Interference  Equations  in  MATHCAD® 

This  program  will  calculate  a  two-dimensional  interference  pattern  achieved 
when  three  laser  beams  intersect.  It  depends  on  the  angle  theta  (incident  angle)  the 
rotation  angle  (phi)  from  a  fixed  reference  plane  x-z  and  the  polarization  of  each 
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beam  in  accordance  to  each  beam’s  plane  of  incidence  (which  is  related  to  the  angle 
phi  also) .  These  calculations  are  oriented  with  the  sample  that  is  being  written  on  to 
be  in  the  x-y  plane.  The  angle  alpha  is  the  angle  of  polarization  with  respect  to  the 
plane  of  incidence  of  beam  3.  Theta,  alpha  and  phi  must  be  in  radians 
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92  := 

93  ■= 
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02  :=  120  •  ^ 

03  :=  240  •  ^ 


a  .  10  •  180 


A  :=  0.3648 


k  :  =  2-f 


The  wave  vectors  are  as  follows: 
(^  \ 


x  :  = 


V  '■= 


1 
0 

v0/ 

/  o  \ 

i 

V0/ 
/  o  N 


z 

0 

V1 

kl 

:=  —k 

k2 

:=  -k 

k3 

:=  -k 

The  E-helds  are  given  by  the  following  equations: 


cos (92)  ■  z] 
cos(93)  ■  z] 
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El  :=  y 


E 2  :=  —  (cos(9 2)  •  cos((p 2))  •  a;  —  (cos(6,2)  •  sin(4> 2))  •  y  +  sin{92 )  •  z 

E3  :=  \/2-  [-(sro(a)  -sin{<j)3)  +  cos(a)  ■ cos{93 )  -cos((j> 3))  ■i+(sm(«)  -cos((j) 3)  — 
cos(a)  •  cos(03)  •  sin(<t> 3))  ■  y  +  (cos(a)  •  sin(93 ))  •  z] 


Now  calculate  the  structure  factors  (Smn  =  Em  ■  En 
(  El  ■  E2  ^ 


S'  :  = 


V 


El  ■  E3 
E2  •  El 
E2  •  E3 
E3  •  El 
E3  •  E2 

Now  calculate  the  reciprocal  lattice  vectors  (Gmn  =  km  —  kn). 
(  kl  -  k2  ^ 

/cl  —  /c3 


/ 


G  := 


k2  —  kl 
k2  -  k3 
k3  —  kl 
k3  -  fc2 


/ 


r(a,b,c )  := 


/a\ 

6 

vc/ 


n  :=  0,1. .5 

/(a,  6,  c)  :=  4  +  E,(^-e!'G"'r(a’M) 


IV  :=  100 
scale  :=  0.01 


This  gives  a  scale  such  that  30  is  0.3  microns 
q  :=  1,2. .N 
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w  :=  1,  2..N 
aq  :=  q  ■  scale 
bw  :=  w  ■  scale 
cw  :=  w  ■  scale 

•  I(ciq,  bw,  0) 

Then  insert  a  plot  and  use  ’JO’  as  the  argument  to  have  plotted. 
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Appendix  F.  SU-8  Processing  Steps 


This  appendix  explains  the  process  used  to  spin  on  and  develop  the  SU-8  photo¬ 
resist. 

F.l  Process  for  the  Spin  on  and  Development  of  SU-8  Photo-Resist. 

The  process  used  in  this  experiment  to  spin-on  the  SU-8  2010  negative  photo¬ 
resist  by  MicroChem  for  a  15  pm  thick  layer  is  detailed  below: 

Spin-on 

The  glass  substrate  was  cleaned  with  Acetone  and  Methonal  and  dried  with 
Air/Nitrogen. 

Then  it  was  baked  at  200°  C  for  5  minutes. 

Then  it  was  placed  in  a  spinner  that  was  set  for  the  following  parameters: 

Ramp  to  500  rpm  at  100  rpm/sec. 

Ramp  to  3000  rpm  at  300  rpm/sec  and  hold  for  30  seconds. 

Soft  bake  to  evaporate  the  solvent  and  density  the  film  at  65°  C  for  1  minute,  then  at 
95°  C  for  2  minutes. 

Post-exposure  process 

After  exposure,  a  post-exposure  bake  is  necessary  for  the  cross-linking  to  occur. 
Heat  on  hot  plate  at  65°  C  for  1  minute,  then  at  95°  C  for  2  minutes.  It  is  important 
to  note  that  since  the  SU-8  is  on  a  glass  substrate,  the  hot  plate  may  need  to  be 
turned  up  higher  to  ensure  the  temperature  of  the  SU-8  reaches  95°  C.  It  must  get 
above  95°  C  for  the  cross-linking  to  occur;  however,  it  can  not  get  hotter  than  115° 
to  120°  C.  To  facilitate  this  requirement,  the  samples  were  put  on  a  105°  C  hot  plate 
for  the  2  minute  time  period.  This  ensures  the  above  restrictions  are  meet.  Allow  the 
samples  to  cool  slowly.  Put  in  SU-8  developer  from  MicroChem,  for  3  minutes  with 
slight  agitation.  The  samples  were  then  rinsed  with  isopropyl  alcohol  and  allowed  to 
dry.  Once  dry,  hard  bake  at  180°  C  for  1  minute  in  order  to  ensure  the  hardness  of 
the  structure. 
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NANO™  SU-8  2000 

Negative  Tone  Photoresist 
Formulations  2002-2025 

a  High  aspect  ratio  imaging 
Near  vertical  side  walls 

m  Near  UV  (350-400nm)  processing 

^5  Improved  coating  properties 

Uniformity  (lower  surface  tension) 
Adhesion 

3  Faster  drying 

Improved  throughput 

SU-8  2000  is  a  high  contrast,  epoxy  based  photoresist  de¬ 
signed  for  micromachining  and  other  microelectronic  ap¬ 
plications,  where  a  thick,  chemically  and  thermally  stable 
image  is  desired.  SU-8  2000  is  an  improved  formulation 
of  SU-8,  which  has  been  widely  used  by  MEMS  produc¬ 
ers  for  many  years.  By  using  a  faster  drying,  more  polar 
solvent  system,  improved  coating  properties  and  higher 
throughput  are  realized.  Film  thicknesses  of  0.5  to 
>200mm  can  be  achieved  with  a  single  coat  process.  The 
excellent  imaging  characteristics  of  SU-8  are  maintained. 
The  exposed  and  subsequently  cross-linked  portions  of 
the  film  are  rendered  insoluble  to  liquid  developers.  SU-8 
2000  has  very  high  optical  transparency  above  360nm, 
which  makes  it  ideally  suited  for  imaging  near  vertical 
sidewalls  in  very  thick  films.  SU-8  2000  is  best  suited  for 
permanent  applications  where  it  is  imaged,  cured  and  left 
in  place 

Process  Guidelines 

SU-8  2000  is  most  commonly  processed  with  conventional 
near  UV  (350-400nm)  radiation,  although  it  may  be  im¬ 
aged  with  e-beam  or  x-ray.  i-line  (365nm)  is  recommended. 
Upon  exposure,  cross-linking  proceeds  in  two  steps  (1)  for¬ 
mation  of  a  strong  acid  during  the  exposure  process,  fol¬ 
lowed  by  (2)  acid-initiated,  thermally  driven  epoxy  cross- 
linking  during  the  post  exposure  bake  (PEB)  step. 


5pm,  10pm  and  20pm  post  arrays  in  a  50pm  thick  film. 


Figure  F.l:  Data  sheet  for  SU-8. 
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A  normal  process  is:  spin  coat,  soft  bake,  expose,  post  ex¬ 
pose  bake  (PEB)  followed  by  develop.  A  controlled  hard 
bake  is  recommended  to  further  cross-link  the  imaged  SU-8 
2000  structures  when  they  will  remain  as  part  of  the  device. 
The  entire  process  should  be  optimized  for  the  specific  ap¬ 
plication.  A  baseline  process  is  given  here  to  be  used  as  a 
starting  point. 

Substrate  Pretreat 

To  obtain  maximum  process  reliability,  substrates  should  be 
clean  and  dry  prior  to  applying  the  SU-8  2000  resist.  Start 
with  a  solvent  cleaning,  or  a  rinse  with  dilute  acid,  followed 
by  a  DI  water  rinse.  Where  applicable,  substrates  should  be 
subjected  to  a  piranha  etch  /  clean  (H^SO^  &  F^O^).  To 
dehydrate  the  surface,  bake  at  200°C  for  5  minutes  on  a  con¬ 
tact  hot  plate  or  30  minutes  in  a  convection  oven.  Adhesion 
promoters  are  typically  not  required.  For  applications  that 
require  electroplating  and  subsequent  removal  of  SU-8  2000 
apply  MicroChem’s  OmniCoat  prior  to  processing. 

Coat 

SU-8  2000  resists  are  designed  to  produce  low  defect  coat¬ 
ings  over  a  very  broad  range  of  film  thickness.  The  film 
thickness  versus  spin  speed  data  displayed  in  Table  1  and 
Figure  1  provide  the  information  required  to  select  the  ap¬ 
propriate  SU-8  2000  resist  and  spin  conditions,  to  achieve 
the  desired  film  thickness. 

The  recommended  coating  conditions  are: 

(1)  STATIC  Dispense:  Approximately  1ml  of  SU-8  2000 
per  inch  of  substrate  diameter. 

(2)  Spread  Cycle:  Ramp  to  500  rpm  at  100  rpm/second 
acceleration.  This  will  take  5  seconds. 

(3)  Spin  Cycle:  Ramp  to  final  spin  speed  at  an  acceleration 
of  300  rpm/second  and  hold  for  a  total  of  30  seconds. 


Product  Name 

Viscosity** 

(cSt) 

Thickness** 

(/i  ms) 

Spin  Speed 
(rpm) 

2 

3000 

SU-8  2002 

7.5 

2.5 

2000 

3 

1000 

5 

3000 

SU-8  2005 

45 

6 

2000 

7.5 

1000 

7 

3000 

SU-8  2007 

140 

8.5 

2000 

12.5 

1000 

10 

3000 

SU-8  2010 

380 

13 

2000 

20 

1000 

15 

3000 

SU-8  2015 

1250 

21 

2000 

38 

1000 

25 

3000 

SU-8  2025 

4500 

41 

2000 

75 

1000 

Table  1.  t  hickness  vs.  spin  speed  data  for  selected  SU-8 
2000  resists. 

**  Approximate 


SU8-2000  Spin  Speed  Curve 


Figure  1.  Spin  speed  vs.  thickness  curves  for  selected 
SU-8  2000  resists. 

Soft  Bake 

After  the  resist  has  been  applied  to  the  substrate,  it  must  be 
soft  baked  to  evaporate  the  solvent  and  densify  the  film. 
SU-8  2000  is  normally  baked  on  a  level  hot  plate,  although 
convection  ovens  may  be  used.  The  following  bake  times 
are  based  on  contact  hot  plate  processes.  Bake  times  should 
be  optimized  for  proximity  and  convection  oven  bake  pro¬ 
cesses  since  solvent  evaporation  rate  is  influenced  by  the 
rate  of  heat  transfer  and  ventilation. 

For  best  results,  ramping  or  stepping  the  soft  bake  tempera¬ 
ture  is  recommended.  Lower  initial  bake  temperatures  al¬ 
low  the  solvent  to  evaporate  out  of  the  film  at  a  more  con¬ 
trolled  rate,  which  results  in  better  coating  fidelity,  reduced 
edge  bead  and  better  resist  -to-subslrate  adhesion.  Refer  to 
Table  2.  for  TWO  STEP  contact  hot  plate  process  recom¬ 
mendations. 


Product  Name 

Thickness 

(fims) 

Prc-hakc 
(fib,  65°  C 

Softbakc 
(fb,  95°  C 

2 

i 

2 

SU-8  2002 

2.5 

i 

2 

3 

1 

2 

5 

1 

2 

SU-8  2005 

6 

1 

2 

7.5 

1 

2 

7 

i 

2 

SU-8  2007 

8.5 

i 

2 

12.5 

1 

2 

10 

1 

2 

SU-8  2010 

13 

1 

2 

20 

i 

3 

15 

1 

2 

SU-8  2015 

21 

i 

3 

38 

2 

5 

25 

1 

3 

SU-8  2025 

41 

2 

5 

75 

3 

9 

Table  2.  Recommended  soft  bake  parameters 


Figure  F.2:  Data  sheet  for  SU-8. 
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Expose 

SU-8  is  optimized  lor  near  UV  (350-400nm)  exposure,  i- 
line  exposure  tools  are  recommended.  SU-8  is  virtually 
transparent  and  insensitive  above  400nm  but  has  high  ac¬ 
tinic  absorption  below  350nm.  This  can  be  seen  in  Figure 
2.  Excessive  dose  below  350nm  may,  therefore,  result  in 
over  exposure  of  the  top  portion  of  the  resist  film,  resulting 
in  exaggerated  negative  sidewall  profiles  or  T-topping.  The 
optimal  exposure  dose  will  depend  on  film  thickness  (thicker 
films  require  higher  dosage)  and  process  parameters.  The 
exposure  dose  recommendations  in  Table  3.  are  based  on 
source  intensity  measurements  taken  with  an  i-line  (365nm) 
radiometer  and  probe. 

Expose  tip:  When  using  a  broad  spectral  output  source,  for 
best  imaging  results,  i.e.  straiglitest  sidewalls,  filter  out 
excessive  energy  below  350nm. 

Catastrophic  adhesion  failure,  severely  negative  sidewalls 
and  excessive  cracking  often  indicate  an  under  cross-link¬ 
ing  condition.  To  correct  the  problem,  increase  the  expo¬ 
sure  dose  and/or  increase  the  post  exposure  bake  (PEB)  time. 


Wavelength  (nm) 

Figure  2.  SU-8  absorbance  vs.  film  thickness 


SU-8  2000  Exposure  vs  Film  Thickness 


Film  Thickness  (um) 


Table  3.  Recommended  expose  dose  processes 

Post  Expose  Bake 

Following  exposure,  a  post  expose  bake  (PEB)  must  be  per¬ 
formed  to  selectively  cross-link  the  exposed  portions  of  the 
film.  This  bake  can  be  performed  either  on  a  hot  plate  or  in 
a  convection  oven.  Optimum  cross-link  density  is  obtained 
through  careful  adjustments  of  the  exposure  and  PEB  pro¬ 
cess  conditions.  The  bake  recommendations  below  are  based 
on  results  obtained  with  a  contact  hot  plate. 

PEB  tip:  SU-8  is  readily  cross-linked  and  can  result  in  a 
highly  stressed  film.  To  minimize  stress,  wafer  bowing  and 
resist  cracking,  a  slow  ramp  or  TWO  STEP  contact  hot  plate 
process,  as  shown  in  Table  4.,  is  recommended.  Rapid  cool¬ 
ing  after  PEB  should  be  avoided. 


Product  N ame 

Thickness 
(a  ms) 

PEB  1 

(if:  65"  C 

PEB  2 
®,95°C 

2 

1 

1 

SU-8  2002 

2.5 

1 

1 

3 

1 

1 

5 

1 

1 

SU-8  2005 

6 

1 

1 

7.5 

1 

1 

7 

1 

1 

SU-8  2007 

8.5 

1 

2 

12.5 

1 

2 

10 

1 

2 

SU-8  2010 

13 

1 

2 

20 

1 

2 

15 

1 

2 

SU-8  2015 

21 

1 

2 

38 

1 

3 

25 

1 

3 

SU-8  2025 

41 

1 

3 

75 

1 

7 

Table  4.  Recommended  post  expose  bake  parameters 


Figure  F.3:  Data  sheet  for  SU-8. 


157 


Develop 

SU-8  2000  resists  have  been  optimized  for  use  with 
MicroChem's  SU-8  Developer.  Immersion,  spray  or  spray- 
puddle  processes  can  be  used.  Other  solvent  based  develop¬ 
ers  such  as  ethyl  lactate  and  diacetone  alcohol  may  also  be 
used.  Strong  agitation  is  recommended  for  high  aspect  ra¬ 
tio  and/or  thick  film  structures.  Recommended  develop  times 
are  given  in  Table  5.  for  immersion  processes.  These  pro¬ 
posed  develop  times  are  approximate,  since  actual  dissolu¬ 
tion  rates  can  vary  widely  as  a  function  of  agitation  rate, 
temperature  and  resist  processing  parameters. 


Product  Name 

Thickness 

(urns) 

Development 

('minutes) 

2 

i 

SU-8  2002 

2.5 

i 

5 

i 

5 

i 

SU-8  2005 

6 

i 

7.5 

1 

7 

i 

SU-8  2007 

8.5 

2 

12.5 

3 

10 

2 

SU-8  2010 

13 

3 

20 

3 

15 

3 

SU-8  2015 

21 

3 

38 

4 

25 

4 

SU-8  2025 

41 

5 

75 

7 

Table  5.  Recommended  develop  processes 


Rinse  and  Dry 

Following  development,  the  substrate  should  be  rinsed  briefly 
with  isopropyl  alcohol  (IPA),  then  dried  with  a  gentle  stream 
of  air  or  nitrogen. 

Rinse  tip:  If  a  white  film  is  produced  during  rinse,  this  is 
an  indication  that  the  substrate  has  been  under  developed. 
Simply  immerse  or  spray  the  substrate  with  SU-8  developer 
to  remove  the  film  and  complete  the  development  process. 
Repeat  the  rinse  step. 

Hard  Bake  (cure) 

SU-8  2000  has  good  mechanical  properties,  therefore  hard 
bakes  are  normally  not  required.  For  applications  where  the 
imaged  resist  is  to  be  left  as  part  of  the  final  device,  the 
resist  may  be  ramp/step  hard  baked  between  150-200°C  on 
a  hot  plate  or  in  a  convection  oven  to  further  cross  link  the 
material.  Bake  times  vary  based  on  type  of  bake  process  and 
film  thickness. 

Removal 

SU-8  2000,  after  expose  and  PEB,  is  a  highly  cross-linked 
epoxy,  which  makes  it  extremely  difficult  to  remove  with 


conventional  solvent  based  resist  strippers.  MicroChem’s 
Remover  PG  will  swell  and  lift  off  minimally  cross-linked 
SU-8  2000.  However,  if  OmniCoat  has  been  applied  immer¬ 
sion  in  Remover  PG  should  effect  a  clean  and  thorough  Lift- 
Off  of  the  SU-8  2000  Material .  It  will  not  remove  fully  cured  or 
hard  baked  SU-8  2000  without  the  use  of  OmniCoat.  Alter¬ 
nate  removal  processes  include  immersion  in  oxidizing  acid 
solutions  such  as  piranha  etch  /  clean,  plasma  ash,  RIE,  laser 
ablation  and  pyrolosis. 

To  remove  minimally  cross-linked  SU-8  2000,  or  if  using 
Omnicoat,  with  Remover  PG,  heat  the  bath  to  50-80°C  and 
immerse  the  substrates  for  30-90  minutes.  Actual  strip  time 
will  depend  on  resist  thickness  and  cross-link  density 
For  more  information  on  MicroChem  Omnicoat  and  Remover 
PG  please  see  the  relevant  product  data  sheets. 

Storage 

Store  SU-8  2000  resists  upright  in  tightly  closed  containers  in 
a  cool,  dry  environment  away  from  direct  sunlight  at  a  tem¬ 
perature  of  40-70°F(4-21°C).  Store  away  from  light,  acids, 
heat  and  sources  of  ignition.  Shelf  life  is  twelve  months  from 
date  of  manufacture. 

Disposal 

SU-8  2000  resists  may  be  included  with  other  waste  con¬ 
taining  similar  organic  solvents  to  be  discarded  for  destruc¬ 
tion  or  reclaim  in  accordance  with  local  state  and  federal 
regulations.  It  is  the  responsibility  of  the  customer  to  ensure 
the  disposal  of  SU-8  2000  resists  and  residues  made  in  ob¬ 
servance  all  federal,  state,  and  local  environmental  regula¬ 
tions. 

Environmental,  Health  and  Safety 

Consult  product  Material  Safety  Data  Sheet  before  working 
with  SU-8  2000  resists.  Handle  with  care.  Wear  chemical 
goggles,  chemical  gloves  and  suitable  protective  clothing 
when  handling  SU-8  2000  resists.  Do  not  get  into  eyes,  or 
onto  skin  or  clothing.  Use  with  adequate  ventilation  to  avoid 
breathing  vapors  or  mist.  In  case  of  contact  with  skin,  wash 
affected  area  with  soap  and  water.  In  case  of  contact  with 
eyes,  rinse  immediately  with  water  and  flush  for  15  minutes 
lifting  eyelids  frequently.  Get  emergency  medical  assistance. 

The  information  is  based  on  our  experience  and  is,  we  be¬ 
lieve  to  be  reliable,  but  may  not  be  complete.  We  make  no 
guarantee  or  warranty,  expressed  or  implied,  regarding  the 
information,  use,  handling,  storage,  or  possession  of  these 
products,  or  the  application  of  any  process  described  herein 
or  the  results  desired,  since  the  conditions  of  use  and  han¬ 
dling  of  these  products  are  beyond  our  control. 


1254  Chestnut  Street 
Newton,  MA  02464 


tel:  (617)965-5511  fax:  (617)965-5818 
email:  mcc@micttichem.com  www.microchem.com 


Figure  F.4:  Data  sheet  for  SU-8. 
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triangular  PhC  structures.  The  fabrication  techniques  and  materials  used  to  create  the  PhC  in  part  determine  the  characterization 
technique  required  to  investigate  the  PBG.  Techniques  used  here  include  the  coupling  of  a  beam  by  means  of  a  prism  into  a  wave-guiding 
medium  and  edge  firing.  Two  commercial  FDTD  modeling  programs  were  evaluated.  Finally,  FDTD  modeling  and  simulation  results  are 
compared  to  the  design  and  experimental  results.  Experimental  techniques  and  results,  along  with  modeling  and  simulation  results,  will  be 
discussed. 
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